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PREFACE 


In the first year of calculus we study limits, derivatives, and integrals of functions 
with a single input, and a single output. The transition to advanced calculus is 
made when we generalize the notion of “function” to something which may have 
multiple inputs and multiple outputs. In this more general context limits, derivatives, 
and integrals take on new meanings and have new geometric interpretations. For 
example, in first-year calculus the derivative represents the slope of a tangent line at 
a specified point. When dealing with functions of multiple variables there may be 
many tangent lines at a point, so there will be many possible ways to differentiate. 

The emphasis of this book is on developing enough familiarity with the material 
to solve difficult problems. Rigorous proofs are kept to a minimum. I have included 
numerous detailed examples so that you may see how the concepts really work. All 
exercises have detailed solutions that you can find at the end of the book. I regard 
these exercises, along with their solutions, to be an integral part of the material. 

The present work is suitable for use as a stand-alone text, or as a companion 
to any standard book on the topic. This material is usually covered as part of a 
standard calculus sequence, coming just after the first full year. Names of college 
classes that cover this material vary greatly. Possibilities include advanced calculus, 
multivariable calculus, and vector calculus. At schools with semesters the class may 
be called Calculus III. At quarter schools it may be Calculus IV. 

The best way to use this book is to read the material in each section and then try 
the exercises. If there is any exercise you don’t get, make sure you study the solution 
carefully. At the end of each chapter you will find a quiz to test your understanding. 
These short quizzes are written to be similar to one that you may encounter in a 
classroom, and are intended to take 20—30 minutes. They are not meant to test every 
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idea presented in the chapter. The best way to use them is to study the chapter until 
you feel confident that you can handle anything that may be asked, and then try the 
quiz. You should have a good idea of how you did on it after looking at the answers. 
At the end of the text there is a final exam similar to one which you would find at 
the conclusion of a college class. It should take about two hours to complete. Use it 
as you do the quizzes. Study all of the material in the book until you feel confident, 
and then try it. 

Advanced calculus is an exciting subject that opens up a world of mathematics. 
It is the gateway to linear algebra and differential equations, as well as more 
advanced mathematical subjects like analysis, differential geometry, and topology. 
Itis essential for an understanding of physics, lying at the heart of electro-magnetics, 
fluid flow, and relativity. It is constantly finding new use in other fields of science 
and engineering. I hope that the exciting nature of this material is conveyed here. 
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CHAPTER 1 





Functions of 
Multiple Variables 


1.1 Functions 


The most common mental model of a function is a machine. When you put some 
input in to the machine, you will always get the same output. Most of first year 
calculus dealt with functions where the input was a single real number and the output 
was a single real number. The study of advanced calculus begins by modifying this 
idea. For example, suppose your “function machine” took two real numbers as its 
input, and returned a single real output? We illustrate this idea with an example. 


EXAMPLE 1-1 
Consider the function 


foy =x +y 
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For each value of x and y there is one value of f(x, y). For example, if x = 2 and 
y = 3 then 
f (2,3) =2° +37 =13 


One can construct a table of input and output values for f (x, y) as follows: 





x|y[fG.») 
ojo 0 
ijol 1 
oji I 
1|1| 2 
112] 5 
2|1| 5 











Problem 1 Evaluate the function at the indicated point. 


1. f(x,y) =x? + y?; Œ, y) = (3, 2) 
2. g(x, y) = sin x + cos y; (x, y) = (0, 2) 
3. h(x, y) =x’ sin y; (x, y) = (2, 2) 


Unfortunately, plugging in random points does not give much enlightenment as 
to the behavior of a function. Perhaps a more visual model would help.... 


1.2 Three Dimensions 


In the previous section we saw that plugging random points in to a function of two 
variables gave almost no enlightening information about the function itself. A far 
superior way to get a handle on a particular function is to picture its graph. We'll 
get to this in the next section. First, we have to say a few words about where such 
a graph exists. 

Recall the steps required to graph a function of a single variable, like g (x) — 3x. 
First, you set the function equal to a new variable, y. Then you plot all the points 
(x, y) where the equation y — g(x) is true. So, for example, you would not plot 
(0, 2) because 0 Z 3 - 2. But you would plot (2, 6) because 6 = 3 - 2. 

The same steps are required to plot a function of two variables, like f(x, y). 
First, you set the function equal to a new variable, z. Then you plot all of the points 
(x, y, z) where the function z — f (x, y) is true. So we are forced to discuss what 
it means to plot a point with three coordinates, like (x, y, z). 
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AZ 








Figure 1-1 Three mutually perpendicular axes, drawn in perspective 


Coordinate systems will play a crucial role in this book, so although most readers 
will have seen this, it is worth spending some time here. To plot a point with two 
coordinates such as (x, y) = (2, 3) the first step is to draw two perpendicular axes 
and label them x and y. Then locate a point 2 units from the origin on the x-axis 
and draw a vertical line. Next, locate a point 3 units from the origin on the y-axis 
and draw a horizontal line. Finally, the point (2, 3) is at the intersection of the two 
lines you have drawn. 

To plot a point with three coordinates the steps are just a bit more complicated. 
Let’s plot the point (x, y, z) = (2, 3, 2). First, draw three mutually perpendicular 
axes. You will immediately notice that this is impossible to do on a sheet of paper. 
The best you can do is two perpendicular axes, and a third at some angle to the 
other two (see Figure 1-1). With practice you will start to see this third axis as a 
perspective rendition of a line coming out of the page. When viewed this way it 
will seem like it is perpendicular. 

Notice the way in which we labeled the axes in Figure 1-1. This is a convention, 
i.e., something that mathematicians have just agreed to always do. The way to 
remember it is by the right hand rule. What you want is to be able to position 
your right hand so that your thumb is pointing along the z-axis and your other 
fingers sweep from the x-axis to the y-axis when you make a fist. If the axes are 
labeled consistent with this then we say you are using a right handed coordinate 
system. 

OK, let's now plot the point (2, 3, 2). First, locate a point 2 units from the 
origin on the x-axis. Now picture a plane which goes through this point, and is 
perpendicular to the x-axis. Repeat this for a point 3 units from the origin on the 
y-axis, and a point 2 units from the origin on the z-axis. Finally, the point (2, 3, 2) 
is at the intersection of the three planes you are picturing. 

Given the point (x, y, z) onecan "see" the quantities x, y, and z as in Figure 1-2. 
The quantity z, for example, is the distance from the point to the xy-plane. 
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Figure 1-2 Plotting the point (2, 3, 2) 


Problem 2 Which of the following coordinate systems are right handed? 





x 
NC 
ey 
e 
NC 
ay 
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Problem 3 Plot the following points on one set of axes: 


1. (,1, D 
2-4 id) 
3. et pos 


1.3 Introduction to Graphing 


We now turn back to the problem of visualizing a function of multiple variables. 
To graph the function f (x, y) we set it equal to z and plot all of the points where 
the equation z = f(x, y) is true. Let's start with an easy example. 
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EXAMPLE 1-2 

Suppose f (x, y) = 0. Thatis, f (x, y) is the function that always returns the number 

0, no matter what values of x and y are fed to it. The graph of z = f(x, y) = Ois 

then the set of all points (x, y, z) where z = 0. This is just the x y-plane. 
Similarly, now consider the function g(x, y) = 2. The graph is the set of all 

points where z = g(x, y) = 2. This is a plane parallel to the xy-plane at height 2. 


We first learn to graph functions of a single variable by plotting individual points, 
and then playing “connect-the-dots.” Unfortunately this method doesn’t work so 
well in three dimensions (especially when you are trying to depict three dimensions 
on a piece of paper). A better strategy is to slice up the graph by various planes. 
This gives you several curves that you can plot. The final graph is then obtained by 
assembling these curves. 

The easiest slices to see are given by each of the coordinate planes. We illustrate 
this in the next example. 


EXAMPLE 1-3 
Let's look at the function f(x, y) =x +2y. To graph it we must decide which 
points (x, y, z) make the equation z = x + 2y true. The xz-plane is the set of all 
points where y — 0. So to see the intersection of the graph of f (x, y) and the xz- 
plane we just set y = 0 in the equation z = x + 2y. This gives the equation z = x, 
which is a line of slope 1, passing through the origin. 
Similarly, to see the intersection with the yz-plane we just set x = 0. This gives 
us the equation z — 2y, which is a line of slope 2, passing through the origin. 
Finally, we get the intersection with the xy-plane. We must set z — 0, which 
gives us the equation 0 = x + 2y. This can be rewritten as y = — $x. We conclude 


this is a line with slope -4. 


The final challenge is to put all of this information together on one set of axes. 
See Figure 1-3. We see three lines, in each of the three coordinate planes. The graph 
of f (x, y) is then some shape that meets each coordinate plane in the required line. 
Your first guess for the shape is probably a plane. This turns out to be correct. We'll 
see more evidence for it in the next section. 


Problem 4 Sketch the intersections of the graphs of the following functions with 
each of the coordinate planes. 


1. 2x +2 Sy 
2. x? +y 
ey 
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SY 
<y 











Figure 1-3 The intersection of the graph of x + 2y with each coordinate plane is a line 
through the origin 


4. 2x? + y? 


5. / x2 + y? 


6 x? — y? 


1.4 Graphing Level Curves 


It’s fairly easy to plot the intersection of a graph with each coordinate plane, but 
this still doesn’t always give a very good idea of its shape. The next easiest thing 
to do is sketch some level curves. These are nothing more than the intersection of 
the graph with horizontal planes at various heights. We often sketch a “bird’s eye 
view" of these curves to get an initial feeling for the shape of a graph. 


EXAMPLE 1-4 

Suppose f(x, y) = x? + y?. To get the intersection of the graph with a plane at 
height 4, say, we just have to figure out which points in IR? satisfy z = x? + y? and 
z = 4. Combining these equations gives 4 = x? + y?, which we recognize as the 
equation of a circle of radius 2. We can now sketch a view of this intersection from 
above, and it will look like a circle in the xy-plane. See Figure 1-4. 

The reason why we often draw level curves in the x y-plane as if we were looking 
down from above is that it is easier when there are many of them. We sketch several 
such curves for z = x? + y? in Figure 1-5. 

You have no doubt seen level curves before, although they are rarely as simple 
as in Figure 1-5. For example, in Figure 1-6 we see a topographic map. The lines 
indicate constant elevation. In other words, these lines are the level curves for the 
function which gives elevation. In Figure 1-7 we have shown a weather map, with 
level curves indicating lines of constant temperature. You may see similar maps in 
a good weather report where level curves represent lines of constant pressure. 
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(b) 


Figure 1-4 (a) The intersection of z = x? + y? with a plane at height 4. (b) A top view 
of the intersection 


EXAMPLE 1-5 

We now let f(x, y) = xy. The intersection with the xz-plane is found by setting 
y = 0, giving us the function z = 0. This just means the graph will include the 
x-axis. Similarly, setting x = 0 gives us z = 0 as well, so the graph will include 
the y-axis. Things get more interesting when we plot the level curves. Let’s set 
z =n, where n is an integer. Solving for y then gives us y = 7. This is a hyperbola 
in the first and third quadrant for n > 0, and a hyperbola in the second and fourth 
quadrant when n < 0. We sketch this in Figure 1-8. 


Figure 1-5 Several level curves of z = x? + y? 
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Figure 1-6 A topographic map 
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Figure 1-7 A weather map shows level curves 
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Figure 1-8 Level curves of z = xy 
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Problem 5 Sketch several level curves for the following functions. 


1. 2x + 3y 
2. x*4-y 
3. /x? + y? 


4, x1— y? 





; 7 


Problem 6 The level curves for the following functions are all circles. Describe 


the difference between how the circles are arranged. 


Lx ty 
2. yx? + y? 
1 
4. sin(x? + y?) 


1.5 Putting It All Together 


We have now amassed enough tools to get a good feeling for what the graphs of 
various functions look like. Putting it all together can be quite a challenge. We 


illustrate this with an example. 


10 Advanced Calculus Demystified 











(b) 


Figure 1-9 Sketching the paraboloid z = x? + y? 


EXAMPLE 1-6 

Let f(x, y) = x? + y?. In Problem 4 you found that the intersections with the 
xz- and yz-coordinate planes were parabolas. In Example 1-4 we saw that the 
level curves were circles. We put all of this information together in Figure 1-9(a). 
Figure 1-9(b) depicts the entire surface which is the graph. This figure is called a 
paraboloid. 


Graph sketching is complicated enough that a second example may be in order. 


EXAMPLE 1-7 
In Figure 1-10 we put together the level curves of f (x, y) = xy, found in Example 
1-5, to form its graph. The three-dimensional shape formed is called a saddle. 


Problem 7 Use your answers to Problems 4 and 5 to sketch the graphs of the 
following functions: 


. 2x +3y 
x?- y 

ye ee i 
Kx y 


vy? 


AR 5 Mon 
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Figure 1-10 Several level curves of z = xy piece together to form a saddle 


1.6 Functions of Three Variables 


There is no reason to stop at functions with two inputs and one output. We can also 
consider functions with three inputs and one output. 


EXAMPLE 1-8 
Suppose 


fy, z =x +xy +yz 


Then f (1, 1, 1) = 3 and f(0, 1,2) = 4. 


To graph such a function we would need to set it equal to some fourth variable, 
say w, and draw a picture in a space where there are four perpendicular axes, x, 
y, z, and w. No one can visualize such a space, so we will just have to give up 
on graphing such functions. But all hope is not lost. We can still describe surfaces 
in three dimensions that are the level sets of such functions. This is not quite as 
good as having a graph, but it still helps give one a feel for the behavior of the 
function. 


EXAMPLE 1-9 
Suppose 


fe, yz) =x +y +z 


To plot level sets we set f (x, y, z) equal to various integers and sketch the surface 
described by the resulting equation. For example, when f (x, y, z) = 1 we have 


12x +y +z 
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This is precisely the equation of a sphere of radius 1. In general the level set 
corresponding to f (x, y, z) = n will be a sphere of radius y/n. 


Problem 8 Sketch the level set corresponding to f (x, y, z) = 1 for the following 
functions. 


1. fy j= +yz 
2. f(x,y, 2) =x? -y — 2? 


1.7 Parameterized Curves 


In the previous sections of this chapter we studied functions which had multiple 
inputs, but one output. Here we examine the opposite scenario: functions with one 
input and multiple outputs. The input variable is referred to as the parameter, and 
is best thought of as time. For this reason we often use the variable f, so that in 
general such a function might look like 


c(t) = Cf (t), 8) 


If we fix a value of t and plot the two outputs we get a point in the plane. As f 
varies this point moves, tracing out a curve, C. We would then say C is a curve that 
is parameterized by c(t). 


EXAMPLE 1-10 
Suppose c(t) = (cos f, sint). Then c(0) = (1, 0) and c (3) = (0, 1). If we continue 
to plot points we see that c(t) traces out a circle of radius 1. Indeed, since 


cos? t + sin? t = 1 


the coordinates of c(t) satisfy x? + y? = 1, the equation of a circle of radius 1. In 
Figure 1-11 we plot the circle traced out by c(t), along with additional information 
which tells us what value of t yields selected point of the curve. 


EXAMPLE 1-11 

The function c(t) = (cos t?, sin t?) also parameterizes a circle or radius 1, like the 
parameterization given in Example 1-10. The difference between the two param- 
eterizations can be seen by comparing the spacing of the marked points in Figure 
1-11 with those of Figure 1-12. If we think of f as time, then the parameterization 
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Figure 1-11 The function c(t) = (cost, sin t) parameterizes a circle of radius 1 


depicted in Figure 1-12 represents a point moving around the circle faster and 
faster. 


EXAMPLE 1-12 


Nowletc(t) = (t cos t, t sint). Plotting several points shows that c(t) parameterizes 
a curve that spirals out from the origin, as in Figure 1-13. 





Figure 1-12 The function c(t) = (cos t°, sint?) parameterizes a circle of radius 1 in a 
different way 
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Figure 1-13 The function c(t) = (t cost, t sint) parameterizes a spiral 


Parameterizations can also describe curves in three-dimensional space, as in the 
next example. 


EXAMPLE 1-13 

Let c(t) — (cost, sint, t). If the third coordinate were not there then this would 
describe a point moving around a circle. Now as ft increases the height off of 
the xy-plane, i.e., the z-coordinate, also increases. The result is a spiral, as in 
Figure 1-14. 





Figure 1-14 The function c(t) = (cos f, sin t, t) parameterizes a curve that spirals 
around the z-axis 
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Problem 9 Sketch the curves parameterized by the following: 
1. (t,t) 

2 Ug) 

3. (0,1) 

4. (t,t) 

5. (cos2t, sin 3t) 

Problem 10 The functions given in Examples 1-10 and 1-11 parameterize the same 


circle in different ways. Describe the difference between the two parameterizations 
for negative values of t. 


Problem 11 Find a parameterization for the graph of the function y = f (x). 


Problem 12 Describe the difference between the following parameterized curves: 


1. c(t) = (cost, sint, t?) 
2. c(t) = (cost, sint, 1) 
3. c(t) = (t cost, t sint, t) 


Quiz 
Problem 13 


1. Determine if the coordinate system pictured is left or right handed. 
Az 





2. Let fx y) = e 
a. Sketch the intersections of the graph of f (x, y) with the xy-plane, the 
xz-plane, and the yz-plane. 
b. Sketch the level curves for f (x, y). 
c. Sketch the graph of f (x, y). 


3. Sketch the curve parameterized by c(t) = (2cost, 3 sint). 


This page intentionally left blank 


CHAPTER 2 





Fundamentals of 
Advanced Calculus 


2.1 Limits of Functions of Multiple Variables 


The study of calculus begins in earnest with the concept of a limit. Without this 
one cannot define derivatives or integrals. Here we undertake the study of limits of 
functions of multiple variables. 

Recall that we say lim f(x) = L if you can make f(x) stay as close to L as 
you like by restricting x to be close enough to a. Just how close “close enough” is 
depends on how close you want f(x) to be to L. 

Intuitively, if lim f (x) = L we think of the values of f (x) as getting closer and 


closer to L as the value of x gets closer and closer to a. A key point is that it should 
not matter how the values of x are approaching a. For example, the function 


X 


f= > 


[x] 
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does not have a limit as x — 0. This is because as x approaches 0 from the right 
the values of the function f (x) approach 1, while the values of f(x) approach —1 
as x approaches 0 from the left. 

The definition of limit for functions of multiple variables is very similar. We say 


lim X,y)—L 
T » 


if you can make f (x, y) stay as close to L as you like by restricting (x, y) to be 
close enough to (a, b). Again, just how close "close enough" is depends on how 
close you want f (x, y) to be to L. 

Once again, the most useful way to think about this definition is to think of the 
values of f(x, y) as getting closer and closer to L as the point (x, y) gets closer 
and closer to the point (a, b). The difficulty is that there are now an infinite number 
of directions by which one can approach (a, b). 


EXAMPLE 2-1 
Suppose f (x, y) is given by 


X 


fA, y= 





x+y 


We consider lim x,y). 
odis d y) 


First, let’s see what happens as (x, y) approaches (0, 0) along the x-axis. For all 
such points we know y = 0, and so 


2 Z j 


fœ = TG x 





Now consider what happens as (x, y) approaches (0, 0) along the y-axis. For all 
such points we have x = 0, and so 


x 0 
x+y y 





fœ, y= 


We conclude the values of f(x, y) approach different numbers if we let (x, y) 


approach (0, 0) from different directions. Thus we say m o f (x, y) does not 
x,y) (0,0) 


exist. 


Showing that a limit does not exist can be very difficult. Just because you can 
find multiple ways to come at (a, b) so that the values of f(x, y) approach the 
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same number L does not necessarily mean li md (x, y) = L. There might be 


some way to approach (a, b) that you TM t tried that gives a different number. 
This is the key to the definition of limit. We say the function has a limit only when 
the values of f(x, y) approach the same number no matter how (x, y) approaches 
(a, b). We illustrate this in the next two examples. 


EXAMPLE 2-2 
Suppose f (x, y) is given by 


zo 
Tic SB 


As we let (x, y) approach (0, 0) along the x-axis (where y — 0) we have 


xy 0 
=. .=6 
fœ, y) = su 


Similarly, as we let (x, y) approach (0, 0) along the y-axis (where x — 0) we have 


xy 0 
uà 


= —0 
x+y y 


fA, y= 


But if we let (x, y) approach (0, 0) along the line y = x we have 


F = xy A od 
T x-y 2x? 2 





So once again we find lim f(x, y) does not exist. 
(x, y)— (0,0) 


Our third example is the trickiest. 


EXAMPLE 2-3 
Let 


x’y 


fœ, y) = x yl 


As (x, y) approaches (0, 0) along the x-axis (where y — 0) we have 
x^y 


fœ, y)= say r 


1 
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As (x, y) approaches (0, 0) along the y-axis (where x = 0) we have 


xy 0 


iO) "aa ok = ge 


If we approach (0, 0) along the line y = x we get 
x^y X x 
As x approaches 0 we have 


lim = 
x0 x2 +1 





So far it is looking like perhaps 


x^y 


lim  ———— 
(,)2.(0,0) x4 + y? 
since, as (x, y) approaches (0, 0) along the x-axis, the y-axis, and the line y — x, 


the values of f(x, y) approach 0. But what happens if we let (x, y) approach (0, 0) 


along the curve y — x?? In this case 
x?y x 1 


vty? xx 2 





fox.» = 


We can evaluate the limit of this function as x approaches 0 by dividing the numer- 
ator and denominator by x‘. 


xt 
lim = lim 


= 1 
x>0x24+x4  x201 +1 
x^ 








So again the limit does not exist. 


Problem 14 Show that the following limits do not exist: 


1. 


2 
lin == 
(x,y) (0,0) 30" 


2. 


-3 3 
(x,y) (0,0) * +9 
x+y 


lim 
(x,y) > (0,0) «/x?-+y? 


y? 


lim : 
(x,y) (0,0) VH 


3. 





4. 
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2.2 Continuity 


We say a function f (x, y) is continuous at (a, b) if its limit as (x, y) approaches 
(a, b) equals its value there. In symbols we write 


lim | f(x, y) = f(a, b) 


(x, y) (a, b) 


Most functions you can easily write down are continuous at every point of their 
domain. Hence, what you want to avoid are points outside of the domain, where 
you may have 


1. Division by zero. 
2. Square roots of negatives. 
3. Logs of nonpositive numbers. 


4. Tangents of odd multiples of 7. 


In each of these situations the function does not even exist, in which case it is 
certainly not continuous. But even if the function exists it may not have a limit. 
And even if the function exists, and the limits exist, they may not be equal. 


EXAMPLE 2-4 
Suppose 


Ec 
Vy 


There is no zero in the denominator when (x, y) = (1, 1), so f(x, y) is contin- 
uous at (1, 1). 


fœ, y) = 


EXAMPLE 2-5 
Evaluate 


l xi? 
lim. —— — — 
G.3) 0,0) x2 + y2 4-1 
There are no values of x and y that will make the denominator 0, so the function 


is continuous everywhere. Since the value of a continuous function equals its limit, 
we can evaluate the above simply by plugging in (0, 0). 


xy? 0 


G,»29(0,0 x?-- y? --1 0+1 
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Problem 15 Find the domain of the following functions: 


1. Œ 


. x—y 


2. /y-x? 
3. In(y — x) /x — y 


Problem 16 Consider the function 





x+y? 
f(x, y) = | met C 9) FO, 0) 
i 1 


(x, y) = (0, 0) 


Is f (x, y) continuous at (0, 0)? 


Quiz 
Problem 17 


1. Show that the function 


x sin y 
x? + y? 





fœ, y= 


does not have a limit as (x, y) — (0, 0). 
2. Is the function 
sey (x, y) z (0,0) 


son= |i (x, y) = (0,0) 





continuous at (0, 0)? 
3. Find the domain of the function 


1 
fay) =In—. 
A =y 


CHAPTER 3 





Derivatives 


3.1 Partial Derivatives 


What shall we mean by the derivative of f (x, y) at a point (xo, yo)? Just as in one 
variable calculus, the answer is the slope of a tangent line. The problem with this 
is that there are multiple tangent lines one can draw to the graph of z = f(x, y) at 
any given point. Which one shall we pick to represent the derivative? The answer 
is another question: “Which derivative?” We will see that at any given point there 
are lots of possible derivatives; one for each tangent line. 

Another way to think about this is as follows. Suppose we are at the point (xo, yo) 
and we start moving. While we do this we keep track of the quantity f(x, y). The 
rate of change that we observe is the derivative, but the answer may depend on 
which direction we are traveling. 

Suppose, for example, that we are observing the function f(x, y) = x?y, while 
moving through the point (1, 1) with unit speed. Suppose further that we are travel- 
ing parallel to the x-axis, so that our y-coordinate is always one. We would like to 
know the observed rate of change of f(x, y). Since the y-coordinate is always one 
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the values of f (x, y) that we observe are always determined by our x-coordinate: 
f (x, 1) = x?. The rate of change of this function is given by its derivative: 2x. 
Finally, when x — 1 this is the number 2. 

The above is a particularly easy computation. Given any function, if you are 
traveling in a direction which is parallel to the x-axis then your y-coordinate is 
fixed. Plugging this number in for y then gives a function of just x, which we can 
differentiate. Here's another example. 


EXAMPLE 3-1 

We compute the rate of change of f(x, y) = x?y? at the point (1, 2), when we 
are traveling parallel to the x-axis. During our travels the value of y stays fixed 
at 2. Hence, the values of the function we are observing are determined by our 
x-coordinate: f(x,2) = 8x?. The derivative of this function is then 24x”, which 
takes on the value 24 when x is one. 


What if we wanted to repeat our computations, with different values of y? It 
would be helpful to keep the letter “y” in our computations, and plug in the value at 
the very end. Notice that when we plugged in a number for y it became a constant, 
and was treated as such when we differentiated with respect to x. If we leave the 
letter y in our computations we can still treat it as a constant. 


EXAMPLE 3-2 

Let f(x, y) 2 x +xy + y?. We wish to treat y as a constant, just as if we had 
plugged in a number for it, and take the derivative with respect to x. Recall that the 
derivative of a sum of functions is the sum of the derivatives. So we will discuss 
the derivatives of each of the terms of x + xy + y? individually. 

There is no occurrence of y in the first term, so it is particularly easy. Its derivative 
is just one. 

The second term is a little trickier. Since we are treating y as a constant this is of 
the form const-x. The derivative of such a function is just const. So the derivative 
of xy is just y. 

Finally, the quantity y? is also a constant. The derivative of a constant is zero. 
Hence our answer is 1 + y. 


Notice in the above example that if we thought of x as constant, and y as the 
variable, then the derivative would have been very different. We need some notation 
to tell us what is changing and what is being kept constant. We use the symbols ar 
to represent the partial derivative with respect to x. This means x is considered a 
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variable, every other letter is a constant, and we differentiate. Similarly, the notation 
2f. is the partial derivative with respect to y. 


ày 
EXAMPLE 3-3 
Again we let f (x, y) = x + xy + y?. We compute 


9 
af ky 
dy 


There are various other ways to think of partial derivatives that are useful. One 
is completely algebraic. Recall that the derivative with respect to x of a function 
f (x) of one variable is defined as the limit 


Hd ou fG + Ax) — f(x) 
— im 


dx v Ax—0 Ax 





The partial with respect to x is defined similarly. Just remember that y is kept 
constant, so that f(x, y) really becomes a function of just x. We then apply the 
above definition to get: 





Of — i fœ + Ax, y)-— fœ, y) 
— = hm 
OX Ax>0 Ax 


The partial with respect to y is defined similarly 





Of qus GE mA») e FO) 
I— im 


dy  Ay>0 Ay 


There is yet another way to think about the partial derivative. We began this 
section by claiming that the derivative will still represent the slope of a tangent 
line. In the following figure we see the graph of an equation z = f(x, y). Through 
the point (xo, yo) in the xy-plane there is also drawn a vertical plane parallel to 
the xz-plane. The intersection of this vertical plane with the graph is a curve. 
The slope of the tangent line to this curve is exactly the value of the partial 
derivative with respect to x at (xo, yo). To see the partial derivative with respect 
to y we would have a similar picture, where the vertical plane is parallel to the 
yz-plane. 
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(Zo; Yo) 


Problem 18 Compute a (2, 3) and AQ. 3) for the following functions. 


I. x 4- xy 
2. xIny 


3. X /XY 


Problem 19 Compute af (x, y) and 2 (x, y) for the following functions. 


I. x^y* 

2. 5 
Problem 20 For the function f (x, y) = —x + xy? — y? find all places where both 
af 


9f 
2s and ay are zero. 


3.2 Composition and the Chain Rule 
3.2.1 COMPOSITION WITH PARAMETERIZED CURVES 


Suppose we have a parameterized curve ġ (t) = (x(t), y(t)) in the plane. That is, 
for a given value of t we are given the numbers x(t) and y(t), which we visualize 
as a point in the plane. We can also take these two numbers and plug them in to 
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a function f(x, y). The result is the composition f(@(t)). Notice that only one 
number goes in to this function, and only one number comes out. 


EXAMPLE 3-4 
Let f(x, y) 2 x? + y?. Let (t) = (t cost, t sint). Then 


f($(t)) = f(tcost,tsint) = t? cos? t +t’ sit = 1? 


There are various ways to visualize the composition. One is to imagine the graph 
of z = f(x, y) in three dimensions. Then draw the parameterized curve $ (f) in 
the xy-plane, and imagine a vertical piece of paper curled up so that it sits on this 
curve. Now mark where the paper intersects the graph of f (x, y), and unroll it. The 
result is the graph of f (o (t)). 


Unroll 


ge SN 
































Since f ($(t)) is a function of one variable we can talk about its derivative just 
as if we were in a first term calculus class. But what we want to do here is relate 
it to the derivatives of f (x, y), x(t), and y(t). The formula we end up with is the 
multivariable calculus version of the “chain rule” 


_ of dx | Of dy 


EXAMPLE 3-5 
Let f (x, y) and $ (t) be defined as in the previous example. We wish to compute 
£f (t)) when t = £. Note first that 


JT 3r TT 
«G-( 12 4) 
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To use the chain rule we must compute the partial derivatives of f(x, y) at this 
point. First note that 


0 
oe y) = 2x, and 
Ox 


à 
A. (y) e 2y 
y 








Thus, 
af [ 3m m NEL 
: = , and 
əx 12 12 6 
af [ 3n m oT 
dy\ 12°12) 6 


We also need the derivatives of x(t) and y(t) when t = A 
x'(t) 2 cost — t sint, and 
y'(t) = sint + t cost 


Thus, 





Finally, we have 


(Ie. 22 








6// ax dt ðydt 
NE EL: V3 T "nm MEE; 
| 6 2 12 6 12 12 
Ta 


Now that we have gone through the pain of using the chain rule to compute the 
derivative, it should be pointed out that the same answer could have been found 
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much faster directly. We know from the previous example that f ($(t)) = t?. So the 


derivative function is 21. Plugging in t = 4 then immediately gives (2) (z) = 7 
Problem 21 Let ¢ (t) = Qt, t?). Suppose you don't know what f (x, y) is, but you 


know a (2, 1) = 6 and ar Qo, 1) = —1. Compute £ fa). 


Problem 22 Suppose you don't know what $(t) = (x(t), y(t)) is, but you know 
$D = (1,3), xQ)2-2, and y'(2)=1. Let f(x,y) 2 x?y. Compute 
arf (2). 


3.2.2 COMPOSITION OF FUNCTIONS OF MULTIPLE 
VARIABLES 


We now look at the idea of composition with more complicated types of functions, 
as in our next example. 


EXAMPLE 3-6 
Let x(u, v) = uv and y(u, v) = u? + v?. Suppose f (x, y) = x + y. Then we may 
form the composition f (x(u, v), y(u, v)) as follows. 


f(x, v), y(u, v)) = x(u, v) + y(u, v) = uv +u? + v? 


Notice that the result is a function whose input is a pair of numbers, u and v, and 
whose output is a single number. Hence, we may talk about the partial derivatives, 
with respect to u and v, of the function given by composition. The result is given 
by another variant of the chain rule. As in the previous example, we will begin with 
the functions x(u, v), y(u, v), and f (x, y). The following formulas give the partial 
derivatives of the composition f (x (u, v), y(u, v)). 


af  Ofüx | df ay 
ðu dx du dy Ou 
af af ax df ay 
2s Bx de” By ov 





At this point it is mathematically meaningless to think of terms like “dx” as 
entities in themselves that can be canceled, but thinking this way may help you 
remember the above formulas. Note that the formula for the partial derivative of f 
with respect to u, for example, has two terms. If you cancel the “quantity” dx in 
the first term, and dy in the second, you are left with ar in both. 
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EXAMPLE 3-7 

Let f(x, y) = x? + xy. Suppose we don't know x(u, v) or y(u, v), but we know 
x(1,2) = 3, y(1,3) = —2, x (1, 2) = —1, and 2 (1, 2) = 5. Then we may use the 
chain rule to compute Fc, 2). To do this we will need to know af and x at the 
point (x(1, 2), yd, 2)) = (3, —2). 


9f 3, —2) = (2)(3) + (-2) = 4, 9f a, —2)=3 
Ox dy 


We now employ the chain rule: 
of _ofdx  Əf dy 
ðu Ox du Əy ðu 
= (4)(-1) + 6)5) 
= 11 





Problem 23 The following table lists values for a function f (x, y) and its partial 
derivatives. 


aD -3 | 2 
(1,2) 5 1| 1 


(2, 1) 2 | o| 7 
9,5) .. —1 2| 3 
Q,3)| M 1|—1 
(3,2) 2 1| 0 











Let x(u, v) = uv and y(u, v) =u + v’. Find the following partial derivatives at 
the indicated points. 


1. 1. at (u, v) = (1,2) 
2. 1. at (u,v) = (2, 1) 


Problem 24 Let f(x,y) = sin(x + y), x(u,v) =u +v, and y(u, v) —-u — v. 
Find the following quantities when (u, v) — (5s T). 


L. IUS y) 
af 
2x 
af 
As 
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3.3 Second Partials 


In the previous sections we saw that the partial derivatives of a function f (x, y) are 
also functions of x and y. We can therefore take the derivative again with respect 
to either variable. 


EXAMPLE 3-8 
Let f(x, y) 2 x?y?. The partial derivatives are oF = 2xy? and af = 3x?y?. We 


can take the partial derivative again of both of these functions with respect to either 


variable: 
9 (0 ð (0 
ae. of = 2yP, L— af = 6xy? 
Ox \ ox dy \Ox 
9 (0 9 (a 
rn =6xy, — 9f = 6x?y 
ax \ dy dy \ dy 

As is customary, we adopt the following shorthand notations for the second 

derivatives: 
3 (If) _ PF 9 (af\ PF 
dx \dx)] üx2' dy NOx)  dyax 
ð (9f 9f d (IF) _ PF 
ax Voy) àx8üy dy \day)  8y2 


ae 








02f 








The quantities =~ and ;-7- are called the mixed partials. The above example 
illustrates an amazing fact: Under reasonable conditions the mixed partials are 
always equal! The "reasonable conditions" are only that the mixed partials exist 
and are themselves continuous functions. The proof of this goes back to the limit 


definition of the partial derivative. 


Problem 25 Find all second partial derivatives of the following functions. 


1. f(x, y) =xy 
2. fo, y) =x - y? 
3. f(x, y) = sin(xy?) 


m$ x 
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Quiz 


Problem 26 Let f (x, y) 2 x?y + x3y?. 


Í; 


2. 
3. 


Find af and 2. 
If b(t) = (t,t — 1), then what is f (b(t))? 
Suppose you don’t know what w(t) = (x(t), y(t)) is, but you know w(2) = 


(1, 1), 4x (2) — 3, and 2» (2) = l. Find the derivative of f(w(t)) when 
t —2. 


. Suppose x and y are functions ofu and v, x(u, v) = u? + v, and y(1, 1) = 1. 


What would a have to be when (u, v) = (1, 1), if £ = 12? 


CHAPTER 4 





Integration 


4.1 Integrals over Rectangular Domains 


The integral of a function of one variable gives the area under the graph and above 
an interval on the x-axis called the domain of integration. For functions of two 
variables the graph is a surface. We will interpret the act of integration as that 
of finding the volume below the surface, and above some region in the xy-plane. 
Eventually, we will examine how to do this with very general-shaped regions. We 
begin by looking at rectangles. 


i (a, b) 





- 


x 





Let R be the region of the x y-plane pictured above. Suppose we want to find the 
volume of the region in IR? which sits above R, and below the graph of z = f (x, y), 
as in the following figure. We do this by following these steps. 
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1. Begin by choosing a grid of points {(x;, y;)} in R, so that the horizontal and 
vertical spacings between adjacent points are Ax and Ay, respectively. 


2. Connect these grid points to break up R into rectangles. Note that the area 
of each rectangle is Ax Ay. 


sy 








| — (zi; V;) 








Ay{ Zi 





























SY 


roe 


Ax 





3. We now draw a box of height f (x;, y;) above each such rectangle to get a 
figure which approximates the desired volume. The volume of each box is 
its length x width x height, which is f (x;, yj) Ax Ay. 
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4. We now add up the volumes of all of these boxes to obtain the quantity 
3 f Gi, ypAxAy 
j i 


5. Finally, we repeat this process indefinitely, each time using a grid with Ax and 
Ay smaller and smaller. The corresponding figures that we get approximate 
the desired volume more and more. In the limit we get what we’re after, 
which we denote as f f f (x, y) dx dy: 

R 


J| 16 ctn tim, Jim Y fes spans 
8 | 


Basic properties of summation and limits allow us to rearrange the above equation 
as follows: 


J| 16 cr im Y [im Y reos Ay 
R J ! 


The quantity in the brackets Jim 5; f Gi, yj) Ax is exactly the definition of 
what you get from f(x, y) when you fix y and integrate x, as in the following 


example. 


EXAMPLE 4-1 
Let f(x, y) = x? + y?. If we fix y = 3 then this is the function f(x) = x? + 27. 
We may now integrate this function over some range of values of x, such as [0, 2]: 


2 
. 1 
Jim dof Gi yp Ax = a +27 dx = ae + 27x 
' 0 


2 





8 
=~ +56 
x=0 3 


If we don’t substitute the value 3 for y, but we still think of y as a constant, very 
little changes: 


2 
8 

=_42y3 
gre 





2 
1 
Jetra Ix + Px 
3 x=0 
0 


In the above example, notice that when we think of y as a constant and we 
integrate with respect to x, our answer is a function of y. We may now integrate 
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this new function with respect to y. This is precisely what we are instructed to do 
in the limit we obtained for ff f (x, y) dx dy. 
R 


| f(x, y) dx dy = lim am, Doe nas] Ay 


i Ay—0 j 

= ij 

ua [fends Ay 
J Lo 

b a 

= [ tena dy 
0 0 

EXAMPLE 4-2 


Let f(x, y) = xy’. Suppose Q is the rectangle with vertices at (1, 1), (2, 1), (1, 4), 
and (2, 4). To find the volume under the graph of f and above this rectangle we 
wish to compute 


4 2 


MELLES fa dy 


Q 1 1 


We work inside the brackets first. To do this integral we must pretend y is a 
constant. 





We now have 


4F2 T i 
2 2 3 
dx | dy = | -yd = 
Jif» y J» y= 3 
1 Li 1 


Note that in the definition of f f f (x, y) dx dy we could have rearranged the 
R 





terms so that the integral with respect to y was inside the brackets. The fact that the 
answer does not depend on what order you integrate is called Fubini's Theorem. 
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EXAMPLE 4-3 

Let S be the rectangle with vertices at (1, 0), (2, 0), (1, 1), and (2, 1). We integrate 
the function f(x, y) = ex. We compute the integral of f (x, y) over S by doing the 
integral with respect to y first: 


2r1 
f [él fea 
5 1 Lo 


EP 





Il 
To 
— 
= 
Nn 

Mie 
| 
= 
mea 
> 





—x —x 1 2 
= —xe e 2% l 
o3 2 3 
|. e e 2 


To save space the integral PE Jg f(x, y)dxjdy is often written as 
b pd 
fo f. fœ, y) dx dy. 


Problem 27 Compute the following: 


13 

1. f f x t xy? dx dy 
0 2 
11 

2. f [ x2? dy dx 
AO 
22 

3. f f cos(x + y) dx dy 
0 0 


Problem 28 Let R be the rectangle in the xy-plane with corners at (—1, —1), 
(—1,0), Q, —1), and (2,0). Find the volume below the graph of z — xy and 
above R. 
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Problem 29 Let R be the rectangle with vertices at (0, 0), (1, 0), (0, 1), and (1, 1). 
Find a formula for f f x" y" dx dy. 
R 


Problem 30 Let V be the volume below the graph of e^. and above the rectangle 
with corners at (0, 0), (2, 0), (0, 3), and (2, 3). Find the area of the intersection of 
V with the vertical plane through the point (1, 1) which is parallel to the xz-plane. 


4.2 Integrals over Nonrectangular Domains 


In the previous section we saw how to compute the volume under the graph of a 
function and above a rectangle. But what if we are interested in the volume which 
lies above a nonrectangular area? For example, suppose now that R is the region 
of the x y-plane that is bounded by the graph of y = g(x), the x-axis, and the line 
x = 1. We would like to find the volume of the figure V which lies below the graph 
of z = f(x, y) and above the region R. 














Inthe previous section we computed volume by breaking up the region in question 
into small boxes. Our strategy here is to cut it into “slabs” parallel to the yz-plane. 
To specify the location of each slab we must give a value of x. Hence, the volume 
of each slab, calculated as an integral with respect to y, is a function of x. Adding 
the slabs up is just like integrating with respect to x. 

We follow these steps: 


1. Choose points {x;} in the interval [0, 1]. 


2. Compute the area A (xj) of the intersection of V with the plane x = x;. We 
do this by plugging x; in for the variable x in the function f(x, y) and 
then integrating with respect to y. Notice that y ranges from 0 to g(x;) on 
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this domain. So the requisite area is given by 


gi) 
nu J fend 


glxi) 











3. The volume of a thin slab is then given by 


gGi) 
AG ihe J Ce 
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4. Adding up the volumes of all of the slabs thus gives 


8 (xi) 


Y | fany) dy Ax 
E 
5. Finally, the desired volume is obtained from this quantity by choosing smaller 
and smaller slabs. This is equivalent to taking a limit as Ax tends toward O. 


gi) 
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eaves f(x, y) dy dx 
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The key to understanding this formula is the limits of integration. Each slab is 
parallel to the y-axis. The area of the side of the slab is thus computed as an integral 
with respect to y. But the range of values that y can take on depends on what the 
value of x is. Hence, the limits of integration of the inner integral depend on x. 


EXAMPLE 4-4 
Let R be the region in the xy-plane bounded by the graph of y — x?, the x-axis, 
and the line x — 1. We compute the volume below the graph of z — xy? and above 


R as follows: 
ig 
Volume = | fe dy dx 
0 0 





EXAMPLE 4-5 

We now use the above ideas to compute a more complicated volume. Let Q be the 
region of the xy-plane bounded by the graphs of y = x? and y = 1 — x?. We wish 
to determine the volume that lies below the graph of z = x? + y and above Q. 
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The above figure depicts the region Q. Notice that for a fixed value xo of x the 
values of y range from x? to 1 — x?. This tells us the limits of integration for the 
inner integral. Now notice that the smallest and largest possible values of x are 
where the graphs of x? and 1 — x? coincide. To find this we solve 


x21-x? 


Which implies 
2x? =1 
And hence, 


2 
MM 
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This tells us the limits of integration for the outer integral. We now compute the 
desired volume 
y2 


1—x? 


Volume = f [ 9*0 


42 x? 





42 Advanced Calculus Demystified 


vs 


1 1 
xê + ian 3x + 





X 
2 


w| — 
Js 


lS 


Sometimes the problem is set up so that integrating with respect to x first is more 
natural, as the following example illustrates. 


EXAMPLE 4-6 

Let S be the region of the xy-plane bounded by the graph of y = x?, the y-axis, 
and the lines y = 1 and y = 2. We will integrate the function f(x, y) = x over this 
domain. 

The domain S is pictured in the following figure. Notice that for a fixed value 
yo of y the range of values that x can take on goes from 0 to ,/yo. This tells us the 
limits of integration for the inner integral. The smallest value yo can be is 1 and the 
largest value is 2. These are the limits for the outer integral. 


Ay 








We now compute 
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Often you can set up an integral in two ways: one so that dx comes first, and one 
so that dy comes first. Sometimes you'll find that only one of these ways makes 
the problem accessible. We illustrate this in the next example. 


EXAMPLE 4-7 
Let P be the region of the xy-plane bounded by the graph of y — x, the x-axis, and 


the line x = 1. We wish to integrate the function f (x, y) = e-* over P. 





The region P is depicted above twice. Notice from the figure on the left that if 
we fix a value yo of y then x can range from yọ to 1. The smallest possible value 
for yo is O and the biggest is 1. This tells us the limits of integration when we set 


up the integral with dx first. 


1 1 


[feo dx dy 


0 y 


If, on the other hand, we fix a value xo of x then y can range from 0 to x, as in 
the figure on the right. The smallest possibility for x is 0 and the largest possibility 
is 1. This tells us that the integral can also be set up as 


] x 
| fe dy dx 
0 0 
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It is only possible to do the second of these integrals, which we compute as 
follows: 


X 


1 
fe dydx = J ye * 
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(by u-substitution) 





Problem 31 Evaluate: 


2 


1. 2xy? dx dy 


oO 


EU 


ete ay Ot 


ety dy dx 


Problem 32 Let T be the region of the xy-plane bounded by the graph of y = 


x? — x — 2 and the x-axis. Integrate the function f (x, y) = x? over T. 


Problem 33 Let R be the region of the positive quadrant of the x y-plane bounded 
by y = x and y = x?. Set up two different integrals for f (x, y) over R. 


Problem 34 Evaluate the following integral by switching the order of integration: 


zx? m 
J [ o» dy dx 
0 x 


4.3 Computing Volume with Triple Integrals 


The cylindrical figure pictured below has a base with area A and has height 1. The 
volume is therefore the product of these quantities, namely A - 1 — A. But the top 
surface of this figure is the graph of the equation z — 1. The integral of this function 
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gives the volume below the graph, which we just found out was A. In other words, 
if you want to find the area of a region R of the xy-plane, you just need to compute 
the integral of the function f (x, y) = 1 over R: 


Area(R) — Hi 1 dx dy 
R 





EXAMPLE 4-8 

Suppose R is the region below the graph of y = g(x), above the x-axis, and between 
the lines x — a and x — b. Then we can find the area of R by evaluating the double 
integral: 


b g(x) 
Area(R) m J 1 dy dx 
a 0 
But this just reduces to 
b b 
foe dx = [eo dx 


which should, of course, be familiar from first term calculus as the area under the 
graph of y = g(x). 


Just as one can compute area by evaluating the double integral of the function 
f(x, y) = 1, one can also compute volume by a triple integral of the function 
f(x,y,z) = 1. The tricky part usually involves finding the limits of integration. 
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EXAMPLE 4-9 
Let V be the region of IR? bounded by the planes y = x, z = x + y, z = 0, y = 0, 
and x = 1, as pictured below. 





M 


We set up a triple integral to compute the volume of V. The order of integration 
will be dz dy dx. To find the limits for the innermost integral (the one with respect 
to z) we fix x and y and observe that z can range from 0 to x + y. For the second 
integration we just fix x, and observe that y can vary from 0 to x. Finally, the 
smallest value that x can take on is 0 and the largest is 1, determining the limits of 
the outermost integral. We now compute 
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Volume(V) — 1 dz dy dx 
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Problem 35 Use a triple integral to find the volume of the solid bounded by the 
graph ofz = 1 — y^, and the planes z = 0, x = 0, and x = 1. 


Problem 36 Find the volume of the solid bounded by the graph of z = 1 — x — y?, 
and the coordinate planes z = 0 and x = 0. 


Problem 37 Set up a triple integral to find the volume which lies below the 
paraboloid z = 1 — x? — y? and above the xy-plane. 


Problem 38 Write a triple integral which will compute the volume bounded by the 
sphere x? + y? +z? = 1. 


QUIZ 
Problem 39 


1. Evaluate the following integrals: 


23 
a. f f cosQx + y) dx dy 
12 


11 
b. f f A y! dy dx 
0 x 


2. Set up an integral for the volume which lies between the cone 4/1 — x? — y? 
and the xy-plane. 
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CHAPTER 5 





Cylindrical and 
Spherical Coordinates 


5.1 Cylindrical Coordinates 


A coordinate system is a systematic way of locating a point in (some) space by 
specifying a few numbers. The coordinate system you are most familiar with is 
called rectangular coordinates. In IR? this works by specifying x, y, and z, which 
represent distances along the x-, y-, and z-axes, respectively. 

Rectangular coordinates are extremely cumbersome when trying to define certain 
common shapes, such as cylinders or spheres. A second way to determine the 
location of a point is to adapt polar coordinates to three dimensions. In this case 
we give a location in the xy-plane by a value of r and 0, and then the height off of 
the xy-plane by a value of z. 
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ge 





Converting from cylindrical coordinates to rectangular coordinates is as easy as 


it is in polar coordinates. The z coordinate is exactly the same in both systems. To 
get x and y you just use the formulas 


x —rcos0 
y=rsing 


EXAMPLE 5-1 


We locate the point r = 2, 0 = c, and z = 7 in rectangular coordinates: 


x =2cos = V3 
m 
= 2sin— = 1 
y 6 


So the desired point is at (v3, 1,7). 


Problem 40 Write the rectangular coordinates of the following points. 


1. (86,2) = Q, £, -1) 
2. (r,0,z) = (0, 2,3) 


3. (r,0,z) = (4, $,0) 
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5.2 Graphing Cylindrical Equations 


To get areal feel for cylindrical coordinates it is helpful to look at the graphs of a 
few simple equations. The graph of something like z = 3 in cylindrical coordinates 
is the same as the graph of z = 3 in rectangular coordinates (it’s the same z!): a 
horizontal plane that hits the z-axis at 3. The next two examples are a bit more 
interesting. 


EXAMPLE 5-2 

Consider the cylindrical equation r = 3. This is the set of all points in IR? which are 
exactly 3 units away from the z-axis. The shape is a cylinder of radius 3 centered 
on the z-axis. 


EXAMPLE 5-3 

We now examine the cylindrical equation 0 = 7. This describes the set of all points 
that make an angle of 7 with the xz-plane. The shape is a plane through the z-axis 
that is halfway between the x z-plane and the yz-plane. 


z 


T 


More complicated equations involve relationships between r, 0, and z, as in the 
next example. 


EXAMPLE 5-4 

Consider the equation z = r°. To visualize the graph we are being instructed to 
find all of the points in R? whose z-coordinate is the square of their distance to 
the z-axis. The most relevant feature of the equation, however, is the lack of the 
coordinate 0. This tells us that no matter what 0 is we will see the same picture. 
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Following this observation we consider the half plane which meets the z-axis at 
some arbitrary angle. Take this plane out of the picture and look at it head-on, as 
in the figure below. As you move left and right in this plane it is the quantity r that 
is changing. As you move up and down it is z that changes. Hence, it makes sense 
to label our axes r and z in this plane. 


z 











In the half plane at the right we have depicted the graph of z = r?, a famil- 
iar parabola. Now, the key to understanding the complete picture is to place this 
parabola in every half plane at the left. The result is the surface that you get when 
you rotate the graph of y = x? (in R2) around the y-axis. It is called a paraboloid, 
and is depicted below (see also Figure 1-9 for a computer-generated picture). 





Problem 41 Describe the shape of the following cylindrical equations. 
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1. 0 
2.0 
3. z 
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4.z=r 
5. 24+r=1 


Problem 42 Sketch the graphs of the following equations. 


l. z=Ł 


2. r? -z =1 
3. z2 —r?=1 


5.3 Spherical Coordinates 


Cylindrical coordinates are very convenient for describing cylinders and other sur- 
faces of revolution. In this section we explore spherical coordinates, which are 
useful for describing many other shapes. 

To locate a point using spherical coordinates you need to know three quantities: 
p (pronounced “Row”), 0, and $ (pronounced “Fee”). The quantity 0 is precisely 
the same as in cylindrical coordinates. o is the distance from the origin. $ is the 
angle with the z-axis, as pictured below. 


z 
© 
/ 
b/p 
/ 


Given the spherical coordinates of a point we would like to be able to convert to 
rectangular coordinates. The easiest quantity to find is z. Consider the triangle in 
the figure below. The side adjacent to the angle $ is precisely z. The hypotenuse is 
p. So cos ó = z: Solving for z then gives 


Zz = pcos ġ 
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z = pcos @ 


yo r 





To find x and y itis helpful to first find the quantity r from cylindrical coordinates. 
Going back to the triangle in the above figure, we see that the side opposite the 


angle $ is precisely r. Hence, sing = z and thus 
r=psing 
We know from the previous section that x = r cos @ and y = r sin 0. Substituting 
for r gives 


x = psin$ cos0 


y = psing sing 


EXAMPLE 5-5 
We find the rectangular coordinates of the point where p = 2, 0 = 7, and o = 


ala 


. m 1 
x-—2sin—cos— = = 
6 3 2 
_u . T J3 
y = 2sin — sin — = — 
6 3 2 


z= 2c0s = = V3 
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Problem 43 Convert the following to rectangular coordinates. 


1. (0,0,0) = (0,2, Z 
2. (0, 6,6) = G, x, 2) 
3. (0,66) = (4, 2, € 


5.4 Graphing Spherical Equations 


We now explore the graphs of several equations in spherical coordinates. 


EXAMPLE 5-6 

The simplest equation to understand is something like = 2. This describes the 
set of all points that are precisely 2 units away from the origin. In other words, this 
is a sphere of radius 2. 


If 0 is a constant then the equation is the same in spherical coordinates and 
cylindrical coordinates; a plane through the z-axis. The previous example shows 
that when p is constant we get a sphere. In the next example we see what happens 
when ¢ is held constant. 


EXAMPLE 5-7 
Consider the equation @ = 7. This describes the points which make an angle of 7 
with the z-axis. This is a circular cone, centered on the z-axis. 


Az 
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As with cylindrical coordinates, things get more difficult when equations involve 
relationships between multiple coordinates. 


EXAMPLE 5-8 
We graph the equation o = $. Notice that the coordinate 0 does not appear in 
the equation. This tells us that the intersection of the graph with every half plane 
incident to the z-axis is the same. (Recall that we encountered a similar type of 
graph when we looked at some cylindrical equations.) This, in turn, tells us that the 
graph is a surface of revolution. 

In the following figure we have sketched one half plane. In this half plane we 
have indicated how a spiral is the graph of p = $. 


W 


To form the complete graph we must now take this spiral and rotate it about the 
z-axis. 


In the previous example we saw that when 0 did not appear in a spherical equation 
the graph was a surface of revolution. In the next example we explore what happens 
when ¢ does not appear in a spherical equation. 


EXAMPLE 5-9 

We now explore the graph of o = 0. We do this again by investigating various half 
planes incident to the z-axis. This time the pictures are not all the same. When 
0 = 0, for example, we are looking at the graph of o = 0, which is just a point. 
When 0 = 5 we are looking at p = 5, which is half of a circle of radius 7. As 0 
increases we see larger and larger half circles incident to the z-axis, as in the figure 
below. 
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Putting it all together is challenging. We have done this for you in the computer- 
generated picture below. The result is reminiscent of a sea shell. 





Problem 44 Describe the graphs of the following spherical equations. 


I. p=0 
2. $255 
3. ó$ —m 
4. psing =2 
5. pcos@ = 2 


Problem 45 Sketch the graphs of the following spherical equations. 


l. p=sing 
2. p = cos 


1 
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Quiz 
Problem 46 


1. Give rectangular coordinates for the following points: 


a. The point with cylindrical coordinates r = 1, 0 = §, a 2. 


ndz = 
b. The point with spherical coordinates p = 2, 0 = Z, and $ = 


AIA 


2. Sketch the graphs of the following: 


a. The cylindrical equation r = cos 0. 


b. The spherical equation 0 = $. 


CHAPTER 6 





Parameterizations 


6.1 Parameterized Surfaces 


In Chapter 1 we encountered parameterized curves. These were functions that 
had one input (the parameter) and multiple outputs. One visualizes such curves 
by drawing the range of the function. In this section we introduce parameterized 
surfaces. These are functions which have two inputs and multiple (usually three) 
outputs. Again, we visualize such functions by picturing their range. This is most 
often some surface in IR?. 


EXAMPLE 6-1 
We explore the parameterization given by 


y(u, v) = (u, v, u? +0) 
Note that for each point in the range of w the z-coordinate is obtained from the 
x and y coordinates by squaring them and adding. In other words, if a point is in 


the range of y then it is also on the graph of z = x? + y?. This is the paraboloid 
depicted in Figure 1-9. 
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EXAMPLE 6-2 
We now look at a slight modification of the previous example. 


y(u, v) = (2u, v, u^ de v?) 


Note that this parameterization is obtained from the previous one by doubling 
the x-coordinate of every point. This has the effect of stretching the graph in the 
x-direction. Intersections with horizontal planes will thus be ellipses, as opposed 
to circles as in the previous example. 


Example 6-1 illustrated another important idea. Suppose we want to represent the 
graph of a function z — f (x, y) as a parameterized surface. Then we may simply 
write 


y(u, v) = (u, v, f (u, v)) 


Note that the z-coordinate of every point in the range is obtained from the x- and 
y-coordinates by plugging them into f. Hence, each point in the range is indeed 
on the graph of f. 

It is not much more difficult to represent graphs in other coordinate systems as 
parameterized surfaces. The trick is to always translate to rectangular coordinates. 


EXAMPLE 6-3 

Suppose we want to write the graph of the cylindrical equation z — r? as a parame- 
terized surface. First, we write down the translation from cylindrical to rectangular 
coordinates: 


x —rcos0 

y=rsing 

Z-—Z 
Now we substitute r? for z: 

x —rcos0 

y=rsing 

peg 


Hence, a parameterization may be written as follows: 


wir, 0) = (rcosO, rsin6, r?) 
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EXAMPLE 6-4 
Suppose we want to find a parameterization for the graph of the spherical equation 


p = coso 
First we write down the translation from spherical to rectangular coordinates: 


x = p sin $ cos0 
y = psin $ sin 


z = pcosó 
Now we substitute cos $ for p: 


x = cos ¢singcos 0 
y = cos $ sin à sind 


Z = cos$ cos 
This gives us the parameterization 
V (0, $) = (cos sin $ cos 8, cos o sing sin 0, cos $ cos $) 


Problem 47 Find parameterizations for the graphs of the following equations. 


lz-x? 
2r=07 
3. p=" 


Problem 48 The following are parameterizations of the graphs of rectangular, 
cylindrical, or spherical equations. Find these equations. 


1. $(p,0) = (psinO cos0, p sin? 0, p cos Q) 
2. $(x, z) = (x, x + xz, Zz) 
3. $(r, 0) = (r cos0, rsinO, sinr) 


T9 c Advanced Calculus Demystified 


6.2 The Importance of the Domain 


When you are trying to specify a particular shape, the domain of the parameteriza- 
tion can be just as important as the function. We illustrate this with an example. 


EXAMPLE 6-5 

Suppose we want to parameterize the part of a sphere of radius 1 that lies in the 
first octant (i.e., the region of R? where x, y, z > 0). The equation of a sphere is 
particularly simple in spherical coordinates, so we start here. Such an equation is 
p = 1. Writing down the translation to rectangular coordinates and substituting 1 
for p then yields the parameterization 


V (0, $) = (sin $ cos 9, sing sin 0, cos $) 


But how do we make sure we just get the part of the sphere we want? The answer 
is to restrict the values of 0 and ¢ that can be plugged into ¢: 


0<0<5, 0<¢< 


R2] 3 


These ranges constitute the domain of $. 


EXAMPLE 6-6 

We parameterize the portion of the graph of z — r? that lies inside the cylinder 
x? + y? = 1. The equation whose graph we are interested in is given in cylindrical 
coordinates, so it is natural to start there. The restriction on the domain can also be 
naturally expressed in cylindrical coordinates, so we do not anticipate any problems. 
In Example 6-3, we found a parameterization for the desired surface: 


wir, 0) = (r cosO, rsin6, r?) 


The points of R? that lie inside the cylinder x? + y? = 1 all have r < 1. Hence, we 
restrict the domain of our parameterization to 
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Problem 49 The function 
V (0, 6) = (sing cos 9, sing sin 0, cos $) 


parameterizes a sphere of radius 1. How would you restrict the values of 0 and $ 
to just get the portion of the sphere where 


Xy Z 
2.x <0,y20,z <0 


Problem 50 Describe the difference between the shapes with the following param- 
eterizations. 


1 Wa, y=, y, yx? +y); 0<x<1,0<y<1 
2. Y(r,0) = (rcosé,rsiné,r);0<r<1,0<0<27 
3. Wr, 0) = (r cos0,r sin, r); 0 <r <1,0<0<7x 











6.3 This Stuff Can Be Hard! 


Parameterizations are extremely powerful tools that can describe very complicated 
surfaces. But with such great power can come great difficulty. There may be no 
obvious way to look at a complicated parameterization and know what shape it 
describes. For example, consider the parameterization 


$ (t, u) = ((4 + cost + 2 sin u) cos(2u), 
(4 + cos t + 2 sin u) sin(2u), sint + 2 cos u) 


O0<t<2nm, O<u<27 


No expert would be able to look at an equation so complicated and immediately 
know what it looks like. This is where computers become an invaluable tool. Plug- 
ging this equation into a standard graphing program reveals the following beautiful 
picture. 
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Despite the apparent difficulty in understanding such complicated functions, 
there are a few tricks that can help. These are basic ways in which you can take 
a parameterization and modify it to suit your needs. For example, if you multiply 
some coordinate by a constant, this has the effect of stretching the parameterized 
shape by that factor. If you add a constant to some coordinate the shape is moved 
in that direction. 


EXAMPLE 6-7 
Consider the following parameterization: 


W(6, $) = (2sin$ cos0 + 1,3 sing sin0 — 2, cos d) 
0x0-x2n, O0<¢d<az 
This looks much like a parameterization of a sphere from spherical coordinates. 
But the x-coordinate was multiplied by 2, and then 1 was added to it. Similarly, 


the y-coordinate was multiplied by 3 and then 2 was subtracted from it. Thus the 
parameterized shape is an ellipsoid with center at (1, —2, 0). 
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Problem 51 Describe how the following parameterizations are related to graphs 
of functions in rectangular, cylindrical, or spherical coordinates. 

L Uy, z) = (^ 4z2y,z- 1) 

2. V(x, y) = Qu x, x t siny) 

3. V(0,z) = (2cos0,3sin0 + 1, z — 1) 


Problem 52 Parameterize the surface whose level curves are the ellipses shown 
below. 





6.4 Parameterized Areas and Volumes 


Earlier in the text we looked at parameterized curves in IR? and R°. This chapter 
began with parameterized surfaces in IR?. Now we look at parameterized areas in 
IR? and parameterized volumes in R°. Often these types of parameterizations begin 
with writing down the translation from some coordinate system to rectangular 
coordinates, and then restricting the domain. 
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EXAMPLE 6-8 
The translation from polar coordinates to rectangular coordinates is given by 


x —rcos0 
y=rsing 
We can use this to parameterize the area that is both inside the unit circle and in the 
first quadrant: 
Wir, 0) = (rcos@,rsin@) 


srzi, 0<0<5 


Note that it is the restriction on the domain that guarantees that the correct area is 
parameterized. 


EXAMPLE 6-9 
Let R be the region below the graph of y = f(x), above the x-axis, and between 
the lines x = a and x = b. Then R is parameterized by 


W(x, t) = (x, tf x) 
a<x<b, red 


In three dimensions one usually begins with rectangular, cylindrical, or spherical 
coordinates. 


EXAMPLE 6-10 
Consider the parameterization 


Wir, 0, z) = (rcosO,rsinO, z) 
l<r<2, 0x0-z2m, 0<z<1 
The function itself is just the translation from cylindrical to rectangular coordinates. 


The restriction on the domain means that we have only parameterized the region 
between cylinders of radii 1 and 2, and between the planes z = 0 and z = 1. 
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Of course, as in the previous section, one can distort and move shapes by multi- 
plying and adding constants to each coordinate. 


EXAMPLE 6-11 
We examine the parameterization 


V (0,0, 6) = (2p sin ġ cos0, psing sind, p cos $) 


T 


952251 Uses, (EOS 


wl 3 


Without the 2 in the first component of the parameterization this would be one- 
eighth of a solid ball. The 2 just stretches it in the x-direction, making it one-eighth 
of an ellipsoid. 


Problem 53 Sketch the areas of the plane parameterized by the following. 


1. V(r,0) = (rcosé,rsin@),0<r<1,0<86,< = 
2. W(r, 8) = (rcosOo,2rsin0),1 xr x2,0x0,x x 








Problem 54 Describe the volumes given by the following parameterizations. 


1. V(p,0,0) = (psiné$cosO, psing sind, pcosQ),.1 < p x 2,0 0 x 2st, 
06-5; 
2. V(r.0O,z) = (rcosé,rsiné,2z+1),0<r<1,0<0<27,0<z<l 








Problem 55 Find a parameterization for the volume which lies below the cone 
z =r, above the xy-plane, and inside the cylinder x? + y? — 1. 
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Quiz 
Problem 56 


1. Find a parameterization for the portion of the cylinder of radius 2, centered 
on the z-axis, which lies below the graph of z = 2r and above the xy-plane. 


2. Sketch the region in R? parameterized by the following: 
V (0,0, 0) = (2p sin ġ cos0, p sing sin0, pcos $) 


where0<p<1,0<9 <5, and0<ġ<m. 





CHAPTER 7 





Vectors and Gradients 


7.1 Introduction to Vectors 


Mathematically, there is very little difference between a vector and a point. In two 
dimensions, for example, a point is really just a pair of numbers. We visualize this 
pair as a dot in a plane. A two-dimensional vector is similar. Technically, it is also 
just a pair of numbers. But this time we visualize it as an arrow (Figure 7-1). This 
arrow does not have a preferred location. The two numbers that specify it only 
give its magnitude (1.e., size) and direction. This can be convenient, since certain 
algebraic operations with vectors can be seen visually by moving them around. 

Some differences between points and vectors are purely cosmetic. For example, 
the numbers that you give to specify a point are called its coordinates. The numbers 
that specify a vector are its components. The main difference between points and 
vectors is that we can do algebra with vectors. This is a topic studied in detail in a 
class on linear algebra. We will need some of the ideas from linear algebra here. 

The first algebraic operation we will study with vectors is addition. This is very 
easy: to add two vectors we simply add its components. 


Copyright © 2007 by The McGraw-Hill Companies. Click here for terms of use. 
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Figure 7-1 The point (1, 2) and the vector (1, 2) 


EXAMPLE 7-1 
We compute the sum of the vectors (2, 3) and (—1, 2) as follows: 


(2,3) 4(=1,2) S(2 = 1,3 4 2) = (1,5) 


We can visualize the operation of addition as follows. To add a vector V to a 
vector W we slide W so that its “tail” (the base of the arrow) coincides with the 
"tip" of V. The result of V + W is then the vector whose tail is at the tail of V and 
whose tip is at the tip of W. 





Negatives are equally easy. If V is a vector then we let —V denote the vector 
whose components are the negatives of the components of V . The picture, of course, 
is a vector which has the same length as V , but points in the exact opposite direction. 
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EXAMPLE 7-2 
If V = (2,3) then —V = (—2, —3). 


Finally, we define subtraction by combining addition and negatives: 
V-W=V+(-W) 
The result of subtraction is therefore exactly what you would expect: the vector 
whose components are the difference between those of V and W. 
The magnitude (or length) of a vector is easily computed by the Pythagorean 


Theorem. We denote the magnitude of a vector V by the absolute value bars | - |. 
Hence, 


(a, b)| = va? +b? 


EXAMPLE 7-3 


42 42 JB. dum. dic 
(GF 


A vector whose length is one is often called a unit vector. 








Another common operation is multiplying a vector by a number. This is referred 
to as scalar multiplication. Like addition and subtraction, this is done component- 
wise. That is, the product of a scalar (i.e., a number) c and a vector V is a vector 
whose components are c times the components of V. That is, 


c (a, b) = (ca, cb) 


Geometrically, scalar multiplication produces a vector which points in the same 
direction, but whose magnitude has been altered by the given scale factor. For 
example, if we multiply a vector by the scalar 2, then we get a new vector, pointing 
in the same direction, which is twice as long. 


EXAMPLE 7-4 

Suppose we would like to find a unit vector that points in the same direction as 
(3, 4). The trick is to multiply it by some scalar, and arrange it so that the resulting 
vector has length one. In other words, we would like to find a number c such that 


Ic (3,4)|=1 
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A quick computation gives us 


le(3, 4)| = |(3e, 4c)| 
= VG tee 
= V25c2 
= 5c 


Hence, 5c = 1, and so c = i. We now multiply the vector (3, 4) by this number 
to get the answer (3, 1). 


Problem 57 For the following vectors V and W compute V + W, —V, and V — W. 


1. V (1,6, W = (6, 1) 
2. V = (0,0), W = (1, 2) 
3. V=(-1,1), W = (1, —2) 


Problem 58 Find the magnitude of the following vectors. 


W 
< 


1. (2, 
2. (1, 3) 
Problem 59 Find a unit vector that points in the same direction as (2, 1). 


Problem 60 Find a unit vector that is perpendicular to the vector (5, 12). 


7.2 Dot Products 


After we have mastered the basic operations of addition and subtraction we move on 
to multiplication. Unfortunately, we do not have anything quite as straightforward. 
It is tempting to define the product of V and W as the vector whose components are 
the products of the components of V and W. This particular operation, however, 
does not turn out to be terribly useful. 

Instead we define the product of two vectors V and W to be the sum of the 
products of the components. This means the result of multiplication is a single 
number, as opposed to a vector. This product is called the dot product, since it is 
always denoted with a dot. In symbols, we write 


(a, b) - (c, d) = ac + bd 
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EXAMPLE 7-5 
(2,3) - (3,4) =64+12= 18 


The usefulness of the dot product comes from the fact that it has a nice geometric 
interpretation. Suppose the angle between vectors V and W is 0. Then, 


V-W=|V||W|cosé 


This follows from some tricky trigonometry as follows. Consider the following 
triangle: 


Now recall the Law of Cosines, which says that if a triangle has side lengths a, b, 
and c, and the angle between a and b is 0, then 


c? = a? + b? — 2ab cos 0 
If we apply this to the above triangle we get 
IV — WP = |V]? + Wf? — 2]|VIIW| cosó 
which we can rearrange to 


| ve qwp -iv - wp 


|V||W| cos 8 
2 





Now suppose V = (v1, v2) and W = (wj, w2). Then, 


VP +IWP-IV—WP optus +w? + w3 — (vi — w)? — (o — w)? 
2 2 
_ 2v, Ww, + 2v2.wW2 
= 


= vw + v2W2 
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And so 
V . W = vw; + vw: = |V||W|cosé 


In retrospect, this formula is very reasonable. Consider two extreme cases: 
Case 1. Suppose we take the dot product of V with itself. By definition, 
V. V = vv 4+ 1 — v? ++ v2 


But in this case 0 = 0, so cos 0 = 1. Hence, 





[V||V| cos8 = |V||V| = ET v2 vi + v) =v ty 


Case 2. The other extreme case is when V and W are perpendicular. The slope 
of the line containing V is m Similarly, the slope of the line containing W is 
a If these lines are perpendicular then their slopes are negative reciprocals, so 
an = —i". Cross multiplying and bringing everything to one side of the equation 
then gives U1U1 + v2W2 = O. But the quantity on the left is precisely the definition 
of V - W. To verify the cosine formula just note that if V and W are perpendicular 


then cos @ = 0, so |V||W|cos@ = 0 as well. 


EXAMPLE 7-6 
The cosine of the angle between vectors (2, 3) and (—1, 2) can be found using the 
dot product 


oe uod 3 eae ee ed 
= |(2, 3}||(—1, 2)| cos 8 
= /744+9/1+4+4cos 6 


65 cos 0 
Hence, 
Bene 4 _ 44/ 65 
| 4/65 65 
EXAMPLE 7-7 


Since the cosine of the angle between two vectors is zero if and only if the vectors 
are perpendicular, the dot product can be used to detect this. For example, 


Q,3)-(-6,4) = —12+ 12 = 0 


so these vectors must be perpendicular. 
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Problem 61 Compute the dot product of the following pairs of vectors. 
1. V = (2,4, W = (—3, 1) 
2. V = (0,7), W = (5, 2) 
3. V = (—2, —1}, W = (6, —3) 
Problem 62 Compute the cosine of the angle between the following pairs of vectors. 


1. V = (2,4), W = (-3, 1) 
2. V = (-2,-1), W = (6, —3) 


Problem 63 Suppose V and W are vectors of length 3. What are the smallest and 
largest possible values for V - W ? What is the smallest possible value of |V - W|? 


Problem 64 Use the dot product to decide if each of the following pairs of vectors 
are perpendicular. 


1. V = (2,4), W = (—4, 8) 
2. V = (2,4), W = (8, —4) 
3. V = (-2,-1), W = (-3, 6) 
4 V =(-2,-1), W = (6, —3) 


7.3 Gradient Vectors and Directional 
Derivatives 


In this section we return to the concept of differentiation of functions of two vari- 
ables. Recall that the partial derivative af was the rate of change of the function 
f (x, y) if you walk in the positive x-direction with unit speed. Similarly, the partial 
derivative -< ar is the rate of change if you walk in the positive y-direction with unit 
speed. We are left with the obvious question, “What is the rate of change of f(x, y) 
if you walk in some other direction?” 

Let V = (a, b) represent a vector which points in the direction of travel. The 
length of V will be one, to reflect the fact that we are walking with unit speed. We 
already know two things: 


0 

Rate of change of f (x, y) in direction (1,0) = af = 197. +0— of 
Ox ET dy 
of of of 





Rate of change of f(x, y) indirection (0,1) = =0—+1 
dy Ox dy 
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It is reasonable, then, to expect 


Rate of change of f(x, y) in direction (a, b) = a— + b— 
x 


EXAMPLE 7-8 
Suppose f(x, y) = xy’. Then the rate of change of f(x, y), at the point (2, 1), in 


the direction (2, By is 


Dn p XH o, js M E 





The rate of change of f(x, y), in the direction of V, is called a directional 
derivative, and is denoted as Vy f (x, y). Hence, we have the formula 


af pot 
Vy f(x, y) — an T e. 


But notice that the right side of this equation can be rewritten as a dot product: 


0 
viens | f T) 


ax’ dy 


We can further shorten this equation by coming up with new notation for (E ; of ). 


Henceforth we will call this vector the gradient of f, and denote it as V f. The 
equation for Vy f then becomes 


Vyf =V- Vf 


EXAMPLE 7-9 
The gradient of the function f (x, y) = xy? at the point (x, y) is 


Vf, y) = Qxy*, 3x^y*) 
The gradient at the point (1, 1), then, would be V f (1, 1) = (2, 3). 
Recall that the rate of change of a function is also the slope of a tangent line to its 


graph, as long as you are traveling with speed one. Here's a nice way to think about 
the situation. Suppose you are climbing a mountain, and you have a good trail map 
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in your hands. Let (x, y) be your coordinates when you locate yourself on the map. 
The function f(x, y) is your elevation at that point. Now turn your body to face the 
direction V (on your map). If you sight up or down so that your gaze just grazes 
the mountainside then you are looking along the tangent line whose slope is given 


7 Acme Trail Map 


7 eT vy 





You are here 


Recall that the dot product is given by the product of the magnitudes of the two 
vectors, times the cosine of the angle between them. If we fix the point we are at 
then V f is a fixed vector. If |V| = 1, then the only way to change Vy f = V -Vf 
is to change the angle between V and V f. The largest this quantity can be is when 
the value of cos 0 is largest. This happens when 0 — 0. We conclude that the largest 
value of Vy f is given by 


Vy f =V- Vf —|VIIVf[cos8 = |V f| 


and that this value is attained when V (the direction we are facing) coincides with 
the direction of V f. In other words, if you are on the mountainside and you want 
to face directly uphill you should point yourself in the direction of the gradient 
vector. When you do this and sight along the mountainside the slope you see is the 
magnitude of the gradient vector. 


EXAMPLE 7-10. 
Let f (x, y) =e” . At the point (x, y) the gradient vector is 


Vf(x,y)- (yer, 2xye™) 
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So, at the point (2, 3) we have V f (2, 3) = (9e!8, 12e!8). The largest slope of any 
tangent line to the graph at the point (2, 3) is then given by 





|(9e!8, 12e!8)| = /(e!8)2 + (12e!8)2 = 15e!8 


What if you were standing on the mountainside and wanted to face the direction 
you would have to travel to keep your elevation constant? In other words, how would 
you find the direction of your level curve? If you were facing such a direction you 
would be looking along a horizontal line, i.e., a line whose slope is zero. The only 
way for Vy f = V - Vf to be zero is for V and Vf to be perpendicular. 


EXAMPLE 7-11 

Suppose again f(x, y) = e*”. In the previous example we saw V f (2,3) = 
(9e!8, 12e!5). A vector which points in a direction perpendicular to this would be 
(4, —3) (Check this!). Hence, this vector is tangent to a level curve at the point (2, 3). 


Problem 65 For each of the functions f (x, y) below answer the following ques- 
tions: 


e Find the gradient vector V f (x, y). 
* Calculate the rate of change of f (x, y) at the point (1, 1) in the direction 
3 4 
(5, 5} 
* Find a unit vector that points in the direction of the maximum rate of change 
at the point (1, 1). 


e Find the largest slope of any tangent line at (1, 1). 


* Find a unit vector that lies in a line tangent to a level curve through (1, 1). 


I. xIny 
2. 2x + 3y 
Vy 4 xy? 


Ww 


7.4 Maxima, Minima, and Saddles 


Ata local maximum or a local minimum of a graph the tangent plane is horizontal. 
Another way to say this is that the slope of every tangent line to a local maximum 
or minimum is zero. But recall that the slope of a tangent line at (x, y) in direction 
V is given by V - V f. The only way for this to be zero for every possible V is if 
V f is the zero vector. 
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EXAMPLE 7-12 
We show that f(x, y) = x?y? cannot have a maximum or minimum at any point 
where x and y are nonzero. The gradient is 


Vf(x,y)-— Qxy?, dy] 


If x and y are nonzero then this is not the zero vector. This tells us that there is 
a direction where the slopes of tangent lines are nonzero, and hence we cannot be 
at a local maximum or minimum. 


Unfortunately, just because the gradient is the zero vector it does not necessarily 
mean that there is a local maximum or minimum. 


EXAMPLE 7-13 
The gradient of f(x, y) = x? + y? and g(x, y) = x? — y? is (0, 0) at the origin. 
The function f (x, y) has a minimum there, while g(x, y) has a saddle. 


In a first-semester calculus class you learned to detect local maxima and minima 
by a second-derivative test. We would like to do the same thing here. The problem, 
of course, is that there are four second partial derivatives! To keep track of all this 
information we often write them in a matrix, as follows: 








af 3? f 
ax? s | 
of Pf 
QyOx 8yl | 


Now we examine this matrix for several functions whose graph is familiar. Each 
of these has a gradient vector equal to zero at the origin. 


1. x? + y?. This function has a local minimum at the origin. The matrix of 
0 
0 2| 


2. —x? — y?. This function has a local maximum at the origin. The matrix of 


second partials is 


d partials i —2 0 
second partials is | ġ 2| 


3. x? — y?. This function has a saddle at the origin. The matrix of second 


tials i 2 0 
parüalsis|o — ^ 
4. xy. This function also has a saddle at the origin. The matrix of second 
partials is : 


1 0 


80 Advanced Calculus Demystified 


The first two give us a clue as to the quantity we would like to look at. Consider 
the product of the upper-left and lower-right entries of the matrix. For the maximum 
and minimum above this quantity is positive and for the first of the above saddles 
it is negative. However, this alone would not be enough to distinguish maxima and 
minima from saddles, as the second of the saddles shows. To compensate we must 
subtract the product of the upper-right and lower-left entries, yielding the formula 


of af of of 
əx? 0y? axdy 8yOx 





However, since the mixed partials are equal we can shorten this to 





af af (2 i 
ax? dy? Ox dy 


This is indeed the right quantity to look at, in the sense that if it is greater than 
Zero you have a maximum or minimum, and if it is less than zero you have a saddle. 
Unfortunately, if it is zero you have no information; you may be at a maximum, 
minimum, saddle, or something much more bizarre. Nonetheless, we will single 
this out as our first test. 








Test 1 





t| 


oy ðxəðy 


_ fef (Spy ; B 
Let D(x, y) = 55 — . Suppose that at some point (xo, yo) we have V f (xo, yo) = 


(0, 0). 


D(xo, yo) > 0 => f(x, y) has a local max or min at (xo, yo). 
D(xo, yo) < 0 => f(x, y) has a saddle at (xo, yo). 


D(xo, yo) = 0 > No information about f (xo, yo). 





EXAMPLE 7-14 

We find the location of all saddles of the function f (x, y) = x? + 2xy + 3y?. First, 
we will need to narrow down the possibilities by finding the critical points. To do 
this we find the gradient. 


Vf = (2x + 2y, 2x + 9y?) 
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Setting this equal to the vector (0, 0) gives us the system of equations 


2x+2y=0 
2x + 9y? 20 


The first equation tells us that x = —y. Plugging this into the second equation 
then gives 


—2y 495? 20 


Either y — 0 (and hence x — 0) or we can divide both sides of this equation by y 
to get 


—24+9y=0 


Solving then gives us y = i. Hence, we have critical points at (0, 0) and (-$. 2). 
To determine which of these are saddles we compute the matrix of second partials: 





ef af 

ðx? axdy | [2 2 
ef Pf =|; | 
QyOx 3y? | 


And so 
D(x, y) = 36y — 4 
We now check each critical point: 


D(0, 0) = —4 < 0 > (0, 0) is a saddle. 


D s 4-02 RAT local | 
ao e4- ——,— | isa local max or min. 
9'9 9'9 a loca 


When D(xo, yo) > O it would be nice to have a second test to determine whether 
(Xo. yo) is a local maximum or a local amu, Such a test can be easily guessed 
from our prototypical oe f(x,y) =x? +y? and f(x, y) = —x? — y?. 


Notice that in both cases fs and 5 rf have the same sign. When this sign is positive 


we have a local minimum and won it is negative we have a local maximum. This 
is precisely our second test. 
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Test 2 








2 \? 
ILeti/D(Gx. = a d (24). Suppose that at some point (xo, yo) we have V f (xo, yo) = 


(0, 0) and D(xo, yo) > 0. 
j : 
3290 yo) > O => f(x, y) has a local min at (xo, yo). 


2 
d yo) < 0 2 f(x, y) has a local max at (xo, yo). 








EXAMPLE 7-15 

Recall the function f(x, y) = x? + a + 3y? from the previous example. We 
found critical points at (0, 0) and (-3 ne $). and determined that at (-$. 2) there 
e : local maximum or a local minimum. To determine which we need only look 


at 74 Ê. Since this was 2, and 2 > 0, we conclude that at this critical point there is a 
local minimum. 


It is important to keep in mind that if D (xo, yo) = 0 then we have no information 
about the nature of f (xo, yo). We illustrate this in the next example. 


EXAMPLE 7-16 
Consider the following functions: 


1. fi, y) =x +y 
2. fy) — y 
3. fo, y) 2 xy? 
In each case the only critical point is at (0, 0) and D(0, 0) = 0. But at (0, 0) in the 


first case there 1s a local minimum, in the second there is a local maximum, and in 
the third there is a saddle. 


Problem 66 Find the local maxima, minima, and saddles ofthe following functions. 


I. xy -2x - 3y — 6 
2 xX — xy +y? 


Problem 67 For the function sin(x + y) show that D(x, y) = 0 at every point 
(x, y). Does this function have maxima, minima, or saddles? 
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Problem 68 /f, for some point (xo, yo) you know D(x, yo) > 0 and IT (xo, yo) > 


0 then show that Lt (xo, yo) > 0. (This tells us that you may use ef (xo, yo) instead 


a2 
of S (xo, yo) when distinguishing maxima from minima.) 


7.5 Application: Optimization Problems 


An important application of the mathematics of the previous section is to optimiza- 
tion problems. You encountered several such problems in a first semester calculus 
class. You can now handle much more complicated situations. We illustrate this 
with an example. 


EXAMPLE 7-17 

We find the largest volume of a box with an open top, and surface area 100 m?. 
First, we let the dimensions of the box be a, b, and c, with c the height. Then the 
surface area is given by 


S.A. = 100 = ab + 2ac + 2bc 


Solving this equation for c gives 


100 — ab 
e —L——— 
2(a +b) 


The volume, of course, is the quantity abc. Substituting for c then gives us 


100 — ab 


Volume = ab ————— 
2(a 4- b) 


Thinking of volume, then, as a function V(a, b) we now search for critical 
points. The gradient is a bit tedious to compute, but with some work you can show 
it simplifies to 





—2a?b? + 200b? — 4ab? —2a2b? + 200a? — 4a? b 
vva» =( ET á ee -— 


4(a +b)? ' 4(a +b)? 


Each component of this vector can only be zero if the numerators are zero, so 
we get two equations: 


0 = —2a?P? + 200b? — 4ab? 
0 = —2a?P? + 200a? — 4a? b 
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We may safely assume that neither a nor b are zero. Dividing the first equation 
by 2b? and the second by 2a? then gives 
0 = —a? + 100 — 2ab (7-1) 
0 = —P? + 100 — 2ab 


Solving both equations for 2ab and setting them equal to each other gives us 
—a?^ + 100 = —b° + 100 


and thus (since we are only interested in positive values of a and b) we may conclude 
a = b. Combining this with Equation 7-1 tells us 


0 = —a? + 100 — 2a? 
Solving this for a reveals 


_ 103 
3 


The fact that a — b then gives us 


104/3 


b= 
3 


Finally, our expression for c in terms of a and b above yields 


5/3 
Q——— 
3 


The volume is thus 


5004/3 


Volume — abc — 9 


7.6 LaGrange Multipliers 


Suppose you wanted to look for the maximum value of f(x) = pe — x, on the 


interval [—2, 3]. This function has two critical points, and the one at x = —lisa 
local maximum. But to conclude that f(—1) = - is the maximum value of f (x) on 
the indicated domain would be incorrect. This is because the maximum is attained 


at one of the endpoints of the domain, namely at x = 3. At this point f (3) = 6. 
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Figure 7-2 The function f(x, y) has a local maximum at the point q, but the 
absolute maximum on the domain D is the point p on the boundary 


The same kind of thing can happen with functions f (x, y) of two variables. Using 
techniques of the previous section we can look for, and analyze, critical points in 
the interior of the domain D of f. But if we are interested in the maximum value 
that f attains on D we may have to look at the points around the boundary of D, 
as in Figure 7-2. In the one-variable case this was easy, since it was just a matter of 
checking two points. When there are multiple variables there are an infinite number 
of points around the boundary of D, so we need some technique other than just 
plugging values into f. 

In Figure 7-3 we see the contours for the function f (x, y) of Figure 7-2, super- 
imposed on the domain D. As is customary, we will denote the points around the 
boundary of D as ðD. (Don't get this confused with a partial derivative.) Notice 
that the maximum and minimum values of f on dD occur at the places where the 
level curves of f are tangent to 0 D. 

The key fact that we will use is that the vectors V f are always perpendicular to 
the level curves of f. So, ata point p = (xo, yo) where a level curve of f is tangent 
to ðD we will see V f (xo, yo) perpendicular to dD. 

The next step is to find a new function, g(x, y), for which ð D is a level curve. The 
vector Vg is always perpendicular to all of its level curves. In particular, V g (xo, yo) 
will be perpendicular to ðD at (xo, yo). See Figure 7-4. 
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Figure 7-3 The contours of f(x, y) and the domain D. The absolute maximum of 
f on D is at the point p where a level curve is tangent to 0D 


Since V f (xo, yo) and V g(xo, yo) are perpendicular to the same curve they must 
be parallel to each other. If vectors V and W are parallel, then one is a scalar multiple 
of the other. In other words, there is a constant A such that V = AW. Putting all this 
together, we find that if a maximum (or minimum) of f occurs at a point (xo, yo) 
of ð D then there must be a constant A such that V f (xo, yo) = AV g(xo, yo). Using 
this information we get a finite number of candidate points at which to look for 
maxima and minima. 


EXAMPLE 7-18 
Let f (x, y) = 2x? + 3y?, and suppose D is the set of points (x, y) in R? satisfying 
x? + y? < 1. We will look for the largest value attained by f (x, y) on 8D. 





Figure 7-4 The vector V f is a scalar multiple of the vector Vg 
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To use the “method of Lagrange multipliers” the first step is to find a function 
g(x, y) for which 0D isa level curve. This is very easy: g(x, y) = x? + y?. Then the 
level curve g(x, y) = 1 is the curve x* + y* = 1, which is precisely the boundary 
of D. 

The next step is to look for points (xo, yo), and values of A, such that V f (xo, yo) = 
AV g(xo, yo). We thus compute 


Vf = (4x, 9y) 
Vg = (2x, 2y) 


So if (4x, 9y?) = ìà (2x, 2y) we have 4x = A2x and 9y? = A2y. The first of 
these equations tells us that either x = 0 or A = 2. 


Case I. x = 0. The boundary of D is the set of points satisfying x? + y? = 1. 
So if x — 0 then y = +1. This means that we must check the points (0, 1) and 
(0, —1) for maxima and minima. 


Case 2. à = 2. Combining this with the equations 9y? = A2y then tells us 9y? = 
4y. So again there are two cases: either y = 0 or 9y = 4, in which case y = 3. 
Subcase A. y — 0. Similar to Case 1 above, the boundary of D is the points 
satisfying x? + y? = 1. If y = 0 then x = +1. So we get two more points to check: 
(1, 0) and (—1, 0). : 
1 


Subcase B. y — 2. Plugging this into x?-- y? = 1 gives us x? + zal. 


Solving for x then gives x = +, Thus we get two more points: (8S, $) and 
[v64 ij 
99^ 


We now have six points at which to look for maxima and minima of f. We do 
this simply by plugging them in: 





f(0,1) 23 T0, d = -3 
V64 4 454 
qe dee (e 9 5) 56 


So the maximum value of f on 0D is at (0, 1) and the minimum value is at 
(0, —1). 


5 


Problem 69 Let f(x, y) = x? + y? + 2y - 1 and D = {(x, y)x? + 4 < 1). 


1. Find the locations of all critical points of f (x, y) on D. 
2. Find the absolute minimum of f (x, y) on D. 
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7.7 Determinants 


In the previous chapter we saw that the dot product was a way to take two vectors 
and produce a meaningful number. In two dimensions there is another such way. 
Suppose V = (a,b) and W = (c,d) are vectors and consider the parallelogram 
formed by two parallel copies of each. 





We wish to determine the area of this parallelogram. Notice in the figure below 
that there are four triangles and two small rectangles inside a larger rectangle. The 
area of the desired parallelogram is found by subtracting these areas. 


4 


b+d 














The area of the large parallelogram is (a + c)(b + d). The area of each small 
rectangle is bc. The area of two of the triangles is jab and the other two is bcd : 
Thus the desired area is 

1 1 
Area = (a + c)(b + d) — 2bc — 2 ab —2 5d 
— ab +ad + bc 4- cd — 2bc — ab — cd 
— ad — bc 


We often organize the information contained in a set of vectors by writing them 
in a matrix. The area of the parallelogram between them is then the determinant of 
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the matrix. We write this as follows: 


a b 
c d 





|= ad — be 


Notice that there is a bit more information in the determinant than just the area, 
since our answer may be negative. The sign of the answer tells us the order of 
V and W. If we go counterclockwise to get from V to W then the sign of the 
determinant is positive. Otherwise it is negative. For this reason we call the value 
of the determinant the signed area of the parallelogram spanned by V and W. 


EXAMPLE 7-19 
We compute the (signed) area of the parallelogram spanned by the vectors (1, 2) 
and (4, 3). 
1 2 
=1-4-2-3=-2 
3 4 


In three dimensions there is also a determinant, but it is slightly more complicated. 
If we take four parallel copies of three vectors then they span a 3-dimensional figure 
called a parallelepiped. You should think of this as a cube that has been stretched, 
skewed, and rotated. The determinant of the matrix that is comprised of these vectors 
is then the volume of this figure. 


V 


A matrix containing the vectors U = (a, b, c), V = (d,e, f),andW = (g, h,i) 
would be written like this 
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To write the determinant of this matrix it is helpful to copy the entire matrix next 
to the original one. Then add up all the "right-slanting" products and subtract all 
the “left-slanting” products. 


bcr = aei + bf g + edh — af h — bdi — ceg 
g ONE 


Once again this answer may be negative, indicating the way in which U, V, 
and W are situated. The sign is positive if, when we sweep the fingers of our right 
hand from U to V, our thumb points in the direction of W. Once again, this is the 
"right-hand rule;" that we encountered earlier. 

A second way to remember how to compute the determinant of a 3 x 3 matrix 
is to reduce it to computing the determinants of three 2 x 2 matrices: 





d e 
g h 


f 


& f 
g i 


h i 











+e] 





o & 8 
=% m 
-— n 
I 
a 


Notice that each term on the left is preceded by one element from the top row 
of the matrix, with alternating signs. The 2 x 2 matrix in each term is the one you 
get from the original matrix by deleting the corresponding row and column. So, for 
example, the first term is preceded by a, since that is the first element of the top 
row. The 2 x 2 matrix in the first term is the one you get by eliminating the row 
and column containing a. 


EXAMPLE 7-20 
We compute the (signed) volume of the parallelepiped spanned by (1,2,3), 
(1, —1, 1), and (0, 2, 2). 


12 
1 -1 1|21(-2-2)-22-0)432—0) 24-4462 —2 
0 


w 


N 
N 


Problem 70 Find the (signed) area of the parallelogram spanned by 
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Problem 71 Suppose V and W are parallel vectors in R2. Let M be the matrix 
whose rows are V and W. Show that the determinant of M is zero. 


Problem 72 Calculate the (signed) area of the parallelepiped spanned by 


I. (1,2, 3, (1,0, 2), and (—2, 2, —3) 
2. (0,1,3), (CL, 2, 1), and (2, 0, —1) 


7.8 The Cross Product 


A challenging problem is to find the area of a parallelogram defined by two vectors 
in three dimensions. If the vectors are V — (a, b, c) and W — (d,e, f) then this 
area is given by 





Area(V, W) = y (bf — ec)? + (cd — af)? + (ae — bay? 


The general form of this equation looks a lot like the Pythagorean Theorem: some 
quantity is the square root of the sum of squares. We can employ this observation 
to make the above formula a bit more compact: 


Area(V, W) = |(bf — ec, cd — af, ae — bd)| 


Note that we started with two vectors, V and W, and ended up with a third 
that was somehow related. This gives us a new Kind of product, called the cross 
product: 


V x W = (bf —ec,cd — af, ae — bd) 


Recall that the dot product was a way to take two vectors and multiply them, with 
the result being a scalar. The cross product, by comparison, is a way to multiply 
and get a third vector. 

The astute reader will notice that each component of the cross product is the 
determinant of a 2 x 2 matrix. We may thus write its formula like this: 


| 


b c 
e f 


a b 
d e 


C a 


f d 


, , 


vx w =| 




















92 Advanced Calculus Demystified 


It is traditional to define the vectors i = (1, 0, 0), j = (0, 1, 0), and k = (0,0, 1). 
Using this we may rewrite the definition of the cross product as 
a b 
d e 


e „Lla c 
= 
f d f 
We introduced a negative sign in the second term and switched the order of the 
columns of the second matrix. Now our formula looks just like a determinant: 


„|b 
V xW =i +k 




















i j k 
VxW-a b c 
d e f 


This last formula is probably the easiest way to remember how to compute the 
cross product of two vectors. 


EXAMPLE 7-21 
We compute the area of the parallelogram spanned by the vectors (1, 1, 0) and 
(1, 2, 3). To do this we first compute the cross product: 


i j 
(1,1,0) x (1,2,3) 941. 1 
1 2 


v ow 


=i- 0) -j8 -0)+k2- 1) 
=3i—3j+k 
Edu) 


The desired area is now the magnitude of this vector: 





\(3, 23, | = V32 + (3? + I2 = V19 
We have already seen that the magnitude of the cross product has some geometric 
significance. What about its direction? To get a clue we compute V - (V x W). 
Suppose V — (a, b, c) and W — (d,e, f). Then, we compute: 
V-(V x W) = (a,b,c) (bf —ec,cd — af, ae — bd) 
= a(bf — ec) -- b(cd — af) + c(ae — bd) 
= abf — aec + bcd — baf + cae — cbd 
=0 
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From this we may conclude that V x W is perpendicular to V. A similar compu- 
tation shows W - (V x W) = 0, and hence V x W is perpendicular to W as well. 
We conclude V x W is perpendicular to the plane containing V and W. 


EXAMPLE 7-22 

We find a unit vector that is perpendicular to both (1, 1, 0) and (1, 2, 3). From the 
previous problem we know that the cross product of these two vectors is the vector 
(3, —3, 1), and that the magnitude of this vector is V19. To get a unit vector which 
points in the same direction we just divide by the magnitude: 








(3,-3,1) [3419 —3419 v19 
X 19 "7 19 ' 19 


Vv 19 
Problem 73 Find the area of the parallelogram spanned by the vectors 


1. (1,2,3) and (—1,0, 1) 
2. (1,1,0) and (1,0, 1) 
3. (1,2,3) and (3, 1, 2) 


Problem 74 Find a unit vector which is perpendicular to each of the following 
pairs of vectors. 


Problem 75 Show that V x W = —W x V. 


Problem 76 Let U, V, and W be three vectors in IR2. Let M be the 3 x 3 ma- 
trix whose rows are these vectors. Show that the determinant of M is equal 
to U - (V xW). 


Problem 77 Show that the area of the parallelogram spanned by V and W is equal 
to |V||W|sin6, where 0 is the angle between V and W. 
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Quiz 
Problem 78 
1. Let f (x, y) be the following function: 
fæ, y) =x? + 3xy 
a. Find the largest slope of any tangent line to the graph of f (x, y) at the 
point (1, 1). 


b. Find the critical point(s). 


c. Compute the value of 





Pf d 
0x? dxdy 
oy GE 
ðyðx dy? 


for each critical point found. 
d. Does the graph have a max, min, or saddle at the critical point(s)? 
2. Let V = (1,2, 3) and W = (1, 1, 1). Find 
a. a unit vector which points in the same direction as V. 
b. the cosine of the angle between V and W. 


c. a unit vector which is perpendicular to both V and W. 


CHAPTER 8 





Calculus with 
Parameterizations 


8.1 Differentiating Parameterizations 


Suppose $(t) = (f(t), g(t)) is a parameterization of some curve in R*. What 
should it mean to differentiate o at t = tọ? Obviously, the answer will involve 
two numbers, f'(ty) and g’(to). The best way to make sense out of this pair of 
numbers is as a vector: 


d — / / 
1,209 = (f (t9), 8’ (to)) 


If this vector is drawn with its base at the point $ (tọ) then it is tangent to the 
parameterized curve. In addition, if t really represents time, and ¢ (t) the coordinates 
of a particle, then the length of this vector tells how fast the particle is moving at 
time fo. 


Copyright © 2007 by The McGraw-Hill Companies. Click here for terms of use. 
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EXAMPLE 8-1 

At time t the coordinates of a car are 9 (t) = (t°, t°). We wish to determine how 
fast, and in what direction, the car is heading at time t = 2. First, we differentiate 
the parameterization: 


P'E) = Qt, 3) 
Thus, 
(2) = (4, 12) 


The direction of the car is thus the direction that this vector is pointing. Its speed is 
given by the magnitude of this vector: 


|(4, 12)| = V4? + 122 = AV 10 


Nothing significant changes in three dimensions, as in our next example. 


EXAMPLE 8-2 
We find the speed of a particle moving along the helix $ (t) = (cos t, sin f, t). First, 
we find the tangent vector: 


p(t) = (— sint, cost, 1) 


To find the speed we compute the magnitude of this vector: 





|(— sint, cost, 1)| = /(— sint)? + (cost)? + 12 = V2 


Note that our answer does not depend on f, telling us a particle moving along the 
spiral by the given parameterization travels with constant speed. 


Differentiating parameterizations of surfaces is a bit more complicated, since 
now there are two variables. For example, consider the general parameterization 


Wu, v) = (f (u, v), g(u, v), h(u, v)) 


We can take the partial derivative of V with respect to u, say, simply by taking the 
partial derivatives of f, g, and h with respect to u. Our answer is again a vector: 


dy (2 ag oh) 


~ lau ðu’ ðu 


ou 


Geometrically, — (uo, vo) is a tangent vector to the surface parameterized by Y 
at the point V (uo, vo). Another such vector is given by 2? (uo, vo). 
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EXAMPLE 8-3 
We find two vectors tangent to the graph of z — x? 4- y? at the point (2, 1, 5). This 
surface is parameterized by 


W(x, y) = (x, y, x? + y) 


The desired point is at V (2, 1). To find two tangent vectors we simply take the 
partial derivatives of V and evaluate at (2, 1). 


dd = (1,0,2x) 
ax d *, s X 
and so, 
ow 
=~ (2, 1) = (1, 0, 4) 
Ox 
Similarly, 
ow 
= -— (0, 1, 3y?) 
dy 
and so, 


aw 
— (2, 1) = (0, 1, 3) 
dy 


We conclude (1,0,4) and (0,1,3) are two vectors tangent to the graph of 
z= x? + y? at the point (2, 1, 5). 


EXAMPLE 8-4 
If (0, $) represents your longitude and latitude, then your position on the Earth's 
surface is given by 


V (0, $) = (R cos $ cos 0, R cos o sin 0, R sind) 


where R represents the radius of the Earth (assuming the Earth is a perfect sphere). 
Note the similarity to spherical coordinates. We explore the geometric meaning of 
the derivatives of this parameterization. 

First, we differentiate with respect to $ 


ow 


36 = (—R sin $ cos 0, —R sin $ sin 0, R cos $) 
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The magnitude of this vector is thus 





— | = y (R sin ọ cos0)? + (R sing sin0)? + (R cos $)? 


= 
dp 





= V Resin? 6 + R? cos? ġ 
=R 


Notice that the result is constant, i.e., it does not depend on your longitude or 
latitude. The significance of this is that if you change your latitude by roughly 
1 radian then you change your position on the Earth by a constant amount. We will 
see in a moment that this is very different than changing your longitude. Let’s try 
differentiating with respect to 0. 


ow 
ue (—R cos à sin 8, R cos $ cos 0, 0) 


The magnitude of this vector is thus 





ow 7 
E = J (R cos $ sin 0)? + (R cos $ cos0)? = R cos $ 


This says that if you change your longitude by 1 radian, then the actual change 
in distance on the Earth's surface is proportional to the cosine of your latitude. 
Now let's take the cross product of the vectors we have found: 


i j k 
ov ow 
36 ap 7 —RcosósinO | Rcos$ cos0 0 
—Rsin$cos0 -—Rsin$sinO Rcosó 


= (R? cos? $ cos 0, R? cos? $ sin@, R? cos @ sin $) 
The magnitude of this vector is 


ow ow 
36 x a = /(R? cos? à cos 0)? + (R? cos? ¢ sin 0)? + (R2 cos ¢ sin $)? 





=R? cos 


This says that if you stake out a parcel of land that is 1 radian of longitude by 
1 radian of latitude, then its area is proportional to the cosine of your latitude. 
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Problem 79 Find a function s(t) which gives the speed of a particle at time t if its 
coordinates at time t are given by (sin t?, cos t°). 


Problem 80 A curve in R? is parameterized by c(t) = (P, tt, Pj. Finda unit vector 


which is tangent to this curve at the point (1, 1, 1). 


Problem 81 Let S be the paraboloid parameterized by 


V (r, 0) = (r cos0,r sin6, r?) 


O<r<2, 0x0-x2nx 


1. Find two (nonparallel) vectors tangent to S at the point Y (1, Z). 


2. Find a vector which is perpendicular to S at the point Y (1, z). (Hint: Use 
the cross product.) 


Problem 82 Let Y(t) = (cost, sint, t). 


I. Find a parameterization ® of the curve parameterized by WV in which 


2. We may think of i as itself a parameterized curve. For every value of 


_ as being based at the origin, and pointing to 


some point in RÌ. We may thus differentiate this parameterization to obtain 


t we picture the vector 


2 : . 2 . 
another vector v. Show that a is perpendicular to LE. (Hint: Compute 
dio 


"gr and show that the dot product with e is zero.) 


: 2 . . f 
3. Find x so that — = KN, where N is a unit vector. (The number x is called 
the curvature of the parameterized curve and N is called the normal.) 


4. Compute B — ue x N and show |B| — 1. (The vector B is called the 
binormal.) 





5. Compute dB 
6. Find t such that æ = —TN. (The number t is called the torsion.) 
7. Show that 

dN d® 


ee m s B 
dt "a 


1 
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8.2 Arc Length 


Suppose C is a curve, in IR? or R3, which is parameterized by the function Y(t), 
where a < t < b. We would like to find the length of C. We do this by sampling 
C at various points, looking at the length of the line segment connecting adjacent 
points, and adding these lengths. The result will be an approximation of the length 
of C. This approximation will get better and better as we sample more and more 
points along C. In the resulting limit, we obtain an integral which represents the 
length. 
The steps are very straightforward. 


1. Choose n points in the interval [a, b]. We will call these points [1;]7. . 


2. Note that V (t;) is a point of C. For each i we connect V (f;,4) to V(t;) by a 
line segment to approximate C. 


3. The length of each line segment is precisely the length of the vector V (t;+1) — 
W(t). 





4. We now add up the lengths of all the approximating line segments to get an 
approximation for the length of C: 


vau) — vq 
i=0 
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A little algebraic trick will be useful in the next step: 


n n At 
2 IV (t1) — YEDI = 2. MG) — VGL - 


i=0 


W(ti41) — VG) a 
At 








5. Finally, we take the limit of this quantity as n — oo. Note that we can also 
view this as a limit as At — 0. But, by definition, 


. W(t) — YG) dv 
lm ————— = — 
At>0 At dt 


Hence, the desired limit is 


lim J 
n— oo 4 


i=0 


Vai) — YG)! 
At 








EXAMPLE 8-5 
Let C be the curve parameterized by 


W(t) = (tcost,tsint), O<t<1 


Then the length of C is given by 


1 
[fa 
t 
0 





1 
a= [ eost — sinr, sint +r eos) dt 
0 





= | eos — t sint)? + (sint + t cost)? dt 


1 
[ 82 
0 


1 1 i 
VIER + Ine + in| 
0 


Ae od 
=~ + zl + v2) 
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Problem 83 Calculate the length of the curve in R? parameterized by 
W(t) = (cost, sinf, t) 
0<zż<1 
Problem 84 Show that the curve in R? parameterized by 


zo o 
W(t) = | cost, sint, 3! 


O0zrtxl 
has the same length as the curve in R? parameterized by 


(tcost,tsint), O<t<1l 


8.3 Line Integrals 


Let C represent some parameterized curve in R?. Let f (x, y) be a function on R?. 
In this section we make some sense of the integral of f(x, y) over C. First, let’s 
take a step backward and recall the steps used to define the integral of a function 
of one variable, f (x), over an interval [a, b]. 


. Choose n points in [a, b], which we denote as {x;}. 
. Let Ax; = Xj41 — Xi. 
. For each i compute f (x;)Ax;. 


. Sum over all i. 
b 


. Define f f (x) dx to be the limit of the sum from the previous step as n — oo. 


Un A U Ne 


The trick to defining the integral of f(x, y) over a parameterized curve C is to 
follow the above steps as closely as possible. As we will see presently, the most 
difficult step is probably the one that seems most innocent—the second. 


1. Choose n points along the curve C, which we denote as {p;}. (Note that p; 
is a label for a point of IR’, and as such will have two coordinates.) 


. Let dp; denote the length, along C, from p; to p;+1. 
. For each i compute f (p;)dp;. 


. Sum over all i. 


nA A U N 


. Define f f(x, y) ds to be the limit of the sum from the previous step as 


C 
n — oo. 
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The resulting integral is called the line integral of f(x, y). Unfortunately, in 
practice the above steps are impossible to carry out. To make things easier we 
employ a parameterization of C to translate the problem to an integral of a function 
of one variable over an interval of R!. To this end, let V (f), where a < t < b, bea 
parameterization of C. We now repeat the above steps as closely as possible, this 
time using V. 


1. Choose n points in the interval [a, b]. We denote these points as {t;}. The set 
(V (t;)) is then a collection of n points on C. 


2. Let At; = t;4; — tj. Let ôt; denote the length, along C, from {W(t;)} to 








{Y (ti41)}. 
3. For each i compute f (V (t;))ót;. A little algebraic trickery yields 
At; 
fU) Sti = f (QV (ti))ôti 
At; 
= f) Ar 
am i At; i 


4. Sum over all i. 
5. The integral f f (x, y) ds is defined to be the limit of the sum from the 
C 


previous step as n — oo. But as n — oo it also follows that At; — 0. Our 
integrand contains the term A As At; — oo this converges to [2 |. Hence, 
our integral has become 


f aw 
[ roe f rv F] dt 
C a 


Note that if f (x, y) = 1 then the line integral of f over C gives precisely the 
length of C. This gives us one way to think about line integrals. Imagine a vast 
plane. At every point of the plane there is a number which represents how fast a 
point moving through that point will travel. If this number is always one, then a 
point moving along in the plane will always travel with unit speed. Its total travel 
time will then be exactly the same as the total distance it travels. But if the number 
at each point isn't always one then its travel time may be different than the total 
distance traveled. The line integral gives the total travel time. 


EXAMPLE 8-6 
We integrate the function f(x, y) = y over the top half of the circle of radius 1. 
This curve is parameterized by 


V(r)—(cost,sint), O<t<az 
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We now compute 


J fe. sds= f reos sinn (sins, cos dt 
C 0 


T 
= f sint sins, cos) dt 
0 


m 


= J sin ty sin? t + cos? t dt 
0 


m 


= f sine a 


0 


= — cos t| 


=2 


Problem 85 Let f (x, y) 2 y?. Compute the line integral of f (x, y) over the curve 
C parameterized by 


l5 
Y(t) = 3 , 0<ż<1 


8.4 Surface Area 


One of the most interesting applications of multiple integrals is the computation of 
surface area. Let R be a rectangle in the x y-plane. The problem is to find the area 
of the portion of the surface parameterized by Y (u, v). Let R denote the domain 
of V. To find the area we follow similar steps that we used to derive the formula 
for the volume under f (x, y). 


1. We begin by choosing a grid of points ((u;, v;)} in R. 

2. For each grid point locate the point V (u;, vj) in IR?. Connecting adjacent 
points produces a bunch of parallelograms that piece together to form an 
approximation of the parameterized surface. 
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3. We now compute the area of each parallelogram. 
Observe that each parallelogram is spanned by the vectors 


V, = W(ui44, vj) — Y (ui, vj) 
and 
V, = Y (u;i, vj4i1) — Y (u;i, vj) 
The desired area is thus the magnitude of the cross product of these vectors: 
Area = | V, x V,| 


We now do some algebraic tricks: 








Au Av 
| Vu X V.| = Vu x Vil- 
Au Av 
V, Vy 
— x AuAv 
Au Av 





4. The surface area is the limit, as Au and Av tend toward zero, of the sum of 
this quantity over all i and j: 


Au ^v 








V, V, 
S.A.= li 
m e | Au ý Av 
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W(uj41,0j)—V (uj vj) 
oui Au 
converges on 77. Similarly for Av — 0. In the end our summation becomes 


the integral: 
sa.- ff 
R 


But notice that as Au tends toward zero the quantity i z 





du dv 





ow ow 
— x <a 
ou 0 


EXAMPLE 8-7 

We will use the surface area formula derived above to show that the area of a sphere 
of radius R is 47 R*. The most natural parameterization of the sphere comes, of 
course, from spherical coordinates: 


V(0, $) = (Rsin$cos0, Rsing sind, R cos d) 


0x0-x2n, 0xo-m 


The partial derivatives of this parameterization are 


ow 

30 ^ (—R sing sin 9, R sin $ cos 0, 0) 

ow 
a = (R cos $ cos0, R cos ¢ sind, —R sin $) 


The cross product of these two vectors is thus 


i j k 
—|-—RsinésinO Rsin$cos0 0 
Rcos$cos0  RcosóosinO —Rsinó 


ov əy 

— x —Q 

90 0o 
= —R? sin? $ cos 0i — R? sin? 6 sin6j 

— (R? sin $ cos $ sin? 0 + R? sing cos $ cos? 0)k 


= —R? sin? 6 cos6i — R? sin? ¢ sin 0j — R? sing cos ok 
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The magnitude of this vector is 


aw aw 
— x — 
00 ð$ 





= (CR sin? $ cos 0)? + (— R? sin? ¢ sin 0)? + (— R? sin $ cos 9)? 





= JR sin $ cos? 0 + R4 sinf ¢ sin? 0 + R4 sin? $ cos? $ 





= R? 'sin* $ + sin? $ cos? $ 


= R?sin$ sin? $ + cos? $ 


= R?^sinó 


The desired surface area is thus computed as follows: 


m 2x 
ow OW 
S.A. = —— x —| d0d 
| [ [sg Sa 440 
0 0 
m 20 


= | [ sno dodo 
0 0 


T 


= 2R? | sing do 


0 
= 27R* (— cos $)lo 
= 27R7(1+ 1) 
= 47R? 


There are lots of special cases where one may want to compute surface area. In 
each of these cases we can derive an appropriate formula. The first case we consider 
is the surface area of the portion of the graph of z = f(x, y) that lies above some 
domain R. Such a graph can be parameterized as follows: 


Yx, y) = (x, y, f, y) 


1 
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The partial derivatives of this parameterization are 


aw 

m lt 

Ox Ox 

aw 

ay (o E A 
oy dy 


The cross product of these vectors is 


ijk 
— xs =|I o E 
"E ort 
üf.. Əf. 

zda Ek 

a y 


The magnitude of this vector is 


aw awl HN (fN 
EG) 


The desired formula for the surface area of a graph is thus 


"e 
- [4G +B) o #0 








ov aw 
— X - dx dy 
Ox dy 








EXAMPLE 8-8 
We will use the surface area formula derived above for the graph of a function to 
once again show that the area of a sphere of radius r is 4z?. An equation for the top 
half of the sphere is z = yr? — x? — y?. We will compute the area of this surface 
and double it to get the total area of the sphere. 

The tricky part of this computation is the limits of integration. Notice that the 
domain of integration is the region inside a circle of radius r in the x y-plane. For a 
fixed x the value of y can range from —4/r? — x? to Vr? — x2. Hence, we may set 





CHAPTER 8 Calculus with Parameterizations 109 


up a double integral to compute surface area as follows: 


2 2 
b) 
j T AG QD eae 
dy 
BE ae) 
r Vr2—x?2 2 2 
TT Gaal 
y2 — x2 — y? 2 — x2 — y2 
Sn EE 


Vex? 
NE +y? +r? x? o y? 

















= dy dx 
r2—x2— y? y 
= JVr2—x2 
r Jae 
$ 
j J ee dy dx 
=F r2—x2 i: =» 
r " ae 
ET ER a 
r? — x? || fae 
=F 


r [so — sin ! (—1) dx 


=r 


r 


=r f max 


i 
eaux. 


= 2nr? 


As this is the area of just the top half of the sphere we double it to obtain 4zr?, 
the correct formula. 


Another special case of a surface area formula is for a surface of revolution. This 
formula is typically covered in a second semester calculus class. We will show here 
how this formula can be derived from the one above. 

Suppose we begin with the graph of a function z = f(x), whereO <a € x < b, 
and get a surface by revolving it around the z-axis. The resulting surface is best 


110 Advanced Calculus Demystified 


parameterized with cylindrical coordinates, since it has the equation z = f(r). 


Wir, 0) = (r cos0, rsin0, f(r)) 


a<xr<b, 0x0-x2m 


The partial derivatives are 


ow 
—— = (cos, sin 0, df 
or dr 


ow 
aa (=r sin @, r cos 0, 0) 


The cross product of these vectors is 


i j k 
ow ƏY 
a ae cos 0 sin 0 g 


—r sin rcos0 0 


d d 
= ae e rsino + (r cos? 0 +r sin? 8)k 
r r 


d d 
= apie i= r sine; +rk 
dr dr 


The magnitude of this vector is 


aw aw dfi? dfi? 
— x —|= € 2 + es +r? 
or 00 dr dr 
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The desired surface area is now computed as an integral: 


b 2x 


panji 


0 


b 27 d 2 
=f f) +1 d0 dr 
dr 
a 0 
b 2 
=o fJ (2) +1 dr 
dr 
EXAMPLE 8-9 


We compute the surface area of a sphere one final time. This time we think of the 
sphere as a surface of revolution and use the formula just derived. As in the previous 
example we will compute the area of the top half and double it. 

To think of the top half of the sphere as a surface of revolution we begin with a 
function which gives part of a circle: 





ow ov 
— x —| dédr 
or 00 


Applying the above formula for the surface area of a surface of revolution gives 
the following integral: 











1 
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R 
2x 
aR | ————— dx 
R? — x2 
0 


0 
= -zR f u du (where u = Re = x?) 
R2 


110 
—2nzRu? 
R2 





= 2x R? 


As this is just the top half of the sphere we must double it to get the correct 
formula, 4x R?. 


Problem 86 Find the area of the portion of the plane z = 2x + 3y that lies above 
the square with vertices at (0, 0, 0), (1, 0,0), (0, 1, 0), and (1, 1, 0). 


Problem 87 Find the area of the portion of the graph of z = yx? + y? that lies 
above the domain 


I R-((x,y)OzxxzxlO0ozy»czl 
2. D2(G, y)? +y x1) 





Problem 88 Derive a formula for the graph of the spherical equation p = 


f (9). 


Problem 89 Find the area of the surface pictured below, parameterized by 


Wir, 0) = (r cos0, r sin0,0) 


-l<r<l, 0x0-xm 
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8.5 Surface Integrals 


Just as we could integrate a function over a parameterized curve, we can also 
integrate functions over parameterized surfaces. The resulting integrals are called 
surface integrals. Here’s a way to think about surface integrals. Imagine that there 
is a fog permeating all of IR?. As the fog may not be uniform there is some function 
f (x, y, z) which gives the density of the water droplets at each point. Now imagine 
we have some surface in IR? cutting through the fog. We wish to find out the total 
amount of water droplets in the surface. If f(x, y, z) = 1, then this will just be the 
area of the surface. If not, then the answer will be the surface integral of f. 

Recall the relationship between the arc length of a curve C parameterized by a 
function W(t) and the line integral of a function f(x, y) over C: 


Arc Length Line Integral 
b b 
Sarl at f f eva» [Sel dt 
The relationship between the surface area of a surface S parameterized by 


a function V (u, v) and the surface integral of a function f(x,y,z) over S is 
the same: 


Surface Area Surface Integral 
[f |S x S| du dv Sf EOY u, v)) |S x Z| du dv 
R R 


The surface integral of f (x, y, z) over the surface S is denoted as 


[fern ds 
S 
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EXAMPLE 8-10 
Let f (x, y, z) = at E z. We calculate the surface integral of f over the cylinder of 


radius 1, centered on ‘the z-axis, and lying between the plane z = 0 and z = 1. 
First, the surface is parameterized using cylindrical coordinates as follows: 


V (0, z) = (cos 0, sin 0, z) 


0<@0<2n, 0<z<1 


To compute the surface integral we will next need the partial derivatives of this 
parameterization: 


ow 
a8 (— sin6, cos 6, 0) 
ow 
= (0, 0, 1) 
az. 


The cross product of these vectors is 


i j k 
ow ow 
— x — = |—sin cos 0! = (cos@, sin@, 0) 
00 Oz 

0 0 1 


Finally, the magnitude of the cross product is then 


ow ow 
— x — | = Vcos?0 4- sin?0 = 1 


00 Oz 


We are now ready to compute the surface integral: 


fresas f f rivo id e lao d: 


2x 


|= - 1 d dz 
cos? 0 + sin? 0 

0 

2 


m 


IES dz 
0 


Il 
"| eei c" 
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Problem 90 Calculate the surface integral of the function f (x, y,z) =x+y+z 
over the portion of the plane z = x + y that lies above the square in the x y-plane 
where0 < x x land0 < y <1. 


Problem 91 Calculate the surface integral of the function f (x, y, z) — z over the 
top half of the sphere of radius 1. 


8.6 Volume 


In this section, we revisit the calculation of volumes of regions in R°. If the region 
is nice enough, we can just compute a triple integral of the function f(x, y, z) = 1 
over the region. However, if the region is complicated we may want to employ a 
parameterization to help calculate its volume. The steps to define the volume by 
way of integrating a parameterization are similar to those used to define arc length 
and surface area. To this end, we suppose WV (u, v, w) is a parameterization of a 
volume V in IR?, with domain Q. We now derive a formula for the volume of V by 
executing the following steps: 


1. Choose a lattice of points {(u;,v;,wx)} in Q. This defines a lattice 
(V (uj, vj, we)} in V. 

2. Connecting adjacent lattice points breaks up V into a bunch of paral- 
lelepipeds. 


3. We compute the volume of each parallelepiped. The edges are given by the 
vectors 


V, = Visi, vj, wk) — W(ui, vj, Wg) 
V, = Y (ui, vj41, wk) — Y (ui, vj, wx) 
Vw = W(ui, vj, Ug41) — W(ui, vj, wx) 


The desired volume is then given by the determinant of the 3 x 3 matrix 
whose rows are V,,, V,, and V,,. We denote this as 


| V, V, Vi | 


The usual algebraic trickery transforms this as follows: 


AuAvAw 
AuAvAw 
B | V, V, Vy 


[Va Vo V] = | Va Vo Vul 








Au Av Aw 


AuAvAw 





rt 
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4. The desired volume is now the limit, as Au, Av, and Aw tend toward 0, of 
the sum of this quantity over all i, j, and k: 























: V, V, Vy 
Volume — lim >. d ecc d AvAvAw 
Au, Av, Aw 4 m Au ^v Aw 
But, as in the computation of surface area, 
Ve _ Visi, vj, wy) — VU, vj, wy) OW 
lim = lim = 
Au>0 Au Au—0 Au ðu 


Hence, after taking the limit the above summation becomes 


QV IY ow 
volume = | f f 
Q 


Lt. du dv dw 
du Ov ðw 








EXAMPLE 8-11 

We show that the volume bounded by a sphere of radius R is inR? . One way 
to do this is to note that the top half of the sphere can be realized as the graph 
of the equation z = y R? — x? — y?. The desired volume can thus be obtained by 


doubling a double integral: 
1 JR2-x? 


Volume — 2 | | y R? — x? — y? dx dy 


-l. Ry 


However, this integral is difficult to evaluate. 
Instead, we employ a parameterization. The sphere is best parameterized in the 
usual way by using spherical coordinates: 


V (p,0, 0) = (opsin$cos0, p sing sinO, p cos d) 
Usps k, 0<0<2mr, 0S O27 


We now compute all three partials of the parameterization: 
ow : : . 
mem (sin $ cos 0, sing sin 8, cos $) 

p 
ow : : : 
36 = (—p sin $ sind, p sing cos 6, 0) 


ow 
ae = (p cos $ cos 0, pcos $ sin 0, —p sing) 
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The integrand that we will need is the determinant of the 3 x 3 matrix whose 
rows are the above vectors: 


sin $ cos 0 sin $ sind cos à 


Y 3Y 3Y 
E ak =|—psinġsin psingcosd 0 


dp 86 ab 
pcos@cos@ pcosdsind —psing 
— —p?sinó (Check this!) 


We now get the desired volume by integrating: 


N 
a 


T 


Volume — J 


0 


aw aw aw 
— — — | dp d0 dd 
dp 90 ao 


3 
>=> 


Fi 


— p? sing dp dé dọ 


Il 
ea 
omr Sees e 


3 





i ; R 
= J = do dọ 
0 pn 
f f —R^sinó 
= | [Ste do 
0 0 





2xR? f 1 
3 J —singddd 
0 


2xR? 
= cos ol 





—47 R? 
3 





The negative sign is an unfortunate artifact of the fact that the determinant of a 
matrix will change sign if you swap two rows or two columns. We will see in later 
sections that this can be used to our advantage. For now, though, just note that a 
volume must always be positive. Hence, we take the absolute value of the above 
answer to recover the correct formula, ix RÈ. 
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Although this method may seem very complicated at first, observe that it has a 
significant virtue: each step in the integration was very easy. This is the advantage of 
using parameterizations to compute volumes. Any volume problem can be solved 
by a double or triple integral. But employing a parameterization may transform 
difficult integrals into easy ones. 


Problem 92 Show that the volume bounded by a cylinder of radius R and height 
h is zt R?h by integrating with a parameterization from cylindrical coordinates. 


Problem 93 Calculate the volume under the cone z = 1 — x? — y and above the 
xy-plane by employing the parameterization 


V(r,0,z) = ((1 — z)r cos 9, (1 — z)r sin 9, z) 


O<r<l1l, 0x0-x2mx, O<z<l 


8.7 Change of Variables 


In Calculus I and II you learn how to integrate a function of one variable over a 
domain in R!. In Chapter 4, we learned how to integrate a function of two variables 
over a domain in R?, and a function of three variables over a domain in R3: 


f [| 79» 
V 


If the domain V has a convenient parameterization then it may be easier to use 
this than to directly perform the integral. This is very similar to the relationship 
between line integrals and arc length, or surface integrals and surface area. As usual, 
we let V denote a parameterization of V, with domain Q. If V is two dimensional 


we have 
ow aw 
J] re ess ut » s du dv 


And if V is three dimensional 


aw IY ow 
f [ [ 19» etaz- f f [ rwan [RS du dv dw 
du ðv ðw 
V Q 
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EXAMPLE 8-12 
Recall the functions 


t —t t =f 
e+e . e —e 
cosht = 5 , sinht= 








These functions conveniently satisfy 


cosh? t — sinh? t = 1 
d RN d =. 
di cosh f — sinh f, di cosht = sinh t 
The first of these identities means that (cosh f, sinh f) are the coordinates of a 


point on the hyperbola x? — y? = 1. Suppose, we want to find the area of the region 
pictured below. 







(cosh T', sinh T) 
Area=? 


z2- y-1 


P d 


(cosh(- T), sinh(- T)) 


This region is parameterized by 
V (r, t) = (r cosht, r sinht) 
O<r<l1, -T<t<T 


The area of a region R of R? is computed by integrating the function f (x, y) = 1: 


Area= f [ 1 deay 
R 


rt 
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However, to calculate this integral directly would entail finding complicated 
limits of integration. Instead we use the parameterization V. To do this, we will 
need the partial derivatives of Y: 


ow : 

— = (cosht, sinh f) 
or 

ow . 

E = (r sinht,r cosht) 


The determinant of the matrix of partials is thus 


ov ow 
or ot 


_ | cosht sinh f 
^ irsinht rcosht 











= r cosh? t — r sinh? t 


=r 


We now compute the area: 


[pool] 


dr dt 


Ins 


r dr dt 


l 
neu m 
ex 
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EXAMPLE 8-13 
We evaluate the integral 


1 1 y1- 
NI / x? + y? + z dx dy dz 
00 0 


Doing this directly is difficult, so we employ a parameterization. The domain 
of integration is the intersection of a cylinder of height 1 and radius 1, centered 
on the z-axis, with the positive octant. This can be parameterized with cylindrical 
coordinates by 


Wr, 0, z) = (r cos0, rsin6, z) 
Usrsh 0<0<5, Uszrsl 


The derivatives of this parameterization are 


aw 
— = (cos 0, sin, 0) 
or 
aw 
= (—r sin0, r cos 0, 0) 
90 
aw 
= (0, 0, 1) 
az 


To do the integral we will need the determinant of the 3 x 3 matrix whose rows 
are these vectors: 


cos 0 sind 0 
—r sinô rcos0 O|=r 
0 0 1 


We denote the function x? + y? + z in the original integral by f(x, y, z). Then 


m 


1 141-» 1 1 

i QV aw aw 
Tf / x? + y? -- z dx dy dz = | [reo 7 or 96 or dr d0 dz 
0 0 000 


3 1 


J fe cos 0)? + (r sin 0)? + z)r dr dz 
0 0 


Il 
e=- 
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i 


1 
[esr dr d0 dz 
0 


Oe ig 


r^ 4 sre did: 





Il 
Co 
Oe ig 

Ale 


Oe ig 
iE = 
"4 = 

N 
a 
D 
a 
N 


Il 
onL 





Problem 94 Calculate the integral of the function f(x,y) = 2x — y over the 
elliptical region in R? parameterized by 


V (r, 0) = (2r cos 0, r sin 0) 


O<r<il, 0x0-m 





Problem 95 Calculate the integral of the function f (x, y, z) = 
V of R? parameterized by 


i i over the region 


Wir, 0, w) = (r coshocos0, r coshosin0, sinh œ) 


l<r<2, 0x0-x2my, -2<w<2 


(The region V is between two hyperboloids, as depicted below.) 
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Problem 96 Let V be the region inside the cylinder of radius 1, centered around the 
z-axis, and between the planes z = Qand z = 2. Integrate the function f (x, y, z) = 
z over V. 


Problem 97 Let R be the region in the first quadrant of R2, below the line y =x, 
and bounded by x* + y? = 4. Integrate the function 


y? 
fœ») =1+5 
X 


over R. 


Quiz 
Problem 98 
1. a. Find a parameterization for the set of points in R? that satisfies the equa- 
tion x — sin y. 
b. Find a unit tangent vector to this curve at the point (0, 0). 
2. Let Q be the region under the graph of f (x, y) = x? + y? and above the 
square with vertices at (0, 0, 0), (1, 0, 0), (0, 1, 0), and (1, 1, 0). 
a. Find a parameterization for Q of the form W(x,y,t) where 
O<x,y,t<1. 


b. Use the parameterization Y to integrate over Q the function 


TON cg d 


3. Let R be the region of R? parameterized by 
$(r,t) =(rcosht,rsinht), O<r<1,-l<t<l 





Integrate the function f (x, y) = x? — y? over R. 
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CHAPTER 9 





Vector Fields 
and Derivatives 


9.1 Definition 


A vector field is simply a choice of vector for each point. So, for example, a vector 
field on R? would have some vector at the point (1, 2), some other vector at the 
point (—1, 1), etc. We often draw vector fields by picking a few points and drawing 
the vector based at those points. 


Copyright © 2007 by The McGraw-Hill Companies. Click here for terms of use. 
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More formally, a vector field is a function from R? to R*. What goes in to this 
function are the coordinates of the point where you are. What comes out are the 
components of the vector at that point. 


EXAMPLE 9-1 


Consider the vector field (xy, x + y). At the point (1, 1) this vector field con- 
tains the vector (1, 2). At the point (2, 4) it contains the vector (8, 6). If we use a 
computer to draw it we would see something like 





—"———————— LNN BAASI x M PE 
DOE PEINE Aa NM SAA ZA A AZZ 
Joes ete eS NONE DAS ADA AS FO 
ee re T UP aul eae, 
ESD pee SNE uw Y 
Lith ets VOLE PELL LASS 
Nth d ft Ph PP tt TÉ 
bt AS lol JOY qot qt et tat tat fs 
TYTTTTTTTITTI SANAA NAIARA 
VVVNVVNVAVNAVAV AV UG RIS d Ty eR e 
NN NN NN NN NN MT 
‘So 20 6.6 5 X ogg pure fe Se Set eSI DN 
DESSLLILLSTISI, ApsipbnEcétdtd dE 
NNNNNNNNN NV UEDA SS 
NNNNNNNNNN YQ YY ZO TTT TTT 
SiS RN, NVR n onm 1 — eee He He 
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EXAMPLE 9-2 
The vector field (sin x, sin y) is depicted below. 


pdt E E D 
2a277INNNNAANN A 
LFFIVNSNNAN I 
ZELI TANNA AA 
AR£Z£EANN NNNM 
ZATIENNNNNN AT 
2e2rZ7tn~wnnnn il 


MUNN LO om HM 
NNN IZZ” 
CRAKS PLES 
NONUMY pourra 
— — — o o— o c 
————dPdWcw oo 


/ NN NNN 
VIN Se eee Lore E 





o yo o o 


pb ELEME. RN 
46 ———————— Goo 
A A7ZzZPINNNNSN 
VE VE FT NNUS, 
ELE EEAN NONONON 
PSAFAFSFUEVNNNN 
I 727-2441 NN NN 
PAAA NAA 

2 3 4 5 
Yo m Nr d ara 
VON NNNN VL S422 
VNNNNN VY] S424 24 
VIN NNNNA LL 44 4 
VY NNNNN VI] S424 24 
VY NNSNNAN I L427 2 
Se SS, ae 
PIENSO. | one 


Problem 99 Sketch the following vector fields. 


(x, y) 
(x, —y) 
(y, —X) 


9.2 Gradients, Revisited 


We have already seen many vector fields, although we did not use this language. 
Whenever we take a function f and compute its gradient V f at a point we get 
a vector. The set of all such vectors is then a vector field, which we now call 


“grad fe 


EXAMPLE 9-3 
Suppose f(x, y) = xy?. Then the gradient of f(x,y) at the point (x, y) is 
Vf (x, y) = (y?, 2xy). If we draw this vector at various values of x and y, we get 
the picture depicted below. 
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NN NON ON NOS OOO "PE 
Ne NS OS Sc Sa Ss SS Moa S NN 
NN NON €w Cw SSE wu uh ut uat PK LS A. 
id Eis —————— 9 4 4f. 
VVVVNVNNNNNY sA ZZL trite 
AEEA‘ az77rttirtttt 
batata f rr tira dy 
Keb ii piu My 1 2 3 4 5 

/tttatttitttti 
PESECE ECC aa TERES 
LARP RAE SAAS ee HPF ELT ED 
VAN ANANNNAN™ sAr I LLA ttt 
VVVNNNNNNSO i ZI ae a Gy a a 
KON ON AAR NS SS mcn eee ZZ ASS SLT 
No WOO WS NSN SO Ies eom gee et ZZ AI IL A 
VMONONOM M M M o S ooo mer m P" POLO LOL 


EXAMPLE 9-4 
Let f(x, y, z) = xy?z). Then Vf = DE 2xyz?, Saye: 


Problem 100 For each of the following functions, f, compute V f. 


l1. fx, y) 2x t y 
2. f(x,y,z) =x +yz 
3. f(x,y,z) =xy+xz +yz 


9.3 Divergence 


In the previous section, we saw that the gradient operator gives us a way to take a 
function f(x, y) and get a vector field. In this section, we explore a way to take a 
vector field and get a function. Eventually, we will see that the value of this function 
at a point is a measure of how much the vector field is “spreading out" there. 

Suppose (f (x, y, z), g(x, y, z), h(x, y, z)) is a vector field V on R?. Then we 
define the divergence of V, *Div V;" to be the function 


of dg Əh 

Ox oy Oz 
Note that the first term is associated with the first component of V, the second term 
with the second component, and the third term with the third component. This, 


and the fact that the terms are being added, should remind you of the dot product. 
This gives us a purely notational way to remember how to calculate the divergence 
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of a vector field. We let V denote the “vector” (à. Š, i . This, of course, is only 
a vector in a notational sense. But if we suspend our disbelief for a moment and 
allow such absurdities, we can write the formula for the divergence of a vector field 


in a very compact way: 


DivV=V-V 


EXAMPLE 9-5 
Let V be the vector field (x? y, X + yz, xy?z3 js Then the divergence of V is calculated 
as follows: 


DivV=V-V 


"zz (x yx + yz, xyz") 


zy) + = + yz) + s ove) 
= Ixy +z +3xy*2" 
Problem 101 Compute the divergence of the following vector fields. 
(y, z, x) 


1. 
2: (XF Ve — Pee) 
A (x? + y2,x2— yl zl 


Problem 102 Let f (x, y, z) be a function. What is V - V f? 


Problem 103 Show that 


V.(FxG)—-G-.(VxF)-F.(VxG) 


9.4 Curl 


There are many useful ways to apply partial derivatives to vector fields. We have 
already seen that gradients give us a way to take a function and get a vector field. 
Then we saw that divergence is a way to take a vector field and get a function. Here, 
we define a way to use partial derivatives to transform one vector field into another. 


First, recall our notational absurdity, V — 2, ak = . In the previous sec- 
2 B r4 


tion, we defined a new operation by using V in a dot product. Here, we define 
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an operation called curl by using V in a cross product. As before, let V = 
(f(x, y, Zz), g(x, y, 2, h(x, y, z)). Then we define 


curl V = V x V 


(S zT (> Si (2 x) 
=(— -= Ji-|—-—)j+(—-—]k 
dy Oz Ox | Oz Ox oy 


(= dg Of dh dg T) 


dy Oz dz dx’ Ox dy 


EXAMPLE 9-6 
Let V be the vector field i y,x+ yz, aye) Then the curl of V is calculated as 
follows: 


curl V = V x V 
j k 
_ | a a a 
7 | ax dy Oz 





x xtyz xy 


Qxyz y, =z, 1 = x 


Later we will see that the curl of a vector field measures how much it “twists” at 
each point. 


Problem 104 Find the curl of each of the following vector fields. 


1. (x +y,x —z,y +z) 
2. (yz, —xz, xy) 


Problem 105 Let f (x, y, z) bea function and V a vector field on R?. Compute 


1. V x (V f) 
2. V. (V x V) 
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Quiz 
Problem 106 


l. Find the curl of the following vector field: 





imm 
x24 y2? x24 y?’ 
2. Find the divergence of the following vector field: 
(x? 4 y2, y2— x2, 0) 
3. Find the gradient of 
f(x, y, z) =x’ sin(y — z) 


4. Find two vector fields whose curls are (0, 0, z). 
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CHAPTER 10 





Integrating 
Vector Fields 


In the previous chapter we introduced the concept of a vector field, and looked at 
various concepts of differentiation in relation to these objects. In this chapter, we 
look at ways to integrate a vector field. From calculus you should suspect that when 
we combine integration and differentiation something amazing will happen, akin 
to the Fundamental Theorem of Calculus. This is the topic of the next chapter. 


10.1 Line Integrals 


We will always integrate a vector field over a domain that can be parameterized by 
a line, surface, or volume. In this section, we look at the first of these situations, 
integrals over parameterized curves. As in Section 8.3, these will be called line 
integrals. However, the definition of the integral of a vector field over a parameter- 
ized curve is a bit different than the definition of a function over the same curve, so 
there is potential for confusion here. 
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Before giving the definition, it may help to give a physical motivation for line 
integrals. Drop a leaf into a flowing stream and follow it. Eventually, you may see it 
encounter a stick caught on some rocks. If the stick is perpendicular to the flow of 
water, the leaf will get stuck. If the stick is parallel to the flow, then the leaf rushes 
by as the force of the flowing water pushes it along. But if the stick is at some angle 
to the flow the leaf may hit it, and then slowly work its way past. 


SV NX 
IDE 
| de, 
VIA 


| 
| | 


| 
\ \ 

One can explain this by looking at the force the water can exert on the leaf. The 
leaf can only travel in a direction parallel to the stick, which we can represent as a 
vector, L. The water flowing under the stick exerts some force on the leaf, which can 
also be represented as a vector, W. If these vectors are perpendicular then the leaf 
does not move, so the net force F on the leaf must be zero. On the other hand, the leaf 
moves the fastest when L and W are parallel, so F must be greatest in this situation. 
It should be no surprise, then, that F is proportional to the dot product, L - W. 


Before proceeding further, we must dispense with an important technicality. In 
the above figure, there were two choices for the vector L. One of these is depicted 
in the figure, and the other points in the exact opposite direction (but still along the 
stick). If the other choice was made, then the value of L - W would have exactly the 
opposite sign. Which is the correct choice? There is no right answer to this. We just 
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have to declare, before we do any problem, which is correct. Such a declaration is 
called an orientation. We will say more about this idea shortly. For now, just note 
that if the wrong choice was made then we can correct things by changing the sign 
of our final answer. 

Now we suppose that somehow the leaf is confined to some parameterized curve, 
C, as the water rushes past. Sometimes the water may be perpendicular to the curve, 
and sometimes it may be parallel. The problem is to evaluate the total force that 
the water exerts on the leaf as it travels along the curve. (If the leaf always moves 
with unit speed then this is called the work done by the water.) 

Let Y(t) be a parameterization of C. Let W be a vector field that gives the 
direction and speed of the water at each point of the stream. The direction that the 
leaf is moving at the point Y(t) is the vector 2v. Let W(W(t)) be the vector of 
the vector field W that is based at the point Y(t). Then the quantity we want to 
evaluate is 


dY 
[woo ir dt 


This integral is called the line integral of W over C, and is often denoted as 


[was 
C 


Our answer, unfortunately, is still incomplete. Recall that at some point we must 
make achoice, called an orientation. This amounts to deciding if the vector pi points 
the “right way” or the “wrong way.” An orientation can be denoted pictorially just 
by an arrow along C, showing the correct choice. But this does not mean that to do 
the integral, we must choose a parameterization whose derivative points the right 
way. If it does not then we just need to change the sign of our final answer. 


EXAMPLE 10-1 
Let W = (xy, x + y). Let C be the oriented curve depicted below. 
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The orientation is given by the direction of the arrow in the figure. The curve is 
parameterized by 


W(t) = (t?,3t), O<t<2 


At the point V(1) one can easily check that e — (2, 3). This vector points in 
the direction opposite to the orientation. This just means that after we integrate we 
will have to flip the sign of our final answer. 

We now calculate the line integral of W along C using W. First, note that 


WOE) = ((0)G0, t + 3t) = (80, t + 32) 


and 


dY 
dt 


aV i 


2 
f wwo: 5 = fee T +3t)- (2t,3) dt 
0 


fo +3t? + 9t dt 
0 


2 


6 s 3 9 » 
xt t =f 
par as 





0 


2 
=64- 
5 


As the orientation disagrees with V we must now flip the sign, obtaining our 
final answer, —642. 


An easy way to give am Direntanon on a curve parameterized by a function WV (f) 
is to simply declare that 4 ~ points in the right direction. In this case, we may go 
ahead and use W to avana the integral of W without any further considerations. 
In this case, we say the orientation of C is the one induced by Y. 
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EXAMPLE 10-2 
Let C be the curve parameterized by 


W(t) = (cost, t), O<t,< 


with the orientation induced by W. 
Let W = (, 0). We integrate W over C. 
First, note that 


1 


Next, observe 


dY 
— = (— sint, 2t) 
dt 


We now integrate 


1 
[wm feo) sinon dt 
cost 
1 
= f -untar 
0 


= —In|secż||ġ 


= —Inv2 


As the orientation of C is chosen to agree with Y we do not need to worry about 
changing the sign of our answer. 


One may wonder about the connection between line integrals of vector fields 
as defined here and line integrals of functions defined previously. Let f denote a 
function on IR" and C a curve in R” parameterized by W(t). Let W be a vector field 
chosen so that at the point Y(t) of C the vector W(W (t)) is tangent to C (Le., it 


138 Advanced Calculus Demystified 


points in the same direction as at) and has length f (W(ft)). Then we have 


= 
ET 


The definition of the line integral of W over C then gives 


dv 
[weas= [mena 
dt 
G 


e 
Ez 


-J fro» 2E dV 
ks 


WOE) = fOO) 








W(t SEN 
= fr oe S 





dt 


ine dw |’ 





= i fOO) id dt 
E dt 
This last equality is precisely the definition of the line integral of f over C. 


Problem 107 Calculate the line integral of (—y,x) over the circle of radius 1, 
oriented counterclockwise. 


Problem 108 Let W = (x y, x27, y + z). Let C be the curve parameterized by 
Y(t =(t?,t,1-1), O<t<1 
with the induced orientation. Calculate the line integral of W over C. 


Problem 109 Let f(x, y) = xy?. Let C be the portion of the parabola y = x? 


parameterized by 
W(t) =(t,0?), -l<t<2 


with the induced orientation. 
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1. Calculate the integral of the vector field V f over C. 
2. Calculate f (9 (2)) — fŒ (—1)). 


Problem 110 Let C be the subset of the graph of y = x? where0 < x < 1 (oriented 
away from the origin). Let W be the vector field (—x4, xy]. Integrate W over C. 


10.2 Surface Integrals 


In Section 8.3 we defined line integrals of functions, and in the previous section 
of this chapter we defined line integrals of vector fields. By analogy we will now 
define surface integrals of vector fields, just as we had defined surface integrals of 
functions. 


SV NN 
b | 


4— <— 
<— 


The physical motivation for surface integrals of vector fields can again be seen 
by looking at a river of flowing water. Imagine a net placed vertically in the river. 
Suppose you wanted to calculate the amount of water flowing through the net at 
some particular moment in time. Let W denote a vector representing the flow of 
water at some point of the net. Let N be a vector which is perpendicular to the net. 
If the net is parallel to the direction of flow, then no water passes through it. In this 
situation, N is perpendicular to W. The most water passes through the net when the 
net is perpendicular to the flow. In this case, N and W are parallel. It stands to reason, 
then, that the amount of water flowing through the net is proportional to W - N. 
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Once again we encounter here a technical difficulty. Why did we draw the vector 
N in the figure the way we did? If the only condition on N is that it is perpendicular 
to the net, then why couldn't N point in the exact opposite direction? The answer 
is that it can! Which choice is “correct” will again depend on a choice, called an 
orientation. We will discuss this more shortly. 

Now let S be a surface (representing the net), parameterized by Y (u, v). The 
vectors oe and a are both tangent to S. Hence, the vector e x ez is perpen- 
dicular, or normal, to S. Let W be a vector field representing the water flow, and 
W(Y (u, v)) the vector of W at the point V (u, v). Then the total amount of water 


flowing through S is given by the integral 


f [wow v))- (= x =) du dv 


We call this the surface integral of W over S, and denote it by 


[v5 
S 


But this formula isn't the whole story. We still must deal with the issue of 
orientations. Otherwise, two people evaluating f s W : dS may get different answers, 
depending on which parameterization of S they use to evaluate the integral. One 
way to give an orientation is to say which way is “up” at each point of S. This 
can be done by giving a vector O which is perpendicular to S at some point. Our 
parameterization Y "agrees" with the choice of orientation if the vector et x zt 
points in the same direction as O. If our parameterization does not agree with 
the specified orientation then we can remedy things simply by negating our final 


answer. 
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EXAMPLE 10-3 
Let W = (y, x, z). Let S be the portion of a cylinder parameterized by 


WV (0, z) = (cos 0, sin 8, z) 
0s6x7, 0<z<1 
4 
The vector (—1, 0, 0) defines an orientation on S at the point (1, 0, 0). We now 


compute the integral of W over the surface S with this orientation. 
First, we compute the partials, 


aw 
ES = (— sin 0, cos 0, 0) 
aw 
0,0, 1 
Oz =í ! 
And so, 
A : k 
aw aw A nud 
— x — = | — sin cos0 0| = (cos0, sin 0, 0) 





Now notice that the point (1, 0, 0) = ¥ (0, 0). The vector L— x za is equal to 
(1, 0, 0) at this point. This is exactly opposite to the specified orientation of S. 
Hence, if we use the parameterzation V to compute the integral of W over S we 
will have to remember to negate our final answer. 

We now integrate 


jn dS — j [woo z): (55 =) d0 dz 
1 
di 


(sin 0, cos 0, z) - (cos 0, sin 0, 0) dé dz 


oee 


2 sin 0 cos 0 d0 dz 


Il 
M ex 
Oe aja 


=f f sor dé dz 
0 


142 Advanced Calculus Demystified 


1 
J —— cos 20 
2 


0 


p 


x 
4 


dz 
0 





1 


0 
1 
2 


dz 


NI = 


As the parameterization V disagreed with the specified orientation the correct 
answer is — i. 


Problem 111 Let W — (x, y, z). Integrate W over the unit sphere, with orientation 
given by the vector (1, 0, 0) at the point (1, 0, 0). 


Problem 112 Let W — (xz, yz, 0). Integrate W over the surface S paramaterized 
by 


V (u, v) = (u, v, u? + v?) 


Use the induced orientation. 


Problem 113 Let S be the surface given by the following parameterization: 


V (0, $) = (cos $ cos 8, cos ġ sin 0, sing) 


0x0zx2m, -—<@< 


Ala 
Ala 


(Note that this is not quite spherical coordinates.) 
Integrate over S the vector field 


Problem 114 Let W = (xz, yz, 0). Let S denote the intersection of the unit sphere 
with the positive octant. An orientation on S is given by the vector (0, —1, 0) at the 
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point (0, 1, 0). Compute 


[x vas 
s 


Quiz 
Problem 115 
1. Let C be the portion of the graph of x = y? where 0 < x < 1, as depicted in 


the figure. Integrate the vector field (1, 1) over C with the indicated orienta- 
tion. 





2. Let S be the frustum parameterized by 


Wir, 0) = (r cos, r sin, r), 1 <r x2,0x0 < 2z 





Integrate the vector field (+, —1, z) over S with the induced orientation. 
g TT 
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CHAPTER 11 





Integration Theorems 


A centerpiece of the calculus is the “Fundamental Theorem.” This important result 
says that when you take a function, differentiate it, and then integrate the result 
something special happens. In this chapter, we explore what happens when we do the 
same with vector fields. As there are many different ways to both differentiate and 
integrate a vector field this gives us many different theorems. It should be pointed 
out that all of these results (including the Fundamental Theorem of Calculus) are 
special cases of the generalized Stokes’ Theorem, which we will not cover here. 
This approach uses the theory of differential forms, rather than vector calculus, and 
can be found in the book A Geometric Approach to Differential Forms, published 
by Birkhauser, 2006. 


11.1 Path Independence 


The first way we saw vector fields related to derivatives was through the gradient. In 
this section we see what happens when you start with a function f, take its gradient 
to get a vector field, and then perform a line integral along a curve C on this field. 
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Formally, we study integrals of the form 


[v6 


C 


To really understand what is going on we need to go back to the definition of 
a line integral. Suppose C is parameterized by Y(t), and the domain of WV is the 
interval [a, b]. The steps to define the integral of V f over C are then 


1. Choose n evenly spaced points in the interval [a, b]. Call these points {t; }. 
This gives us n point (V (t;} on the curve C. 


. Let At = tj41 — fi. 
. For each i compute V f - q^ at the point V (f;). 


. Sum over all i. 


Un A WN 


. The integral f, c V f : dsis defined to be the limit of the sum from the previous 
step, as n — oc. 


Let us now focus on the quantity computed in Step 3 above. To estimate V f - V 
at a point p we may substitute 


fü) -—JI 4 VADO 


Vf.Vz 
f At 





This approximation gets better and better as At approaches zero. Furthermore, if 
At is small one can also make the approximation 


dY ^ W(ti41) — VG) 
dt tiga — fi 





Putting these together we have 


yf. dY " f QU G+) — f Obi) 
dt At 





Hence, the term Vf - a 7 7. At computed in Step 3 above can be approximated by 
f Qv G4) — f Q6). 


We now shift our attention to Step 4: 


2. vf. T Ar & M f (Wi) - fr) 
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But this last expression telescopes: 
P3 FW Gis) — FG) = FUG) — Fo) = FO) fO a) 


We conclude with the following theorem: 





Suppose V is a parameterization of C with domain [a, b] which agrees 
with the orientation on C. Then for any function f we have 


[vr = f (Y (b)) — f Ql()) 


C 











This should look a lot to you like the Fundamental Theorem of Calculus. There 
is no such convenient name for it. Often it is referred to as the path independence 
of line integrals of gradient fields. 'This is because it says that the result of a line 
integral of a gradient field only depend on the endpoints of the curve, and not the 
path used to get from one endpoint to the other. 


EXAMPLE 11-1 
Suppose f(x, y) — x? + y?. Let C be the top half of the unit circle, oriented 
counterclockwise. We compute f, c Vf - ds. 

All we really have to know is the endpoints of the curve C. The first is (1, 0) and 
the second is (—1, 0). (Which one is which is determined by the orientation on C.) 
Hence, 


[ vitm reso ra =) - = 2 


C 


There is an important application of the independence of path of line integrals of 
gradient fields that may be familiar. Recall that a line integral of a vector field says 
something about how much work you have to do to move an object along a curve C 
in the presence of a force W. Let's say you want to know how much work you have 
to do against the force of gravity to get a heavy package up a mountain. Suppose 
the mountain is represented by the graph of a function f (x, y). That is, at the map 
cooridnates (x, y) the function f (x, y) gives you your altitude. Then the force of 
gravity at the point (x, y) is proportional to the vector V f (x, y) (i.e., the steeper 
the mountain, the harder you have to work to get up it). Let's suppose you have 
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identified your route on a map. You start at the point (xo, yo), follow some curve 
C, and end up at (x1, y1). Then the work you have to do to overcome gravity is 


| v dics YO ye fn od 


C 


Notice that the result is just the difference in elevation between your beginning and 
ending point, and doesn't matter what path you take to get from one to the other! 





EXAMPLE 11-2 


Let W be the vector field (y*z*, 2xyz?, 3xy?z?). We would like to integrate W over 
the curve C parameterized by 


Ha 0, teres 


with the induced orientation. 
This integral can be done directly, but the astute reader will notice that 


W = V(xy?z) 


Hence, we may use the independence of path of line integrals: 


[vr-as= FOD) — f (0) = f(8,4,2) — fC, 1, 1) = 1024 — 1 = 1023 


C 
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Problem 116 Let W = (1, 1, 1). Suppose C is some curve that goes from (1, 0, 1) 
to (1, 1, 1). Calculate the integral of W over C. 


Problem 117 Let C be the curve pictured below. Let W — (sin y, x cos y). Calcu- 
late fe W - ds. 


^ 


y 


tola 





| > 


2 





Problem 118 Suppose f(x, y) = x? — 4x +4+ y? +2y + 1. Let C be a curve 
that starts at (2, —1) and ends at some other point. Show that i? Vf - ds is (strictly) 
larger than zero. 


Problem 119 Suppose C is a closed curve, i.e., one whose beginning and ending 
points are the same. Show that the integral of any gradient field over C is zero. 


Problem 120 Deduce the Fundamental Theorem of Calculus from the indepen- 
dence of path of line integrals of gradient fields. (Hint: Begin by letting C be a 
curve in R! from a to b.) 


Problem 121 7n this section we saw that the work you must do to carry a package 
halfway up a mountain is just proportional to the difference in your starting and 
ending altitudes, and does not depend on the path you take. But it certainly seems 
like if you started at the bottom of a mountain, went to the top, and then came down 
to the halfway point, you'd be doing more work than if you just went straight up to 
the halfway point. Explain. 


11.2 Green's Theorem on Rectangular 
Domains 


The next theorem that relates derivatives, integrals, and vector fields is special to 
R?. Recall that the curl of a vector field 


V= (f(x, yz. 8x, y, z), AG, y, z)) 
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in R? is defined to be the quantity 


i j k 
yayalla D 
ax dy Oz 
f gh 


Although the curl is only defined for vector fields in IR?, we can use it to define 
a natural (and important!) operation on vector fields in R°. First, given a vector 
field W = (f(x, y), g(x, y)) on R? we can easily create a vector field on R3: 


(f(x,y), g(x, y), 0). We will abuse the notation and continue to call this new 
vector field W. We now compute its curl: 





i j k 
VxW- ð a a 
Ox dy Oz 


fo,» g(x,y) 0 
0g X. üf X. 0g Of 
(C) C) GE) 


G-p 
Ox oy 
dg of 


V x W| = —- = 
BM Ox oy 


Hence, 


So, if W is a vector field on R° then the result of |V x W| is a function on R?. 
This function acts much like the derivative of W. We now ask, what happens when 
we integrate the result of such a derivative? That is, if Q is some domain in R2, 
then what can we say about 


ag ð 

J IV x WI dedy f BE LU ed 
Ox oy 

Q Q 





For simplicity, we will assume here that the domain of integration Q is a rect- 
angle. In the next section, we will explore more general domains. 

To understand what is going on here, we will have to go all the way back 
to the definition of a multiple integral. Recall the steps to define the value of 
Jf h(x, y) dx dy: 

Q 
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1. Choose a lattice of evenly spaced points (x;, yj) in Q. 
2. Define 


Ax = Xi X, Ay-—yja-yj 


3. For each i and j compute h(x;, yj) Ax Ay. 
4. Sum over all i and j. 
5. Take the limits of the resulting number as Ax and Ay approach zero. 


To define the value of f f as = a dx dy we follow the above steps, with h(x, y) = 
i ) 


98 _ 3f. Hence, in Step 3 above we compute 


ox dy’ 
0g of 
c po yASA 
(2 je yj) X AY 


If Ax and Ay are small, we may make the following approximations: 


0g a, BOs Yj) — gU. yj) 
ar Ue Xp I 








Ax 
and 
of o fi yj) — FO. yj) 
z-Qi. yj) 
dy Ay 


Hence, the quantity computed in Step 3 can be approximated by 
dg of 
L——— i, YXDAXA 
(£ x je yj) Ax Ay 


n (&GitL V) = 80i y)  fGi yj) — f Gi. yj) dedo 
Ax Ay 





= (¢(xi41, y) — B(x, »)) Ay ^ (fF Gi, yj) — f Ga. yj) Ax 


We now assume Q is the rectangle pictured in Figure 11-1. As in the figure, let 
L, R, B, and T denote the Left, Right, Bottom, and Top sides of this rectangle, 
with the indicated orientations. 

We now move on to Step 4. In this step, we are instructed to compute the sum of 
the values in Step 3. Using our approximation of each term, we can approximate 
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Ym 


Yor 








Figure 11-1 The rectangle, Q 


this sum by 
5o Guess yj) — Gi y Ay — CF Oi, yj) — FO, ys) Ax 
i=0 j=0 
= X Sg Gini, yj) — Gi. yj Ay — yy Uo Yj+1) — f Ga, yj) Ax 
j=0 i=0 i=0 j=0 


= ) (8 @n, yj) — (xo, X Ay — 9 | CF Gis Ym) — f Gs Yo)) Ax 
j=0 


i=0 


= M gn y)Ay — 9 1800 Y Ay — Y. Fis yn) ^x + Y SCi, yo) Ax 
j-0 j-0 i=0 


i=0 


Note that each term in the final expression is a sum involving values of f or g 
at points on the sides of Q. For example, in the first sum we see the expression 
8&(Xn, yj). The number n is a constant, and represents the largest index of the 
x-values. Hence, a point of the form (x,, yj) must lie on R, the Right side of Q. 

We are now prepared to consider Step 5, where we take limits as Ax and Ay go 
to zero. We do this for each term in the above sum. For example, consider the limit 
of the first term: 


m 


m 
li )Ay = li D^ 
alim o2 86 o) y= fim, 2 60s YAY 


Ym 
2(Xn, y) dy 


yo 
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But this integral is the same as the integral of W over R with the orientation 
indicated in Figure 11-1. To see this, first parameterize R by 


V(y)-— Ga, Y) YOLLY L Ym 


Then, 


and so 


J 
Ym 


[w-as= [uo 36 6955) «8:4 
R yo 
Ym 


gx, y) dy 


yo 


The other three terms in the sum from Step 4 give similar integrals over the 
remaining sides of Q. We are now ready for the conclusion: 


at dy Eus p e | ie 
dy 
L T B 


The sides R, L, T, and B of Q, with the orientations as indicated in Figure 11-1, 
when taken together are referred to as its boundary. The usual notation for this is 
9 Q. Using this notation we can write our conclusion much more succinctly: 


[Parma poa 


This final equation is known as "Green's Theorem.” 














EXAMPLE 11-3 
Let Q be the rectangle in the plane with corners at (0, 0), (2, 0), (0, 3), and (2, 3). 
Let W = (xy, x? + y?). We will use Green's Theorem to evaluate the integral of 
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W over 0Q. 


[v [ [rx -x dvds 
dQ Q 

3 2 

= f [59 

0 0 
3 
[ro 
0 
6 


Green’s Theorem enables us to see the geometric significance of the value of 
ôg at at a point (xo, yo). Let Q denote a very small rectangle around this point. 


oe — oF 
Af is roughly constant, and equal to its 


Then ai each point of Q the value of 7 LS — a 
value at (xo, yo). Hence, 


PSP dy ~ (Shen s - 28 a ens ») | «^ 
dy Ox 


Q 


0 0 
= = tay yo) — a yo) | Area(Q) 
Ox oy 
But Green’s Theorem says 


[Parma perma 


Putting these together gives us 


0 0 
(Eo yo) = A ») Area(Q) ~ fw -ds 
X dy 
aQ 
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Or, 
dg of 1 / 
— (X0, — ——(xo, A —————— W . ds 
a 0, Yo) ay 0. Yo) Keno} 
aQ 
So, the function g — — is a measure of the “circulation” of W around each 


point. There is a way to experience such a function for yourself. Many people enjoy 
going “tubing”—floating down a river in an inner tube. Suppose you are tubing and 
decide you want to stop somewhere to enjoy the scenery, so you drop an anchor. 
Then you notice that the water on your left is rushing past you faster than the 
water on your right. What happens? Your inner tube starts to turn. This turning is 


a measure of the strength of E = 
x dy 


_ —_ 


— 
a a rcs 
AJ 
—— —— € 
id xs —— —  —  ——— 


Problem 122 Let Q be the rectangle ((x, y)|O < x x 1,0 < y x 1}. Use Green's 
Theorem to evaluate the integral of (-y*, x?) over ðQ. 





Problem 123 Let C be the curve pictured below. Show that the integral of W — 
(y, x) over C does not depend on b. 
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11.3 Green’s Theorem over More General 
Domains 


Although we have only demonstrated the validity of Green’s Theorem for rectan- 
gular regions, it holds in much more generality. For other shaped regions Q the 
expression 0 Q should generally be interpreted as the “edge” of Q, with a counter- 
clockwise orientation. But there is an important technical restriction on the types 
of regions Q for which this is valid. 

We begin by examining what happens when we look at Green’s Theorem applied 
to two neighboring rectangles. 








Ty To 
Lı Qi Rill Q Ro 
Bı Bə 


We apply Green's Theorem to the rectangle Q; U Q2: 


[5-8 o [8 


Q1UQO» 
x as+ |w. ds 


dQ) 8 Q» 


But the integral of W along R, will cancel with the integral of W along L», so 
the last integral equals 


W -ds 
9(Q1U Q2) 


as Green's Theorem would predict. 
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The cancellation of the integrals over A, and L» illustrates an important phe- 
nomenon. For example, suppose Q is the following region. 


We can use Green's Theorem to analyze the integral of og — ar over Q by break- 


ing it up into rectangles. The boundary of each rectangle then gets an orientation. 
When we add up the integral of W over the boundary of every rectangle, there is 


much cancellation, as in the following figure. 





We end up with the integral of og = a being equal to the integral of W over the 


squares depicted in the following figure, with the indicated orientations. Hence, we 
say that this is the appropriate boundary of Q. 
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In general, we may use the following rule-of-thumb to figure out the orientation 
on each loop of its boundary: 





If Q isa connected region in R° then the “outermost” loop of its boundary is 
oriented counterclockwise and all other loops of the boundary are oriented 
clockwise. 











EXAMPLE 11-4 

We use Green's Theorem to integrate the function x? + y? over the inside of the 
unit circle. If we denote this region as Q, then the boundary of Q is the unit circle 
itself, with a counterclockwise orientation. This can be parameterized in the usual 
way by 


V(t)-—(cost,snt), O<t<2z 


To use Green's Theorem we must find functions f and g so that 


dg əf 
L—uI-a-yziy-(-x) 
Ox oy 


A suitable choice for g(x, y) can be found by integrating y? with respect to x, 
yielding the function xy. Similarly, we may find f (x, y) by integrating —x? with 


respect to y, yielding —x?y. So 


(f(x,y), g(x, y)) = (—x7y, xy?) 
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Finally, we integrate 


[| xe y? dx dy — [ -ds 
Q dQ 

Qn 

= f- cos? t sint, cos t sin? t) . ad dt 

dt 

0 
Qn 

= f- cos? t sint, cos t sin? t) - (— sint, cost) dt 
0 
Qn 

= [55 t sin? t + cos? t sin? t dt 
0 
Qn 

= f 200s? sin’ dt 
0 
Qn 

= f psn? dt 

2 

0 
An 


1 
= f pna 
4 


0 


Li 4d. on 
= — | —u — —sin2u 
4\2 4 





0 


w| a 


Problem 124 Let W = (—y?, x”). Let o be the region in R? parameterized by the 
following: 


d(u,v) = Qu —v,u + v) 


where | < u < 2and 0 < v < 1. Calculate f W -ds. 
ðo 
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Problem 125 Calculate the area enclosed by the unit circle by integrating some 
vector field around its boundary. 


Problem 126 


I. Suppose W = (f(x), g(x)) is a vector field which is defined everywhere 
except at (0, 0). If og — a = 0 then show that the integral of W along every 
circle centered on the origin, oriented counterclockwise, is the same. 

x24 


2. If W = (as, igi) then show that the integral of W along every circle 
centered on the origin, oriented counterclockwise, is the same. 


3. Calculate the integral of ( over the unit circle. 


-y 2a 
x24 y2 , x24 y? 
Problem 127 Let o be the region parameterized by 


o(r,0) = (rcosé,rsind), O<r<1,0<0<zZ 





Suppose W = (x?, e"). 


1. Use Green's Theorem to show that f W . ds = 0. 
do 


2. Let C be the horizontal segment connecting (—1,0) to (1,0). Calculate 
JeW- ds. 

3. Use your previous answers to determine the integral of W over the top half 
of the unit circle (oriented counterclockwise). 


11.4 Stokes' Theorem 


In the previous section, we saw that if W is a vector field in R2, then we can view 
|V x W| as a kind of derivative. When we integrated this "derivative" we saw 
something special happen, namely Green's Theorem: 


f[ vw ddy- [ W-ds 
Q dQ 


We now move our attention to RÌ, and explore a similar phenomenon. Sup- 
pose now W is a vector field in IR?, and S is a surface. Then we wish to explore 
f s(V x W) - dS. A reasonable guess, based on our experience from the previous 
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section, would be 


[vxwas- f was 
S as 


This turns out to be the case, and is called Stokes’ Theorem. We will not prove 
it here, as the proof is extremely similar to that of Green’s Theorem. The strategy, 
once again, is to choose a lattice of points in S. This breaks up S into a bunch 
of parallelograms, and we can chase through the definition of f s(V x W): dS on 
each. At every parallelogram, we see that we get the same as the integral of W over 
the boundary of the parallelogram. But, because of orientation considerations, the 
integrals over the boundaries of neighboring parallelograms cancel. The result is 
the integral of W over the boundary of S. 

One potential complication in using Stokes' Theorem is determining the bound- 
ary of the surface S in question. To get the proper orientation on dS you need 
to know the orientation of S. Recall that this is often given by an outward- 
pointing normal vector, v. To get the orientation on the boundary, we use the 
“right-hand rule.’ To do this point the thumb of your right hand in the direction 
of v. Your fingers will then curl in the sense that determines the orientation on the 
boundary. 
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EXAMPLE 11-5 
Let S denote the top half of the unit sphere, with orientation given by the normal 
vector (1, 0, 0) at the point (1, 0, 0). We use Stokes’ Theorem to integrate the curl 
of the vector field (—y, x, 0) over S. 

First, note that Stokes’ Theorem says that the answer will be the same as the 
integral of (—y, x, 0) around 0S. The boundary of S (with proper orientation) is 
parameterized by 


W(t) = (cost, sint,0), O<t<27 


Thus, we may integrate 


fo x (-y,x))-dS= fixo) -ds 
S as 


Qn 

= f- sint, cost, 0) - (— sinf, cost, 0) dt 
0 
2x 


sja 


0 
= 2r 


EXAMPLE 11-6 
Let S denote the portion of the paraboloid z = 2 — x? — y? that lies above the plane 
z = 1, with an orientation determined by an upward pointing normal. Let W = 
(cos z, sin z, 0). We will use Stokes’ Theorem indirectly to find f s(V x W). dS. 
First, let D be the disk in the plane z = 1 bounded by the unit circle, with 
orientation given by an upward pointing normal. Then 0S = dD. Stokes’ Theorem 
says that the integral of V x W over both D and S is equal to the integral of W 
over 0S. So, to get an answer to the original problem we may evaluate the integral 
of V x W over D instead of S. 
To do the integral, note that 


V x W = (—cosz,sinz, 0) 
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So, on the plane z = 4 we have V x W = (— cos 1, sin 1, 0). A parameterization 
for D is given by 


V(r0)-—(rcosO,rsin0,1, O<r<l, 0x0 -x2mx 
The reader may check that 


aw ov 
a oH 
3 * 39 99 


We now integrate 


1 Ax 


fo x w)-ds= [ fi-cos1,sin1,0) - (0,0, r) dr dé 
D 0 0 
=0 
Problem 128 Let W = (xy, xz, y). Let S be the surface parameterized by 
W(r,0) = (O,rcosé,rsind), O<r<l, 0x0 x2x 
with the induced orientation. Calculate fe x W).dS. 


Problem 129 Suppose W = (f (x, y), g(x, y), 0) and S is a region of R? that lies 
in the xy-plane. Show that Stokes’ Theorem applied to W and S is equivalent to 
Green's Theorem. 


Problem 130 Let S be the portion of the cylinder x? + y? — 1 that lies between 
the planes z = O0 and z = 1, with orientation given by the normal vector (1, 0, 0) at 
the point (1, 0, 0). Let W — (—yz, xz, 0). Calculate the integral of V x W over S. 


Problem 131 Jf W is a vector field defined on all of R?, then show that the 
integral of V x W over the unit sphere is zero. 
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11.5 Geometric Interpretation of Curl 


Just as Green’s Theorem gave us a way to interpret the function 2 — af geometri- 
cally, we can use Stokes' Theorem to give a geometric interpretation of the curl of 
a vector field. Let W = ( f, g, h) be a vector field on IR? and p a point of RÌ. Our 
goal is to understand the meaning of the vector V x W(p). 

Let D be a small, flat disk centered on p. If D is small enough, then V x W is 
roughly constant at every point of D. Let V (u, v) denote a parameterization of D, 


with domain R, in which 


ow ow 


— x =1 
Ou Qv 








In particular, this implies 


ov OW 
Area( D) = Pes x pm du dv 
u y 
R 


aj ae 


— Area(R) 


Since D is part of a plane the vectors normal to D are all parallel. Hence, our 


assumption that the magnitude of 2*. x & is constant implies 
ou ðv 
ow ow 
ðu əv — 


for some fixed unit vector N which is normal to D. 
We now examine the integral of the curl of W over D: 


ov ow 
[vos ff vx wow (Sx =) du dv 
ðu ðv 
D R 


= f | vx Wt) -N duds 
R 


sv x wip ff du dv 
R 
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= (V x W(p))- NArea(R) 
= (V x W(p))- NArea(D) 


Recall that the dot product of two vectors is the product of their magnitudes times 
the cosine of the angle between them. Hence, if we choose D so that N points in 
the same direction as V x W(p) we get 


(V x W(p)):-N =|V x W(p)| 


and hence 


fo x W) -dS ~ |V x W(p)| Area(D) 
D 


To go further we must appeal to Stokes’ Theorem: 


[ovxwy-as= [was 
aD 


D 


Putting these together then gives 


[weas x |V x W(p)| Area(D) 


aD 


or, 


1 
Vx Ww x ————— | W-.d 
IV x W(p)! um | s 
ðD 


Our conclusion is that the magnitude of V x W at the point p is a measure of how 
much W circulates around it. The direction of V x W is the same as the direction 
of N, which was chosen to be perpendicular to the plane in which this circulation 
is greatest. In m seneg, the curl of a vector field is the three-dimensional version 

8 


of the function =~ — 2 which appears in Green’s Theorem. 
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EXAMPLE 11-7 
Let W = (—y, x, 0). Then 


ij k 
ð ə ə 
Vx w=] — — —|]= (0,0, 2) 
Ox dy Oz 
—y x 0 


This is a vector that points up. Notice that a plane perpendicular to an upward 
pointing vector is the plane that contains the "circulation" of W. In contrast, the 
circulation in a plane which contains V x W would be zero. 


oT SESS A + 


foras 
ore -i TS 

NK P e ^e Mw 
ES XEX 
NF) Ne 


—_> 
H 


Problem 132 Among all circles C in R? centered at the origin with radius 1, 
suppose the circulation J c W : ds is greatest when C lies in the xz-plane. Suppose, 
furthermore, that when C is such a loop, 


[v5 
C 


Estimate V x W at the origin. 


Problem 133 Suppose some vector field has the property that the direction of every 
vector is up, and in any vertical plane parallel to the xz-plane the vector field is 
constant. Then what direction would the curl of this vector field point? 


11.6 Gauss' Theorem 


We know that one way to "differentiate" a vector field is to take its curl. It is then not 
surprising that the integral of the curl of a vector field should be special. Another 
way to "differentiate" a vector field in IR? is to take its divergence. In this section, 
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we explore what happens when we integrate the divergence of a vector field. To 
this end, suppose 


W = (fœ, y. z), 8, y. z), hG y, z) 
Recall the definition of divergence: 


ð à oh 
Deve ke 
Ox doy Oz 


The result is a function on IR?. We may thus integrate this function over 


volumes V: 
[ [ [vw 
V 


As in the previous sections, we might guess that there is a relationship between 
this and the integral of W over the boundary of V: 


f [[v wa [was 
y av 


This equality is in fact true, and is known as Gauss’ Divergence Theorem. 
The proof is again similar to the proof of Green’s Theorem. We choose a three- 
dimensional lattice of points in V, and approximate f f f V - W dx dy dz as a sum 

V 














over the points of this lattice. The lattice breaks up V into little cubes, and we find 
that the integral of V - W over each cube is approximately equal to the integral of 
W over the boundary of each cube, with suitable orientations. But faces of cubes 
inherit opposite orientations from neighboring cubes, so in the sum all that is left 
are the faces of the cubes on the boundary of V. 

To properly orient ð V , we simply choose a normal vector that points "out" of V. 


EXAMPLE 11-8 

Let W — (x, y, z). We would like to find the value of the integral of V - W over the 
volume V bounded by the unit sphere. According to Gauss' Theorem, this is equal 
to the integral of W over the unit sphere S. 


rt 
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To evaluate this, we first parameterize the unit sphere in the usual way with 
spherical coordinates: 


V (0, $) = (singcos 0, sing sin 0, cos $) 


0x0-x2n, O0<¢d<z 


Now we compute 
ov ow . : : : . 
36 x P = (— sin $ sin 9, sin $ cos 0, 0) x (cos $ cos 8, cos o sin 0, — sin $) 
= (— sin? $ cos 0, — sin? $ sin 8, — sing cos $) 
To check orientations, note that at V(0, 7) = (1, 0, 0) this vector is equal to 
(1, 0, 0), which does indeed point out of V. Hence, we do not need to worry about 


negating the value of an integral that is computed using Y. 
Finally, we integrate: 


f f [v wem 
V 
= [was 
E 


na 2x 


= | [wing cos0, sing sind, cosg) 
0 0 


x (— sin? $ cos0, — sin? $ sin, — sing cos $) d0 dob 
-T 2x 


= f | ^s o — sing cosg do a6 
0 0 


-T 2x 

= f | -sing do do 
0 0 

=2n | sinodo 


0 
= —4r 
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EXAMPLE 11-9 
Let W — (0,0, £x We integrate V - W over the volume V which is inside the 
cylinder x? + y? = 1, above the plane z = 0, and below the plane z = 1. 

To use Gauss' Theorem we will have to parameterize each "side" of V. First, 
the cylinder, C: 


V(0,z) —(cos0,sin0,z), 0Ox0-x2xz, O<z<l 
Then the bottom, B: 
V (rn0)—(rcosO,rsin0,0, Oxrzl, 0<0<2r 
And finally the top, T: 
V,(r0)-—(rcosO,rsn0,1), O<r<l, 0<0<2r 
Now we must check orientations by computing normal vectors: 


ow ow : . 
—— x —— = (—sin0, cos0, 0) x (0,0, 1) = (cos 0, sin 0, 0) 
90 Oz 


This vector points out of V. 


ov. OW : : 
— — x — = (cosQ0,sin0, 0) x (—rsin0, r cos0,0) = (0,0, r) 
or 90 


This vector points “up,” which is into V. We will have to remember to negate the 
value of any integral that is computed using W. . 
OW, QW, 
x A 
or 00 





= (cos0,sin0,0) x (—rsin0,r cos0,0) = (0,0,r) 


This vector again points “up,” but at the top this is pointing out of V. 
Gauss' Theorem says 


f f [v wat [ was 
V av 


= / W -dS 


CUBUT 


= [was [ w-as+ f w-as 
C B T 


rt 
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We compute each of these integrals individually: 
1 27 


[weas= | [00.2 x (cos 0, sin 0, 0) d0 dz 
C 0 0 
0 


2z 1 


[weas=-f [00.0.0 x (0,0, r) dr dé 
B 0 0 


=0 


2m 1 


[weas= f [ (0.0.1) oor dr ao 
T 0 0 


2n 1 
[re 
0 


d 


Hence, 


|| [wma 
V 


Problem 134 Integrate W = (x?yz, xy?z, xyz?) over the boundary of the unit cube 
in R?. 


Problem 135 Let W = (0,0, e). Calculate the integral of V -W over the ball 
bounded by the unit sphere. 


Problem 136 Let W = (x2, y?, 0). Let V be the region between the cylinders of 
radii 1 and 2 (centered on the z-axis), in the positive octant, and below the plane 


z — 2. Calculate 
J W. dS 
av 
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Problem 137 The surfaces C and D are defined by 


1. C is the graph of the cylindrical equation r = cos 0 in R?, where 0 < 0 < x 
and 0 «x z « I. 


2. D is the set of points in the plane z — 1, which are within i of a unit away 
from the point G, 0, 1). 


Let W be the vector field (0, xyz, 0). Calculate Sep W - dS. 


Problem 138 Suppose V - W = 0, and S, and S» are oriented surfaces with the 
same oriented boundary. Show that 


[weas= [was 
Si S2 


(For simplicity you may assume that Sı and S5 only meet in their boundary.) 


11.7 Geometric Interpretation of Divergence 


Recall that Green's Theorem gave us a geometric interpretation of z = of and 
Stokes’ Theorem gave us a geometric interpretation for V x W. We now use Gauss’ 
Theorem to give a geometric interpretation of V - W. This will be completely 
analogous. 

Let W denote a vector field on R°. Let p be a point in IR? and let B be a small ball 
centered on p. If B is chosen small enough then the function V - W is approximately 


equal to V - W(p) on all of B. Hence, 


|| [v mast vm f [ [ ada 
B B 


= V.W(p)Volume(B) 


Now, Gauss' Theorem says 


|| [vw [was 
B OB 
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Combining these we get 


V -W(p)Volume(B) ~ fw -dS 


0B 


or, 


1 
V. f — s 
Wop) ws" do 
8B 


The quantity f. ap W : dS is a measure of the net amount of W which leaves the 
ball B. So, if the same amount of W enters and leaves B then this integral will 
be zero. One way to think about this is that he W - dS, and hence V - W(p), is a 
measure of how much W “spreads out” at p. 


EXAMPLE 11-10 

The vector fields pictured below have more “leaving” each point than “entering,” 
and so have positive divergence. For example, the second one pictured might be 
something like (x, y, z), whose divergence is 3. The third one might be (0, 0, z), 
whose divergence is 1. Can you think of a vector field like the first one? 
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EXAMPLE 11-11 
The vector fields pictured below have zero divergence. The first one is a constant 
vector field. The second might be something like (y, —x, 0). 


Problem 139 Suppose B is a ball of radius 1 centered around the origin. Let W 
be a vector field and suppose 


f f [yw 
B 


Estimate the value of V - W at the origin. 


Quiz 
Problem 140 


1. Let C be any curve in R3 from (0, 0,0) to (1, 1, 1). Let W be the vector field 
(y?z2, 2xyz?, 2xy?z). Calculate Jo W . ds. 


2. Let o be the region parameterized by the following: 


6G) = Gv wy, 1<u<2, Iz» $2 


fa — lnx) -ds 


ðo 


Calculate 
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. Let C, and C» be curves given by the following parameterizations (with the 


induced orientations): 


Ci: (t) =(,0,0), 2xnztx4x 


C2: Y(t) =(tcost,tsint,0), 2x xt «4x 


Show that for any vector field W such that V x W = (0, 0, 0) the following 


is true: 
[weas= [vw 
C C; 





Final Exam 


Problem 141 


1. Let f (x, y) be the following function: 


a. 


b. 


fœ, y)=xy+x-—2y+4 


Sketch the intersections of f (x, y) with the xz- and yz-planes. 
Find the critical point(s) and compute the value of 








3? f of 
ax? Oxdy 
orf Pf 
dyax  dy2 


Can you say if the graph has a max, min, or saddle at the critical point(s)? 


. Find the slope of the tangent line to the graph of f (x, y), in the direction 


of (1, 2), at the point (0, 1). 


. Find the volume under the graph of f (x, y), and above the rectangle in 


the xy-plane with vertices at (0, 0), (1, 0), (0, 2), and (1, 2). 
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. Parameterize the portion of the graph of z = 4 — x? — y? that lies above the 


xy-plane. 


. Parameterize the curve that lies on a sphere of radius I such that 0 = q. 


4. LetW = ixz^, 0, xz’). Calculate V - W and V x W. 


. Let V be the volume in the first octant, inside the cylinder of radius 1, and 


below the plane z — 1. Integrate the function 


foy 22y/l1I42x?4y' 


over V. 


. Let C be the curve parameterized by the following: 


p(t) = (2cost,2sint,t?), O<t<2 


Integrate the vector field (0, 0, x? + y?) over C. 


. Let F be the vector field (0, —z,0). Let P be the portion of a paraboloid 


parameterized by 


$ (r,0) = (r cos0, r sin0, r?) 


O<r<1, 0<0<5 


Calculate f F. dS. 
P 


. Let W = (yz, xz, xy). Let C be the curve parameterized by 


4t 4t . 4t T 
$(t) —|-—cost, —sint, —), O<t<— 
T T mU 4 


with the induced orientation. Evaluate Je W - ds without integrating! 


. Let S be the can-shaped surface in R? whose side is the cylinder of radius 1 


(centered on the z-axis), and whose top and bottom are in the planes z = 1 
and z = 0. Let W = (0,0, z2). Use Gauss’ Theorem to calculate Js W -dS. 





Answers to Problems 


Chapter 1: Functions of Multiple Variables 
Problem 1 

1. f3,2) 232-229-8217 

2. g(0, 5) = sin(0) + cos(5) =0+0=0 

3. h(2, Z) = 2? sin(Z) = (41) = 4 

Problem 2 

(b), (d) 

Problem 3 
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Problem 4 


^z Ag ty 
































RT 
e 



































180 Advanced Calculus Demystified 


Problem 6 


1. The circles get larger as z increases, but they begin to “bunch up.” 

2. The circles get larger as z gets larger, and their spacing does not change. 

3. The circles get smaller as z gets larger, eventually becoming very close to the 
origin. 

4. At z = 0 the level set consists of infinitely many circles. For |z| > 1 there are 
no level curves. 





Problem 7 

4. 9 5. i | 

Problem 8 

i »€ 
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Problem 9 
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Problem 10 


The function (cos t”, sin t°) describes a point moving clockwise around the circle, 
and slowing down until t = 0. (It then reverses direction when t becomes positive.) 
In contrast, the function (cos t, sint) describes a point moving counterclockwise 
around the circle, at uniform speed, no matter what f is. 


Problem 11 
c(t) = (t, f (t)) 
Problem 12 


1. The curve spirals around the z-axis, rising faster and faster, so that there is more 
and more space between successive coils. 


2. The curve spirals around the z-axis, descending, becoming closer and closer to 
the x y-plane. 


3. The curve spirals up around the cone z = yx? + y?. 


Chapter 1 Quiz 
Problem 13 
1. Left handed 


2.(a) The intersection with the xy-plane looks like the x-axis. The intersection 
with the xz-plane also looks like the x-axis. The intersection with the 
yz-plane is the line z — y, pictured here. 


| 





Answers to Problems 


(b) 











"ES 








Chapter 2: Fundamentals of Advanced Calculus 


Problem 14 


1. Along the y-axis we have 
x2 


lim | ———s 
(x,y) (0,0) x2 + y3 


Along the x-axis 
xj 


lim |———; 
(x.y) (0,0) x2 + y? 


2. Along the y-axis 


k^y 


lim | ——— 
(.y)9 (0,0) x3 + y’ 


Along the line y = x 


xy 


lim | ———-— 
(x,y) > (0,0) x3 + y3 


= liñ = 
(x.y)>(0,0) y3 


x2 


= tlm 5 = 
(x,y)— (0,0) X 


— lim = 
(x, y) (0,0) y3 


x3 1 


lim — =- 
(x.y)>(0,0) 2x3 2 
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3. Along the y-axis 


fai. Fas Pay 


O00) J/x2-E y2 Œy) (0,0) y 


Along the line y = x 


x+y 2x 


lim ————-— lim —_=¥V2 
(x.y) 0,0) J/x2 + y2 60.0 /2x 


4. Along the y-axis 


xy 


lim ———- lm — = 
(.»)2(0,0 x? + y? (x,y) 0,0) y? 


Along the curve y = x? 


xy? 6 





lim = lm ==] 
(x.y) 0,0 x? Hy? (x,y) (0,0) x? + x$ 


Problem 15 


l. y#x 2. y 2 x? 3. The domain is empty. 


Problem 16 


Yes. For small values of t the function sin(t) is approximately equal to t. The value 
of x? + y? is the square of the distance from (x, y) to (0, 0), which is small for 
points near the origin. 


Chapter 2 Quiz 
Problem 17 


1. On the x-axis we know y = 0, and so the function f(x, y) = 0 (as long as x is 


not also 0). Similarly, on the y-axis we know x = 0, and so again the function 
f (x, y) = 0. But when x = y we have 


x sinx 


Pob = oan 
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To evaluate the limit of this as x — 0 we use L’ Hopital's rule twice 








. xsinx . sinx+xcosx 
lim = lim ———— — — 
x0 2x? x0 4x 

2cosx —xsinx 
= lim 
x0 4 
1 
2 


Hence, the limit does not exist. 
D = — ]. Hence, the function 
f (x, y) = 1 for every value of x and y. The function f (x, y) is thus a constant 


function, which is continuous everywhere. 








3. The first problem we may encounter is a zero in the denominator, in which case 
x — y? = 0. The second problem is taking the log of a nonpositive number, in 
which case = < 0. But this second situation implies x — y? < 0, and hence 
includes the first. 

If x — y? < 0, then x < y*, Hence, the domain of f (x, y) is all the points 
(x, y) such that x > y?. 





Chapter 3: Derivatives 
Problem 18 
a 
1. £(2,3) = 4, £(2,3)=2 
3 2013, 0 ue? 
3. 9 Q, 3) = V6 + NS, 9f o, 3) — É 


2? ay 
Problem 19 
L amy, faa 241 Fos 
Problem 20 
ar = —] + y? and af = 2xy — 2y. Setting these both equal to zero gives the system 
of equations 
0=-1+y’ 


0 = 2xy — 2y 
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The first equation tells us y? = 1, or y = +1. If y = 1 then the second equation 
becomes 0 = 2x — 2, so x = 1. If y = —1 then the second equation becomes 
0 = —2x + 2, so x = 1. Hence, the solutions are (1, 1) and (1, —1). 

Problem 21 
First, note that $ (1) = (2, 1). So 


A eq)» - Le, iret dnd 
af af 

— 1) = 2,1) —-1 
a0) = 2.1 


Now, note that x(t) = 2t, so A = 2. Similarly, y(t) = t?, so 2 — 2t and 
2 (1) — 2. Finally, we compute 


_ of of dy 
C Fo) = Temi A 60770) 





=6-2+(-1)-2 
= 10 


Problem 22 


First, note that t = 2xy and a = x?. Now we compute: 





T soo = = A gaye pet QD Q) 
_ of of 
1,3 2) F 1,3)— (2 
à = (2) ay $ (2) 


dx 2 dy 
CADO) GTI: VO) 


= (6(-2) + (D 
=-l1 
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Problem 23 
1. First, note that x(1, 2) = 2 and QUU j= 2 ence, if (u, m — (1,2) then 
(x, y) = (2,5). Next, note that & =v and 2 = 1, and so xq, 2) =2 and 
ay 
> (1, 2) = 1. We now compute: 
Ta, 2) = Ta, 52a, 2) + Ta, 52a, 2) 
= (2)(2) + G)) 
ET 


2. x(2, 1) =2 and y(2, 1) "E: Hence, if (u, v) = Q, 1) then (x, y) = (2,3). 
Next, note that x =u and 2 — 2v, and so & (2, 1) = 2 and 2 x; (2, 1) = 2. We 
now compute: 





Ho, 1)= Te, aže, n+ Ze, 320, 1) 
= (0) + (-D)Q) 
=0 
Problem 24 
1. First, note that x(5, 7) = 5 +7 = — and y(5, 7) = 5 —z = —5. Now, 
3m mx . (n m . 
fau y)-f (=. -i) — sin (F — z) = sin(x) = 0 
2. First, note that af œ, y) = ar 2 so af (am. -5) = = cos(z) = —1. Simi- 
larly, 5 af -(x, y) = ede +y), so 2 vos -5- = cos(x) = —1. Finally, note that 
A = 1 and 2 = 1. Hence, 
of _ dof a Of Oy ,- z AD 
am à yr 1)(1) + (—1)(1) = —2 
3. Note that Es = | and I — —],so 
O) ofa 0 
f Fx .Cpgec-pcp-0 





ð ax dv | dy dv 
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Problem 25 
l. SP o. $i 
&? ° 8y — 
9? 2 
OF o. 35.5 
ax dy dy? 
a £e. SR. 
ax2 — ^" By 
a2 2 
P£oy eE 
ax dy dy? 
a2 32 
= a = —y'sin(xy^, ae = 2y cos(xy?) — 2xy? sin(xy^) 
3? f of 2 2.2 2 
— = 2y cos(xy? y= 2xy? sin(xy? J; = 2x cos(xy^) — Ax^y* sin(xy^) 


dx dy ay? 


Chapter 3 Quiz 
Problem 26 


1. Of = 2xy + 3x7y? and 3 = x? + 2x3y. 


2. FOA) = CFE- DEPE- D 
= tft —1) +f — 2t +1) 
uj ex G =f 


3. First note from Question 1 that V, 1) — 5 and Za, 1) = 3. Now, 


df(w(t)) — af dx | of dy 
dt ax dt Oy dt 
= (55)(3) + 3)() 


= 18 
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4. Since x(u, v) =u? + v, it follows that x(1, 1) = 2. Hence, when (u, v) = 
(1, D), (x, y) = (2, 1). Now, note 





of 
—(2,1)=44+12=16 
x ) + 
à 
9f o, 1) =4+16= 20, and 
dy 
ax 
—(1,1)=2(1) =2 
x ) (1) 
Now, we have 
12 zi 
ou 
LOU OP QU Oy 
|» xu  Oyüu 
dy 
= (16)(2) + 20— 
(16)(2) + m 
And so, 
dy 
12 = 32+ 20— 
T ou 
from which it follows that 
y ou 
ou 
Chapter 4: Integration 
Problem 27 
1 3 T 
1. [ [r+ aay =f e+ xy? dy 
2 2 x=2 
0 2 0 
^9 9 
= | 2+ 2y = 2 = 9" d 
[5*» y ay 
0 
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x2y? dy dx 
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5 
+ = dy 


5 
2.22 


i 











NJA 


dy 
x=0 


sin(x + y) 





E SN 


sin ( 


T 
— COS (F+y) — cosy 


T . 
5 ty) -siny dy 


zu 





2 
0 
7 


— COS T —c057. + cos = + cos 0 


=2 


Answers to Problems 


Problem 28 


0 2 


Volume = NE dx dy 


==] 


0 

4 
=f tes} 
—1 





x=—1 





Problem 29 


1 1 
Volume = ff x"y "y" dx dy 
0 0 
E: ey m 
0 
fe 
0 


1 
= (n-Dm«D 
= 1 
~ (n+1)(m+1) 


1 


dy 
x=0 








1 
m+1 





0 
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Problem 30 


To find the answer just set y = 1 and integrate with respect to x 


a 1 
"m e 





2 
Area = |a L-—g* 
0 








Problem 31 
1 y 1 y 
l. Jf dx dy = [er dy 
0 0 0 a7 
1 
= [ory dy 
0 
1 
= fr dy 
0 
Lal 
xdi 
8' |g 
_ 1 
E 


2x 


2 2 
2. J J e"? dy dx = J eo 
0 0 


x 








a CHE 
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Problem 32 
The graph of y = x? — x — 2 is depicted below. 





ovy 


-15 -1 -0.5 0 0.5 1 1.5 25 
x 





The x-intercepts of the graph are x = —1 and x = 2. Fixing a value of x in between 
these two numbers, y will range from x? — x — 2 to 0. This tells us how to set up 
the integral 


2 0 2 
0 
J J x? dy dx = J x y| a a dx 
—1 x?—x—2 -1 


2 


= Í =x? (x? —x-—2)dx 


zi 


2 
= f -t+ 420 dx 
-1 





"M I 5 I 4 2 3 

= ud Ti is T 
32 16 16 1 

= == Se 
5 4 3 5 4 3 
33 15 18 


574 "3 
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Problem 33 
The region R is depicted below 











* The range of possible x-values is 0 to 1. If we fix a value of x between these 
numbers, then y can vary from x? to x. So the integral can be set up as 


l1 x 
f | 1e 5o 


0 x3 


* The range of possible y-values is also 0 to 1. If we fix a value of y between these 
numbers, then x can vary from y to ./y. Hence, the integral can be set up as 


Jy 


1 
f | 19» 
0 


y 


Problem 34 


The original limits of integration tell us the shape of the region over which we 
are integrating the function sin(y*). The outside limits tell us that the domain of 
integration lies between the vertical lines x = 0 and x = z?. Fixing a value of x 
between these numbers, y can vary from 4/x to x (the inside limits of integration). 
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This tells us the shape of the domain of integration is 














Now notice from the figure that the possible range of y values is from 0 to z. 
Fixing a value of y between these numbers, x can range from 0 to y?. Hence, we 
may rewrite the integral and evaluate as follows: 


y? x 


22) 
J [ o» dx dy — / x sinò) o dy 
0 0 


0 


j y^ sin’) dy 


0 


T 





: O°’) 
— — COS 
ae 


1 3 
=o z^ — cos 0) 


la 3) 
— — COST 
3 


Problem 35 





T 
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The solid is pictured above. To find its volume, note that if we fix x and y then z 
varies from 0 to 1 — y?. If we fix y then x varies from 0 to 1. Finally, the range of 
values that y can take on goes from —1 to 1. Hence, we have the triple integral: 


1-y? 


1 1 1-) 
f f [omm 
-10 0 


dx dy 





Il 
| | | 
—~> o_O. 

oN 

an 

| 

< 

N 

> 

a 

< 





Problem 36 


The volume can be computed with a triple integral as follows: 


1 1-3 1-x—y? y? 


1 1— 
Wi | 1dzdrdy= f [ 1-1-3 dx dy 
-1 0 0 -1 0 

1 2 


1 2 2 iid 
L—jxzcx^ =x 
J p ow 
=j 


1 


Va. iy 
[50 y) dy 


dy 





0 
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Problem 37 
| JSI-y? 1-x2-y? 
J J 1 dz dx dy 
—1 -i-y 0 
Problem 38 


Th —X1-z? —y 1—z?—x? 


Chapter 4 Quiz 
Problem 39 
3 


dy 
2 





2 3 2 
1 
1.(a) J [ 5e + y) dx dy = J 5 sin(2x + y) 
1 2 1 


2 


1 
a | sin(6 + y) — sin(4 + y) dy 
1 


1 
x (cos(4 + y) — cos(6 + y))|? 


1 
= z (eos 6 — cos 7) 


(b) The trick with this integral is to reverse the order of integration: 


] y 


D] i 
f [i9 f [enn 
0- x 0 0 

1 


= | VTF žy 
0 


Ju du 


Il 
Te 
NI = 
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Qi-1) 


2. The domain of integration is a disk of radius 1 in the xy-plane. Hence, if we fix x 
then y can vary from —4/1 — x? to //1 — x2. The biggest and smallest possible 
values of x are —1 and 1. Hence, the desired volume can be computed by the 
integral 





1—x? 


1 
J J V1 — x? — y? dy dx 
zi 


/1—x? 


Chapter 5: Cylindrical and Spherical Coordinates 
Problem 40 


1. (V2, V2, —1) 2. (0,0,3) 3. (2,2/3, 0) 


Problem 41 


1. The xz-plane. 2. The yz-plane. 3. The xy-plane. 
4. A right circular cone centered on the z-axis. 5. A sphere of radius 1. 


Problem 42 





Problem 43 
1. (0,0,0) 2. (—3,0,0) 3. (V6, 6, 2) 
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Problem 44 


1. A single point at the origin. 
2. The xy-plane. 

3. The negative z-axis. 

4 


. A cylinder of radius 2, centered on the z-axis. (Recall that o sin $ = r, so this 
is just the graph of the cylindrical equation r = 2.) 


5. A horizontal plane at height 2. (Recall that p cos $ = z, so this is just the graph 
of the rectangular equation z = 2.) 


Problem 45 


4 


e 
l. 2. lip mm 


Chapter 5 Quiz 


Problem 46 
X J/3 
1. = 0 = — = — 
(a) x =r cos cos 6 5 


iud (wr d 
=r sin = sm- — = 
z 6 2 








z—z-—2 
m m X243 v6 
= 0=2 =2 = 
(b) x = p sin $ cos sin " cos 6 5 5 5 
x v21 4/2 
= 0 —-2 i = 
y = p sin ġ sin uum sin 6 22 5 


z= peos = 2002 = 2 = J/2 


2.(a) In polar coordinates the equation r = cos 6 describes a circle of radius } n 
centered on the point (1 5,0). In cylindrical coordinates, we see the same 
picture no matter what z is. Hence, the graph is the cylinder pictured below. 
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(b) Consider the sphere of radius 1. On this sphere the points where 0 = $ forms 
a curve which starts at the top (where z — 1), ends at the bottom (where 
z = —]), and passes through the point where y = 1. Now we see the same 
picture no matter what ¢ is, so we see this curve on a sphere of every radius, 


as in the figure below. 
4 





Chapter 6: Parameterizations 

Problem 47 

1. y, y) = @, y, x”) 

2. w(0, z) = (8? cos 0, 9? sind, z) 

3. WO, 9) = (0? sing cos 6, 6? sing sin 0, 0? cos $) 
Problem 48 

1.¢=0 2,.y=x4+xz 3. z2=sinr 
Problem 49 

1.0«0«7, O<¢<2 25200, 50m 
Problem 50 


1. A portion of a cone that lies above a square. 
2. A portion of a cone that lies above a circle. 


3. A portion of a cone that lies above half of a circle. 
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Problem 51 
1. Start with the graph of x = y? + z. Now stretch this graph in the y-direction by 
a factor of 2. Finally, translate the entire picture —1 unit in the z-direction. 


2. Start with the graph of z = x + sin y. Now reflect in the plane y = x, switching 
the roles of x and y. 


3. Start with the cylinder r = 1. Stretch in the x-direction by 2 and the y-direction 
by 3, so that cross sections are now ellipses. Now move the central axis 1 unit 
in the y-direction and —1 unit in the z-direction. Note that this last translation 
does nothing to the shape. 


Problem 52 


If the curves shown were circles, with the one corresponding to z = 1 having radius 
1, then the uniform spacing of the level curves would tell us that the shape is a cone. 
Such a cone would be the graph of the cylindrical equation z — r, and would thus 
be parameterized by 


Wir, 0) = (rcos6,rsin6,r) 
O<r<l, 0<0<2xr 
The actual level curves we are given are ellipses. These can be obtained from 


circular level curves by multiplying the x- and y-coordinates by a factor which will 
stretch them the appropriate amount: 


V (r, 0) = (2r cos 0, 3r sind, r) 


O<r<l, 0<0<2r 


Problem 53 
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Problem 54 


1. The region between spheres of radii 1 and 2, and above the xy-plane. 


2. The region inside a cylinder of height 2 and radius 1, between the planes z = 1 
and z = 3. 


Problem 55 


V (r, 0,t) = (rcos6,r sinO, tr) 


O<r<il, 0<0 <27; O<t<l 


Chapter 6 Quiz 
Problem 56 


1. For each point on a cylinder of radius 2, we know r = 2. If such points are also 
on the graph of z = 2r, then z = 2(2) = 4. Hence, the problem is asking for the 
points of a cylinder of radius 2 for which 0 < z < 4. This is most easily done 
by utilizing cylindrical coordinates: 


WO, z) = (2cos0,2sin6, z) 


0x0-zx2m, O0<z<2 


2. If the 2 were not at the x-coordinate then this parameterization would look just 
like spherical coordinates. The restrictions on p, 0, and @ would then specify 
one-quarter of a ball. The factor of 2 then stretches this shape in the x-direction. 
The result is the figure below. 
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Chapter 7: Vectors and Gradients 


Problem 57 
1. V+W =(7,7), —V = (-1, —6), V — W = (-5,5) 
2. V+W = (1,2), -—V = (0,0), V — W = (-1, —2) 


Problem 58 
1. /13 2. V10 
Problem 59 


The magnitude of (2, 1) is 2? + 1? = 4/5. To get a unit vector that points in the 
same direction we just divide by the magnitude: 


1 24/5 48 
qae ed $) 


Problem 60 


The vector (5, 12) lies on a line of slope 2. A perpendicular line will have a slope 
which is the negative reciprocal of this, >. A vector that lies in this line is (—12, 5). 
To get a unit vector that points in the same direction, we divide by the magnitude: 


1 1 12 5 
————24À(-12,5) = —(-12, 5) =(-—, — 
JAAEYSS egg l 13 ib 


Problem 61 


1. 29673) + (90) = -2 
2. (06) + 000) = 14 
3. (—2)(6) + (-DC-3) = —9 








Problem 62 
die [28.35 m V2 
1 cos = = == 
.44C3. | — 204/10 10 
-— —2, —1) - (6, —3) =9 -Ó 
(-2, —1}||(6,—3)| 4/54/45 15 
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Problem 63 


V. W —|VIIW|cos0 = 9cos@. So V - W is largest when the value of cos @ is 
largest. Since the largest possible value of cos 0 is one (when 0 = 0), the largest 
possible value of V - W is 9. Similarly, the smallest possible value for cos 0 is —1, 
so the smallest possible value for V - W is —9. If V and W are perpendicular then 
cos 0 = 0, so the smallest possibility for |V - W| is zero. 


Problem 64 

1. No, V-W —44z0. 2. Yes, V. W =0. 

3. Yes, V. W = 0. 4. No, V -W = —9 Z0. 
Problem 65 

l. e Vf(x, y) = (Iny, 5) 


VEC, D- (È, $) = (0,1) - 3,2) =F 

The direction of the maximum rate of change at (1, 1) is V f (1, 1) = (0, 1). 
This is already a unit vector. 

IVF, | =1 

Such a line is perpendicular to V f (1, 1) = (0, 1). Such a vector is (1, 0). 
Vf(x, y) = (2, 3) 

Vf, D- (3, $) = (2,3) G, 225 
The direction of maximum rate of change is V f (1, 1) = (2, 3). A unit vector 
that points in this direction is found by dividing by 4/13, the magnitude of 


this vector. The result is (2713 3413, 


13. 13 
IVf(, D| = v13 
A line perpendicular to the one containing V f (1, 1) = (2, 3) would have 


slope =. A vector in this line is (—3, 2). A unit vector pointing in the same 


direction is (-34B, 2413). 


V f (x, y) = Qxy y), x? + 3xy’) 

Vf, 1)- 2, $) = 3,4) (6,225 

The direction of maximum rate of change is V f (1, 1) = (3, 4). A unit vector 
that points in the same direction is (à, 2). 

IVF, Dl =5 

A line perpendicular to the one containing V f (1, 1) = (3, 4) would have 
slope i. A vector in this line is (—4, 3). A unit vector pointing in the same 
direction is (-3, 2). 
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Problem 66 


1. First, we find the critical points by looking for places where V f(x, y) = (0, 0). 
Setting the vector V f (x, y) = (y + 2, x — 3) equal to (0, 0), we see that x must 
be 3 and y must be —2. So the only critical point is at (3, —2). Now note that 
D(x, y) = —1 for all (x, y). In particular, D(3, —2) = —1 < 0, so the critical 
point corresponds to a saddle. 

2. Setting V f (x, y) = (3x? — y, —x + 2y) equal to (0, 0), we see that the coor- 
dinates of any critical point must satisfy the system 


3x7 -y =0 

—x+2y=0 
The second equation tells us x = 2y. Plugging this into the first equation gives 
3(2y)? — y = 0. Solving for y then gives solutions at y = 0 and y = 5- Plug- 
ging these numbers back into x = 2y gives us x values of 0 and d So the 


coordinates of the critical points are (0, 0) and G, 5). 
We now compute 


D(x, y) = (6x)(2) — (—1)* = 12x — 1 


We conclude D(0, 0) = —1 < 0, so (0,0) is a saddle. On the other hand, 
DG, b) = 2 — 1] = 1 > 0, so there is either as maximum or a minimum at 


(1, 1). As ££2, 4) = 1 > 0, it must be a minimum. 





Problem 67 

All second partials are equal to — sin(x + y), so 

efe fN 

zi s — (24) = sin? (x + y) — sin? (x + y) 20 


Inspection ofthe graph, pictured below, shows that there are infinitely many maxima 
and minima, but no saddles. 


D(x, y) = 
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Problem 68 


Recall that D(x, y) = e e y^ Since (Ë 


dy? : y 
to zero, Es i — (ŻE » -)? can only be positive if ££ 32 i z is positive. But if the product 


of two numbers is positive, and you know the first of those two is also positive, then 
the second must also be positive. 





£)? is always biggerthan or equal 


Problem 69 


1. First we look for critical points in the interior of D by setting the partials equal 
to zero: 


af 
—=2x=0 
m X 
of 
— =2y+2=0 
dy 
From the first equation x = 0, and from the second y = —1. So (0, —1) is a 


critical point. (To find out what type we look at the matrix of second partials: 
| 2 0 


=4>0 
T 


Since un = 2 is also greater than 0, the point (0, —1) represents a minimum.) 
We now look for minima on the boundary of D using the method of LaGrange 
multipliers. The boundary of D satisfies the equation g(x, y) = x? + x =k 
Hence, we must look for points (xo, yo) that satisfy this and 


V f (Xo. Yo) = AV g (x0, yo) 
Note that V f = (2x, 2y +2) and Vg = (2x, 2 Hence, we get two equations: 
1 
2x =X(2x), 2y+2=d G») 


The first equation tells us either à = 1 or x = 0. If A = 1 then the second 
equation becomes 


y 
2 22-— 
yu 2 
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Solving for y then gives y = —#. However, if we try to find corresponding 


values for x by plugging this into x? + = | (the equation for the boundary 
of D) we get imaginary answers. 

If x = 0, then the equation x? + = = ] gives us y = +2. So on the boundary 
of D the only potential minima are at (0, 2) and (0, —2). 


2. To check which critical point represents the minimum of f we just plug them in: 


Tul -1)=1-2-1= -2 
f(0,2)=444-1=7 
(OD 24-4212) 


So the smallest value attained by f(x, y) on D is —2, and this happens at (0, —1). 


Problem 70 

1 3 16 23 
T [ec 3 [23]=0 
Problem 71 


Since V and W are parallel there is a number k such that V = kW. Soif V = (a, b) 
then W — (ka, kb). We now compute 





a b 
= akb — bka = 0 
ka kb 
Problem 72 
1 2 3 
1 1 0 2);=10-4)-2(-34+4)+3(2-0)=-4-2+4+6=0 
—2 2 -3 
0 3 


$ 1$ 1[|=0—2-0-10-2+30-4=0+1-12=-11 
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Problem 73 


i jk 
1L |1 2 3|210,—4,2) = /4 16 4 2 2/6 
-1 0 1 





2 |1 1 O|[2]ti; 1, -D| S / 1-141 — «3 


ij 
3. |1 2 
3 1 2 





o zx 


—|0,7,—5)| = V1 +49 + 25 = 54/3 


Problem 74 

1. First we compute the cross product: 
i j k 
1 2. O}= (2,-1,-1) 
1 1 1 





The magnitude of this vector is /4 + 1 + 1 = 4/6. So a unit vector that points 
in the right direction is 


(2,-1,-1) [v6 6 46 
J/6 |. 413^ 6" 6 





2. The cross product is 
i j k 
1 1 2/2(-L1,0) 
1 1 1 


The magnitude of this is /2. So a unit vector that points in the right direction is 


(=1,1,0) (2 42 o) 





J2 = 2^ 2^ 
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Problem 75 
Let V = (a, b, c) and W = (c, d, e). We now compute V x W and W x V: 


i j k 
VxW=ļa b c|=i(bf —ce) —j(af — cd) + k(ae — bd) 

d f 

i j k 
WxVe=ld e f|=i(ce —bf) —j(cd —af)+k(bd — ae) 

c 


Problem 76 
Let U = (a,b,c), V = (d,e, f) and W = (g, h, i). Now compute: 


























i j k 
U-(Vx W)=(a,b,c)-|d e f 
g h i 
d d 
= (a,b,c) - i fi es f| : 
h i g i g h 
_ le f d f d e 
CP ui JH i zi h 
a b c 
=|d e f 
g h i 


Problem 77 


Inspection of the following figure reveals that |W| sin 0 is the height of the paral- 
lelogram, and |V | is the length of the base. The area is just the product of the base 
and the height. 


W |W] sin 8 
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Chapter 7 Quiz 
Problem 78 
1.(a) The gradient of f (x, y) at (x, y) is 
V f (x, y) = (2x + 3y, 3x) 
The largest slope at (1, 1) is the magnitude of the gradient 
IVFd, DI = 20) + 305, 30))] 
= (5, 3) 
= /25+9 
= v34 


(b) To find the critical points we set each partial derivative equal to 0: 
2x+3y=0, 3x =0 


From the second equation x = 0. The first then implies y = 0. So the only 
critical point is at (0, 0). 
(c) We compute 


ef Pf 
dx? — axay f j| 
= = —9 
Pf Pf 3 0 
aya dy? 


(d) Since the determinant of the matrix of second partials is negative, the critical 
point (0, 0) is a saddle. 


2.(a) The magnitude of V is 


(1, 2, 3)|) = V12 +22 +32 = V14 


Hence, the desired unit vector is 


V (1,2,3) dc V14 | 











IV 14 14' 7" 14 
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(b) Since V - W = |V||W| cos 9, it follows that 





V.W 
os 0 = 
IVIIW| 
Aden a 
/144/3 
| 6 
| 4.42 
| N42 
EE. 
(c) First we find the cross product: 
i j k 
VxWe=]1 2 3 
1 1 1 


—i2-3)-jü —3) - k(1 — 2) 
al 2 audi 


The magnitude of this vector is 





|—1,2, 21) = V1 +4+1= V6 


Hence, a unit vector that points in the desired direction is 





(-1,2,-1)_/ V6 v6 v6 
J/6 | 63" 6 


Chapter 8: Calculus with Parameterizations 
Problem 79 


The speed is the magnitude of the derivative of the parameterization: 





|(2t cos t?, —2t sint?)| = V/ At? cos? t2 + At? sin? t? = 2t 


211 


212 Advanced Calculus Demystified 


Problem 80 


A tangent vector is given by the derivative of the parameterization: c'(f) = 
(3:2, 49, 51^). The point (1, 1, 1) = c(1). A tangent vector at this point is then 
c'(1) = (3, 4, 5). The magnitude of this vector is 4/9 + 16 + 25 = 54/2. A unit 
tangent vector is thus 








(3,4,5) [342 242 V2 
5/2 \10’ 5'72 


Problem 81 


1. The derivatives of the parameterization are 


ow : 
aos (cos 0, sin@, 2r) , and 


ow 
*» = (—r sin 6, r cos 0, 0) 


If r = 1 and 0 = 4 these vectors are (8, 1,2) and (—2 A, 0). 


2 2° 
i j k 

2.) 2 pee 11 
EF 

Problem 82 


1. First we compute 








dY ; 
— | = |(—sint, cost, 1)| 
dt 
= vsin? t + cos? t + 1 
—AX2 
1 _ v2 


Hence, if we multiply Y by y; “a We get the required parameterization: 


(t) = (ss V2 unt, 2) 
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2. First, note that 





“IS 


= sin f, COS f, 
dt 2 2 


d (T $) 


It now follows that 


d*® /2 ALAS. 
—— —(—-——ocost, ——-— sint, 0 
dt? 2 2 


We now compute the dot product 


do e [Xs 2 n Uo +) (- V2 E . 
COS f, cost, —-—— sint, 0 
2 2 
1 
cg 
=0 








2 
sin f cos f — ; costsint +0 


3. We compute the magnitude of m 


2 2 
= (2) + (Fam) +0 


1 
— cos? t + —sin?t 
2 * 2 





d? 
dt? 











Il 
rm 


SIS Nile 


Hence, if we let V = J2472 = (— cost, — sint, 0) then N is a unit vector and 


deo /2 


d? 2 


So, the required x is x 
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4. We compute 








B aT N 
=— x 
dt 
J2. V2 /2 . 
=(— sin f, COS f, x (— cost, — sin f, O) 
2 2 2 
i j k 
= —*Y? sint 2 cost s 
—cost —sint 0 
2 2 2 
= ( — sint COS f, 
2 2 2 


We now compute the magnitude of B: 


|B| = (<P) + (2) + (2) 


/3 "UM 
-an — COS = 
2 2 2 








=1 





5; — = (2 cost, 2 sint, 0) 
6. Since E = ee COS f, Y? sin t, 0) and N = (— cost, — sinf, 0) it must be that 


— v2 
T=: 


7. First, notice 





dN 
a = (sint, — cost, 0) 


We now do some calculation: 


TT. V2 ost 2 
"3 





inf, 
2 2 


+O (Pn V2 3 


ao gH 
—k — + tB = — 
dt 


qo erc. OS f, — 
2 


Answers to Problems 215 


1, 1 1 l. 1 1 
—(-—sinf,—-—cosft, ——)-F(—sinf, ——cosft, — 
2 2 2 2 2 2 


= (sint, — cost, 0) 
dN 
"dr 
Problem 83 
The length of Y(t) is computed as follows: 


1 
J dY 

dt 
0 





1 
ar = [ sint, cost, 1) dt 
0 


1 


= [ eos? re sin? +1 at 
0 





E. 


= vii 
=v 


1 
0 


Problem 84 
The length of W(t) is calculated as follows: 


1 1 


dY . 
——| dt = [ |(—sint,cost,t)| dt 
dt 

0 0 





1 
= | Votre sin? +7 dt 
0 





1 
= | TFP ar 
0 


This is precisely the integral we obtained for the length of (tcost,fsint) in 
Example 8-5. 
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Problem 85 
1 
[fen ds = J fe owe ma 
C 0 
1 
= f evel 
0 
2 
= [iia 
= 4 uau 
1 
1 TH 
= -u 
6 | 
-loi 1) 
|. 6 
Problem 86 
id 
Area= | | JOP Glas dy 
0 0 
11 
= | [visas dy 
0 0 
14 
Problem 87 
1 
af \? af \? 
Area = 2 — l dx d 
ef e 
0 
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2. Note that the surface over the domain D is a surface of revolution. In cylindrical 
coordinates it has the equation z = r. So the area is given by 


1 
| 2 
Area=2n f (£) +1dr 
dr 
0 


1 


=2r | VP+ ar 
0 
1 


=r | rar 


0 


= V2n 


Problem 88 


A parameterization for the surface is given by 


V(0, p) = Cf (9) sing cos 6, f (9) sing sind, f($) cos $) 
The partials of this are 


ed (— f ($) sind sin 0, f (6) sin $ cos 0, 0) 


o0 
ow — 
ap ^ (UO sing + FG) cos 4) cos, 

(f (9) sind + f (4) cos 9) sind, f'($) coso — f (4) sin) 


The cross product is 


» x = = (f ($) f ($) sin $ cos $ cos 0 — f (à)! sin? $ cos 0)i 
-FCf (9) f' (0) sing cos ¢ sin — f (9)? sin? $ sin)j 
—(f (9) f (9) sin’ $ + f(b)’ sing cos ¢)k 
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With some work, the magnitude of this vector is 











On ow 
EX eto sino F'O) + f (0 
And so, 
2x 
- | J FO) sin by (FO)? + f (09 dO do 
0 
E f FO sindy F'O)? + FO do 
Problem 89 


m 


Ov oY 
pe ei^ d0 dr 
0 


i 
[ 
| 
—— e 


k 
cos Ó no 0| dO dr 


—r sinô rcos0 1 


Oe 


1 

j 

l1 x 

= f | VTF a0 ar 
-1 0 


1 


=x | Viera 
RS, 





24-1 
— 2 2 In Mer (An integral table helps here.) 
2 \Vv2-1 


Problem 90 


The portion of the plane specified is parameterized by 


W(x, y) = (x,y, x + y) 
0<x<1l, 0<y<l1 
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We now compute 


1 


1 
[foro =f fs (W(x, |e x e| ax ay 


1 


1 
=f [iem 1,—1, 1)| dx dy 
0 


0 


T d 
V3 | | 2x+2y dy 
0 0 


1 


= v3 | x? + 2xy|, dy 
0 


1 
=V3 [1+2yay 
0 


— 43 (y 4 y); 
= 24/3 


Problem 91 


A parameterization is given by 


V (0, $) = (sin $ cos 0, sing sin 0, cos $) 
O<0<2n7, 0<¢<5 


Now calculate: 
5 2x 


| (e pode JJ fo, o» um 9" | dodo 


$2z 


= n cos | (— sin? $ cos 0, — sin? $ sin 0, — sin $ cos $) |d0d 
0 
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cT 
> 2m- 


= f [ cos@sing ao ag 
0 0 


= 2x f cos sing do 


Cis o 


x | sin2ġ dd 


ll 
H 


Problem 92 
A parameterization is given by 
V(r,0,z) = (rcos0,rsin0, z) 


O<r<R, O0Ox0-x2m, Ozz-h 


The derivatives of this are 


hU . 

— = (cos 6, sin 6, 0) 

or 

ow . 

— = (—r sin, r cos0,0) 
o0 

ov 

— = (0,0,1 

X ( ) 


We now compute volume by integrating: 


h 2x R| cos@ sind 0 
volume = | f f —rsinÜ rcos@ 0| dr d0 dz 
00 0 0 0 1 
h 2n R 
= | f f -dr doaz 
00 0 


h 2x 


R? 
=| |F ea 
0 0 
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= fr dz 


© 


= zx R?h 


Problem 93 


First, we calculate the derivatives of the parameterization: 


ow 
— ((1 — z)cos@, (1 — z) sin0, 0) 
ow . 
$3," (—(1 — zr sin6, (1 — z)r cos 8, 0) 
ow 
— = (—rcos0, —rsin0,1 
iE 
Now we integrate: 

^ff aw ov ov 

Volume = f | — — —| dr dé dz 

or 0 oz 
0 0 


(1 — z) cos 0 (1 — z) sin0 


a 


Il 
dinner 
ST 


= om bp oS onr CO, 


—rcosé —r sin 0 


(1 — z)°r dr dé dz 


Fg, a 
Mec 


3 


d0 dz 





(1— z)? 
2 


2 


a 
N 


Il 
E 
os 
= 
| 
a 


pes 





Il 

a 
2T 
N 

| 

N 

N 
+ 
w| = 
N 

W 
Nm 


© 


wl a 


—(l1—z)rsinO (1—z)rcos0 0| dr dé dz 
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Problem 94 


Letting R denote the elliptical region given by the parameterization W, and 
f(x, y) 2 2x — y, we have 


[25 v c f [rom o» dr d0 
R 


m 


. 2cos0 sin 0 
= J J [2(2r cos 0) -rsing]| dr dO 


= 


—2rsinO rcosO 


e 


0 


(4r cos 0 — r sin0)) (2r) dr d0 


Il 
oS 
oo 


8r? cos 6 — 2r? sinÓ dr dé 


Il 
ea 
ex L 














f 8 2 
m -r)cos0 — 2r?sin0| dé 
3 3 0 
0 
f 8 2. 
= [5556 - $ sino a 
3 3 
0 
8. 2 " 
= -sin + — cos 0 
3 3 0 
_ —4 
3 
Problem 95 
i 2 2m 2 à 
ow ow ov 
—— — dx d a= | | | dr d0 dw 
Ilr " ae & 3 3e 
V -20 1 
2 2m 2 cosh o cos 0 cosh o sin 0 0 
=|] |= J —rcoshosinO r cosh æ cos 0 dr d0 dw 
* 0 1 cosh oO | + sinhwcos@ rsinhwsin@  cosho 
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co 


Il 
o —À 


2 
cere 
— — — dr d0 do 
1 


27 


0 
2 27 
f | 559 
= — COS 
2 w 
-2 0 


2 2x 


=] f coho ad do 


-2 0 


Il 
o —À 


2 
[e dr d0 do 
1 


2 
d0 do 
1 





2 
= 3x f costo dw 
23 


= 37 sinh o^, 


= 67 sinh 2 


Problem 96 


The region V is easily parameterized using cylindrical coordinates: 

Psi(r, 0, z) = (r cos 0, r sin, z) 
To use this parameterization to evaluate the integral, we will need the determinant 
of the matrix of partial derivatives of Psi: 


cos 0 sind 0 
—rsinOÓ rcos0 O|=r 
0 0 1 


224 Advanced Calculus Demystified 


We now integrate 


l 
fw dr d0 dz 
0 


Ci 


[ema 
V 


zd0 dz 


NI =|= 


— Ci 


mz dz 


Il 
N o 
3 


Problem 97 
The region R is easily parameterized with polar coordinates: 


V (r, 9) = (r cos 0, r sin 0) 
0<r<2, 0<0<7 
The determinant of the matrix of partials of V is r. Hence, 


2 


(r sin) 
fino f ent) ar de 
0 
Lf 
| 


(1 + tan? 0) yr) dr d0 


ecd 


r sec? 0 dr d 


Hs 


2 sec? 0 d 


Il 
Os 


tan 0|o 
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Chapter 8 Quiz 

Problem 98 

1l.(a) Y(y) = (sin y, y) 


(b) The point (0, 0) is V(0). So a tangent vector is given by the derivative of the 


parameterization, T (cos y, 1), at y — O. This is the vector (1, 1). A unit 
tangent vector is thus 





Q5 G1) [4242 
[ap xm» \2°2 


2.(a) W(x, y,t) = (x, y, tQ + y?) 
O<x,y,t<1 


(b) We will need the partials of Y: 


ow 

<= (1,0, 2xt) 
ow 

— = (0, 1, 2yf) 

ow 

uem (0, 0, x? + y?) 


Next we must compute the determinant of the matrix of partials: 
1 0 2xt 
01 t |=x?4+y’ 
O 0 +y 


Finally, we integrate 


bot d 
; t 
[[fmewa-ff[S GE II e ty d dy d 
x? + y? 
Q 0 0 0 
1 
/ 


tx? + ty” dx dy di 


oO 
oN 


226 Advanced Calculus Demystified 


ilia 
3 


t + ty? dy dt 


w| = 


1 
3 


Il 
Wile ex- 


3. First, we compute the partials of $: 


= (cosht, sinh f) 


e|[& |8 


= (r sinht, r cosht) 


Next we will need the determinant of the matrix of partials: 


cosh f sinh f 








rsinht rcosht 


Now we integrate 


1 1 


Tf x? — y? dx dy = J fe cosht)* — (r sinh t?) (r) dr dt 
R 


-1 0 


i od 
= f| [Para 
-1 0 
1 


d 


dt 


BRI 


NI 
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Chapter 9: Vector Fields and Derivatives 























Problem 99 
SORORE uec ee dé UM dovE ae FR os a a wg 
SNNNNNANANA NV GERD TDA A AZZ ee 
SNS NNN NNNNA ANA TD DD AA AD DA a 
SNNANNNNNNAVA TDD DD Ass eee 
VN NNNM DA AZ AA Ae ee 
IAS S SUN NN NT PITS PA AAA ee 
"RE E IE ee ee ee 
zi us dec ue ae "c CE WAT D le iste a So EM E 
| NC DP | Oe Ba SUE DORE 
-——————-c-4g4lN S M———————- 
coe emet ett at at u^ i I e NS Se Se Se DLE 
aeta ma ant YE fe V NN CNN SOS Se SS 
EE ee ZO ah ^ Po ONG NC NO SENS Sce USO 
Ce eee LAL f gE li NN MT 
Cu OO MEE AL ^ odo quU NNN ONG Se Ne ONE 
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Problem 100 
1. (1, D) 2. (1,z, y) 3. (y+tz,x+z,x+y) 


Problem 101 





a a 
(y) + ay Ot y OOo 


gU »tLG-»*gO-1-141-1 


N 
Fle Sle 


a 9 2 2 82 2 0 ay 
— (x? + mae Or +—(z^) =2x — 2y +2z 
| y^) E y^) x) y 


Problem 102 





vvs=v (2,2 2) 
ac Oy’ a 

-8 VELLET 
tarta 


Problem 103 
We let F = (fx, fy, f.) and G = (8x, gy, 8z). Then 


FxG= (fy& = fy TBs fun fg = Ty8x) 


and so 


0 0 
V-(FxG= ay PE m fz8y) + ay i — fxg) 


ð 
T (gy m Fy8x) 

















| 0f, og, Of, 0g, 
te dl ae mS i. 
of. Og x = ee 
E&x + fz -p 
dy dy ay © ay 
à fs 0g, 0f, 0gx 
ar E 


ac ac 
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of. 8f, 9f, f. Of, s 
ð z z Q à ay 


= (8x, Sy: gz): ( 





0g. 0g, Og, Og. Og, Og. 
=( fo fy, fa | nd 2 [S 2 


dy d z a ð 
—G-.(VxF)-F.(VxG) 
Problem 104 
1l. curl V =V x V 
i j k 
à à à 
dx dy az 
X+yY x—z ye 
= (2,0, 0) 


2. (2x, 0, —2z) 











Problem 105 


of of z 
1. V V =Vx(—, >, > 
x (Vf) e ay’ k 

ij k 
—_{/% 2 2 
— | ax ay a 
af af f 
ax ay a 








-(2- 2 Pf gf DE 
^ A8yüz — züy Ozüx ade’ dxOy — Oyàx 
= (0, 0, 0) 





3 vos-viz-i af _ ðh ag r 


dy & & dà ay 

2: (8-8)-2 (89,2 (88r) 
ax Nagy a dy \ae & z \ax ð 
_ Êh Pg h Pf fg Pf 

^ axdy axa — Oydx  ðy&æz — dküx  &zðy 








=0 
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Chapter 9 Quiz 
Problem 106 





i j k 
1. a a 8 
üx ay az = (0,0, 0) 
x2 + y? x? + y? 


0 0 
2. xc ty yt a0 30-22) 


3. Qxsin(y — z), x? cos(y — z), —x? cos(y — z)) 


= ig Of af ] 
4. Suppose W = (f, g, h). Then the last component of V x W is = — ay" This 
Y if et ig y of _ y - 
can equal - if either = = - or ea ees In the former case g(x, y, z) = yInx. 
In the latter case f(x, y, x) = — L, So two vector fields that have the required 


curl are (0, y In x, 0) and (—ż, 0, 0). 


Chapter 10: Integrating Vector Fields 
Problem 107 


The circle is parameterized by 
W(t) = (cost, sint), O<t<2z 


Note that this parameterization agrees with the specified orientation. Now we may 
integrate 


27 dY 
[weas= [ ww ar 
C 0 


27 
= J (— sint, cost) - (—sinf, cost) dt 
0 
2n 
= [sor 4- cos? t dt 
0 
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Problem 108 


dY 
[weas= [ww at 
C 0 


1 
= fP.Pa oh r+ tn) Qn 1-1) at 
0 


af OF 4f = 1 at 


I 
esL 


1 
= E i t 
s 2 3 r 
3 1 1 
LM NE 
5 2*3 
EHE: 
|. 30 


Problem 109 


1. First, note that V f = (y?, 2xy). Now we integrate 


; dV 
[ vr-as= [vr Sat 
" dt 
1 


2 
= fez - (1, 2t) dt 


2 
= | 5t^ dt 
-1 
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m gr. 
-32- (-D 
= 33 


2. FOA) -— fC) = FG, — f(—1, 1) 
23) 
= 33 


Problem 110 


The curve C is parameterized by 


W(t) = (t,t, O<t<1 


The derivative of this parameterization is (1, 21). 


1 
[v9 [en maam dt 
GC 0 
1 


= fa 


j= °> 


Problem 111 
The unit sphere is parameterized with spherical coordinates by 
V (0, $) = (sin $ cos 0, sing sin 0, cos ġ) 
0<0<2m, O0x$ó-zm 
The derivatives of this are 


ow 


a (— sin $ sin 6, sin $ cos 0, 0) 


* = (cos $ cos 0, cos d sin 0, — sin $) 
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And so, 


i j k 
Ov | ov : ; ; 
$»' 3 —singsind  sin$ cos 0 
cos@cos@ cos@sind —sing 


= (— sin? $ cos 0, — sin? $ sin 8, — sing cos ¢) 


Notice that at the point Y (0, =) = (1, 0, 0) the vector A x 2 is (1,0, 0), which 


agrees with the specified orientation. We may thus integrate to find the correct 
answer: 


m x 


[vas f [ wem n (27) do do 
S 
0 


ap 
0 


na 2x 
f [ sin 6, sin ó sind Te a) 646 
= sin cos g, sin sng, COS . — x — 
wo * 
0 0 


-T 2x 


=| | -sing do do 
0 0 


T 


= J —2m sing dd 


0 
= 2x cos lo 


= —4r 


Problem 112 


First, we calculate the partials of the parameterization: 


= (1,0, 2u) 


2/2 ge 


= (0, 1.2) 
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Now, the cross product is given by 


k 
2u | = (—2u, —2v, 1) 


2|z 
e| g 
O = =e 
— o œ 
S 


Finally, we integrate 


10d 
[was f [ won (FF) du dv 
ou Qv 
s 0 0 


1 
[ «e + v2), v(u? + v2), 0) - (-2u, —2v, 1) du dv 
0 


I 
oN 


11 
J J -2u (u + v?) — 2w (u? + v?) du dv 
0 0 


10 
J | m6 - 41202 - 204 du dv 
0 0 


Il 
oe, 
| 
VIN 
| 
| 
ce 
to 
| 
N 
c 
RR 
a 
iso) 


Problem 113 


First, we calculate the partials of the parameterization: 


= (— cos $ sin 0, cos ¢ cos 0, 0) 


= (— sin $ cos 0, — sing sin 8, cos d) 


Sz mE 
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Now, the cross product of these vectors is 


i j k 
oy ov f 
we —cos@sind cosdcosé 0 
—sin$cos0 —singsind cosó 


= (cos? $ cos 0, cos? $ sin 0, cos $ sin $) 


We are now ready to integrate 


E 2x 
1 1 
[weas= || — 0 
cosġcosð cos¢@ sind 
S 0 





AIA 


-(cos? $ cos 6, cos? $ sin 0, cos $ sin 9) dé dọ 


T 20 
J [5599940 
0 


dE 
4 


4z cos $ do 


Il 
iy 


AIA 


= 4r sin$l|*. 
4 
=4nJ/2 


Problem 114 


First, we compute V x W: 


A parameterization for S is given (via spherical coordinates) by 


V (0, $) = (sin $ cos 0, sing sin 0, cos $) 
T T 
0<8<—, 0<ġ<— 
SS STE 
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The partials of this are 


ow 
E (— sin $ sin 6, sin ó cos 0, 0) 
ow 


a (cos $ cos 0, cos à sind, — sin $) 


And so, 
OP ow j í 
o*a —singsind  sin$ cos 0 
cos$cos0  cosósinO —sing 


= (— sin? $ cos 0, — sin? $ sin 8, — sing cos d) 


Now notice that at the point ws, 3) — (0,1,0) the vector w x a is equal 
to (0, 1, 0), which disagrees with the specified orientation. We must therefore 
remember to switch the sign of our final answer. Finally, we can integrate 


[esas f [wee oy (x 2) dé dọ 
d E 0 M 


OE] 
ESE] 


(sin $ sin 0, sin $ cos 0, 0) ne dé dọ 
sin ¢ sin 8, sing cos 8, 0) - | — x — 
© æ 


II 
— 
— 


0 0 


a 


Z 4 
= f J -2sin*ġsing coso dé do 


0 


2 4 
= | [ -sin’ osin20 do ag 


0 0 


c 


x 
4 


do 


0 


sin? $ cos 20 





Il 
OS S 
NI = 
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-$ü — cos? $) sing dọ 


Oe uis 


l 1 ial 
5 cose S99 i 


The correct answer is therefore i 


Chapter 10 Quiz 
Problem 115 


1. A parameterization for C is given by 
v()-(^,0, 0stx1 
The derivative of this is 


dY 


E = (2t, 1) 


Note that at the point V(0) = (0, 0) this vector is (0, 1), which disagrees with 
the orientation indicated in the figure. Thus, we will have to negate our final 
answer. We now integrate 


1 


fonan- fan ana 


C 0 


-fusia 
0 


=-—2 


rt 
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2. The partials of Y are given by 
= (cos 0, sin 0, 1) 


= (—r sin0, r cos 0, 0) 


eZee 


The cross product of these vectors is 
i j k 
— x — =] cos0 sind 1| = (-—rcosé, —r sinO, r) 
—rsinü rcos0 O0 


We are now prepared to integrate 





2 
[t= ,r)-(—rcosé, —rsin0,r) dr dd 
rcos@’ rsin@’ 
1 


A 
sims 
=x |= 
| 
=< |= 
N 
ND. 
1 


—1+1+r°dr do 


Tnv 


N 
3 


—9 20 





Il 
w| = 


wi N 
&, 
D 


ÓO e 
A 
3 


EJ 
Chapter 11: Integration Theorems 


Problem 116 


Notice that W = (1,1,1) = V(x + y+ z). Hence, if we let f(x, y,z)= 
x + y +z then 


[ v fa. 1 0- ra. 023-221 
C 
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Problem 117 
The trick is to notice that W = V f, where f(x, y) = x sin y. Then 


[weas=r(2 =) - £@,0) 
c 


T 
= 2sin — — Osin 0 
Sins sin 


=2 


Problem 118 


First, notice that f (x, y) = x? - Ax +4 + y? c 2y 1 2 (x 22 4 (y+ 1". 
Thus the graph of f(x, y) can be obtained from the graph of z — x? + y? 
(a paraboloid) by shifting 2 units in the positive x-direction and 1 unit in the 
negative y-direction. This puts the “bottom” of the paraboloid at the point (2, —1). 
Hence, if (a, b) is any point other than (2, —1) then f(a, b) > 0. The result now 
follows easily: 


[ vids fe - ro - no fab) -0 
C 


Problem 119 
Suppose the beginning and ending point of C is p. Then 


Jt- ass= to - foo 


C 


Problem 120 


Suppose f(x) is a function. Let V be a parameterization of [a, b] given by 
V(x)—-x, ax<x<b 


By the theorem, 


[vr-as= fb) — fla) 


C 
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We now follow the definitions: 


b 


/ dV 
[viis [ren Sas 


Putting this all together gives 


b 
[ro dx = f(b) — f (a) 


Problem 121 


When you are coming down the mountain the force of gravity is working with you, 
not against you. So the work done against gravity during this portion of the trip is 
negative. This cancels out the extra positive work you did to get to the top of the 
mountain. 


Problem 122 


ð 
LX 
ax 


J (-y?, x!) -ds 


3Q 


] 


onL 


2 5C» d 


2x 4- 2y dx dy 


Il 
oO 
oN 


1+ 2y dy 


Il 
N en 
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Problem 123 


Let D be the horizontal line segment that connects (0, 0) to (a, 0). Then CU D 
bounds a rectangular region Q. Green’s Theorem says 


J oa [ions J (y, x) -ds 


C D CUD 


SEND 


=0 


We conclude 


J 0s = [mas 


C D 
and so the answer only depends on a. 


Problem 124 


By Green's Theorem, 


MED 
[weas= [uo p a 
ðo 


o 


= [2m asy 


o 


To evaluate this integral we will need the determinant of the matrix of partial 
derivatives of $: 
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We now integrate 


[2+ 2y de dy = (2Qu — v) + 2(u + v)) 3 du dv 


o 


Se i 
Tr 


18u du dv 


Il 
e=- 
Tt 


27 dv 


Il 
—- 


Il 
l2 o 
~ 


Problem 125 


The area of any region Q in R? can be calculated by the integral 


[fee 


Green’s Theorem says that this will be equal to the integral of W = (f(x), g(x)) 
around 0Q if 


8 əf a 


& ð 


One suitable choice for W would thus be (0, x). We now integrate this vector field 
over the unit circle, by using the usual parameterization 


W(t) = (cost, sint), O<t<2z 


2x 
[w-as= f sn - (— sint, cost) dt 
aQ 0 


2x 
2 
= f cos t dt 
0 
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27 


E in 2t 
= -t + -sin 
2 4 0 


Problem 126 


1. Let C, and C; be two circles centered on the origin, oriented counterclockwise. 
Let Q be the region between these circles. Then Q has two pieces. One is C1, 
and the other is C5 with the opposite orientation. As a shorthand only we write 
this as C4 — C5. Green's Theorem now says 


o= f foaay 
Q 


z J (F(x), g(a) - ds 


C1—C3 


= f oe aen sds- f ro ae ds 


Ci C5 


And hence, f (f(x), g()) ds = f (f(x), g(x)) - ds. 
Ci C5 


; z i ag _ af _ 
2. Using the previous part, we just have to show = a =O 


à x ð  —y 
axx? +y? dyx2+ y? 


D 1 2x? —1 " 2y? 

x? + y? (x? + yy x2 + y? (x2 + y?)? 
2 2x? + 2y? 

x2 + y? (x? + yy 











=0 


rt 
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3. Note that on the unit circle x? + y? = 1, so, on this circle the integral of W is 
the same as the integral of (—y, x). 


2x 


[weas= [nds f (sin, cost) (sint eosn dt 
C 
C 0 


2x 
= [ 5 + cos? t dt 
0 


2x 
= fa 
0 


= 2r 


Problem 127 
1. By Green’s Theorem, 


[vs f Ze- PLE dx dy =0 
ðo 


2. The curve C is parameterized by 
Y) = (t,0) —-l<t<1 


The derivative of this parameterization is (1, 0), so 


[weas= [euin dt 
C 


Il 
—À 
- 

N 
a 
~ 
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3. Let D be the top half of the unit circle (oriented counterclockwise). Then 
C U D = 0c. So we have 


o- [was 
do 


= | weas 


CUD 
= [weast |w ds 
C D 
2 
-5+/w ds 
D 
Hence, 
2 
fw EET 
D 


Problem 128 


The surface S is a disk in the yz-plane. Its boundary is thus the unit circle in the 
yz-plane. The tricky part is determining the proper orientation on the boundary 
circle. A normal vector giving the orientation of S is determined by 
Qv ov 
— x — = (r,0,0 
TR (r, 0, 0) 
For r #0 this is a vector which is parallel to the positive x-axis. The orien- 
tation on dS is then given by the right-hand rule, as in the following figure. 
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The usual parameterization for the unit circle in the yz-plane would be 
$(t) = (0, cost, sint), O<t<2z 


Notice that % = (0, — sint, cost). At t=0, $(0) = (0,1,0), and “(0) = 


(0, 0, 1). This agrees with the orientation given in the figure. We are now prepared 
to integrate, using Stokes’ Theorem: 


[ovxwy-as= | Was 
S as 


27 
= / (0, 0, cost) - (0, — sint, cost) dt 
0 


2x 
= foi dt 
0 


hd i "u 
per — sin 
2 34 





= 
As the chosen parameterization $ agrees with the orientation on 0S there is no need 


to change sign. 


Problem 129 


Let Q be the region in IR? such that (x, y) is in Q if (x, y, 0) is in S. Then we may 
parameterize S as follows: 


Wix,y)=(, y,0), Q.»eQ 


It follows that Æ = (1, 0, 0) and Æ = (0, 1, 0). Thus 
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By definition 
à, 0 
[vw «s = [ (0.0.2 - 2). 0,0, 1) dx dy 
dc — Oy 
S Q 
à, 9 
= J a af dx dy 
dx Oy 
Q 
Also by definition 


[weas= | rouen sas 
às àQ 


But Stokes' Theorem says 


[x vas- [was 
S as 


So we get 
ae 8 
- Fax dy = [ (fo osa nas 
dQ 


SIR 


Q 
which is Green’s Theorem. 


Problem 130 


The boundary of S is parallel unit circles Co and C; in the planes z = 0 and z = 1, 
respectively. However, the orientations on these circles are opposite, as can be seen 
in the following figure. 

| z 
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A parameterization for Co is given by 
Wo(t) = (cost, sint,0), O<t<2z 
At the point Yo(0) = (1, 0, 0) the derivative — = (— sint, cost, 0) is the vector 
(0, 1, 0). This agrees with the arrow in the picture. 
Similarly, a parameterization for C, is given by 
W(t) = (cost,sint,1), O<t<2z 
At the point Yo(0) = (1, 0, 1) the derivative we = (— sint, cost, 0) is again the 
vector (0, 1,0). This is opposite to the arrow in the picture, so we will have to 


remember to negate the integral over this curve. 
Stokes’ Theorem says 


[ov xw)-as= W-ds 
d CoUC, 


E W-ds+ | W.ds 
Co Ci 


We do each of these integrals separately 


2x 
w-ds= | (0, 0, 0) - (—sint, cost, 0) dt = 0 
Co 0 


2x 
w-as= | (—sinf, cost, 0) - (— sint, cost, 0) dt 
C 0 
2x 
=) dt 
0 
= 20 


But, for orientation reasons, we have to negate our answer, yielding —27r. Our 
final answer is thus 


[ x W)-d8 = -27 +0= -2x 
S 
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Problem 131 


We assume the sphere is oriented with an outward normal (just reverse the picture 
if not). Dividing the sphere S into a northern hemisphere S4, and a southern hemi- 
sphere S. , we see in the following figure that the equator gets opposite orientation 
as the boundary of each. 


As the induced orientations on the equator are opposite, it follows that 


[sae fws 


0S. oS. 


We now employ Stokes’ Theorem 


[ew dS — LE ds+ [ (xw dS 


a 


aS. as. 
=0 


Problem 132 


The direction of V x W will be perpendicular to the plane containing the loop 
with greatest circulation. Since this is the x z-plane, the direction of V x W will be 
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parallel to the y-axis. The magnitude of V x W can be estimated by 
V x Will e c Wa 
x Z ————— - ds 
P^! 77 "Area(D) 
C 


where D is the disk whose boundary is the loop C. Since the radius of C is .1, the 
area of D is z (.1)?. This gives us 


1 1 
— | W.ds = ——~(5) © 15.92 
cm! S x iyi ) 
e 


We conclude that at the origin 
V x W z: (0, 15.92, 0) 


Problem 133 


The circulation of the vector field around any loop in a plane parallel to the xz-plane 
would be zero, since the vector field is constant on such a plane. A vector V that is 
perpendicular to such a loop points in the y-direction. The circulation around any 
horizontal loop would also have to be zero, since the vector field is perpendicular to 
such a loop. A vector W that is perpendicular to such a loop points in the z-direction. 
The curl must then be perpendicular to both V and W, and so must point in the 
x-direction. 

One can also see this algebraically. Such a vector field must look like (0, 0, f(y)). 
The curl of this is 


i j k 

0 ə a 

— — = |=), 6,0 
x » & (f Œœ), 0, 0) 
0 0 fO) 


Problem 134 


Let V be the unit cube. Then Gauss' Theorem tells us 


[was f f [macie 
oV V 
1 1 1 
=f f [9 wacara: 
0 0 0 
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2xyz + 2xyz + 2xyz dx dy dz 


Il 
inia m 
ox 
ex L 


6xyz dx dy dz 


Il 
de 
ox 
ex 


3yz dy dz 


Il 
diim o 
ox 


NIW 
N 
a 
N 


Il 


Problem 135 


By Gauss' Theorem the integral of V - W over the ball B is equal to the integral of 
W over the unit sphere S, with outward-pointing normal vector. The unit sphere is 
parameterized in the usual way by 


V (0, $) = (sin $ cos 0, sing sin 0, cos $) 
0x0-x2n, O0<¢d<az 
As in Example 11-8, 
Qv ov 
—, x p 
Oo 
which agrees with the orientation on S. 
We now compute 


J| [v mata [ w-as 
B S 


m 2x 


= J fo 0, e*95?) . (— sin? $ cos 6, — sin? ¢ sin 0, 
0 0 


= (— sin? $ cos 6, — sin? ¢ sin, — sin $ cos $) 


— sin $ cos $) d0 d$ 
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m 2r 


= | [ -sing cosge™? dé d$ 
0 0 


T 


= 27 — sin $ cos pee? dp 
0 


-1 


= 2n f ue" du 


1 
= 27 (ue" — g^ 
4r 


e 
Problem 136 
The region V is parameterized by 
V(r,0,z) = (r cos 0, r sin, z) 


l<r <2, 03037, 0xzx2 


We use Gauss' Theorem to transform the integral: 


[vase f f [v-wacaso 
av V 
= f | [5 y» ax dy az 
V 


We now use the parameterization to change variables: 


f | [58 ea 
V 


2 3 cosÜ sind 0 


2 
=f f [36e simo) —r sin r cos0 0| dr dé dz 
1 


E 0 0 1 
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3r?(r) dr dO dz 


Ti 


3r? dr dO dz 


Tt 


^| dé dz 





RI 


d0 dz 


E^ 


x co ees e MEN oi iy comes. id 


alg ecan Ee earn Sar 


Problem 137 


First, we observe that the polar equation r = cos @ is a circle of radius i centered 
on the point G, 0). Hence, the surface C is a cylinder of height 2. The surface D 
is a disk, which caps off the top of the cylinder C, like an upside-down can. Let E 
denote the “bottom” of the can. That is, E is the set of points in the xy-plane which 
are within i of a unit away from the point G, 0, 0). We assume an orientation is 
given on E so that C, D, and E together form the (oriented) boundary of V, the 


points inside the “can.” Applying Gauss’ Theorem gives us 


[ [ [v-wardyaz= f was 
V 


C+D+E 


= J W-dS+ | W-as 
E 


C+D 


Now notice that on E the vector field W = (0, 0, 0), so 


[w-as=o 
E 
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Putting these results together gives us 


J| [vw f was 
V 


C+D 


We may thus obtain the desired answer by evaluating the integral on the left-hand 
side of this last equation. The first thing we will need to evaluate this is the 


divergence of W: 


ð 
VW = —xyz = xz 
dy 


Next, we will need to parameterize V. Notice that V is a solid cylinder of radius P 
whose central axis has been translated i of a unit away from the z-axis, in the 
positive x-direction. A parameterization is thus given by 


1 
W(r,0,z) = (: cos 0 + 5 r sin, :) 


0<rz< 0<6<2n, 0<z<1 


, 


Nl — 


The partials of this parameterization are 


ov : 

— = (cos 0, sin 0, 0) 

or 

ow . 

— = (~-r sinO, r cos 0, 0) 
00 

ov 

— = (0,0,1 

x ( ) 


The determinant of the matrix which consists of these vectors is 


cos 0 sind 0 
—r sin rcos0 O|=r 
0 0 1 
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We now use V to evaluate the integral: 


J| |V wen f | [e 
V V 
2 


| (- cos 6 + ;) (z)(r) dr d0 dz 


0 


© 
v» I 


1 
2 


Il 
jac 
oo| 4 

N 

a 

N 


Problem 138 


Sı, and S5 with the opposite orientation, together bound a volume V of R3, 


C) 


Oriented surfaces with S, and S5, with opposite 
the same oriented boundary orientation, bound a volume V 


rt 
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By Gauss’ Theorem 


J| [v mata [ w-as 
V av 


S-52 
= [ w-as- f was 
$i S 


But V - W = O implies f f f V -W dx dy dz = 0. We conclude 
V 


[v [was 
$i Sp 


Problem 139 


1 
V. f : 
Map) ws" ae 
ðB 


1 
—4- 03 


in (1? 


~ 119.36662 


Chapter 11 Quiz 
Problem 140 
1. The key to this problem is to notice that W = V f (x, y, z), where f(x, y, z) = 


xy. So, 
[ [en 
C 


C 
-f,l,1)-7(0,0,0) 
edic 
=1 
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2. Using Green’s Theorem 


k (1, In x) EE a C2) = (0) dr dy 
= f| | -<avay 
x 


o 


To evaluate this integral we will need the partials of the parameterization: 


ed = (v’, 3u?v) 
Op 3 
T (2uv, u`) 


The determinant of the matrix of partials is thus 


v 3u2v 


3 |= = wv? — 6u3v —5u^v? 








2uv u 


We now integrate: 


1 
[[-e- 
x 


= 


2 
| AC Su^v?) du dv 
1 
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3. The curves C, and C; lie in the xy-plane, and are pictured below. 


Ci 


C» 





Let C, denote the curve C, with the opposite orientation. Then the curves C, 
and C; together form the (oriented) boundary of a region R of the xy-plane. 
Stokes' Theorem then says 


| W-ds= fV xW)-as 


C, UC5 R 


But the statement of the problem specifies V x W = (0, 0, 0), so the integral 
on the right (and hence the integral on the left) is 0. Now notice 


J W-ds= | W-as+ [ Weds 
Ci C2 


C; UC5 
= [ w-as- [was 
C; Ci 


Since this is zero, the desired result follows. 
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Final Exam 
Problem 141 


1.(a) 











(b) To find the critical points we set the partial derivatives equal to zero: 


af 

à 7t 

0 

of, 9-0 
ay 


The first equation tells us y = —1 and the second tells us x = 2. So (2, —1) 
is the only critical point. Now we compute 


ef af 
dx? axðy " j| 
= 2-1 
er sS 1 0 
oy ay? 


Since this is negative the critical point (2, —1) corresponds to a saddle. 


(c) First, we find a unit vector that points in the desired direction: 
(1,2) _ (1,2) _ [v5 25 
K,2] v5 5° 5 
The desired slope is the directional derivative in this direction: 
p ($ 24/5 


TM 





= S. 22 d 
5 5 


—24/5 
5 
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2 1 


(d) Volume = | f ay-+x—2y-+4dx dy 
0 0 
1 


dy 
0 





2 
J 55e pened 
= =X —X = £2VX X 
gp 7 > j 
0 





2 
1 1 
= | -y+--2y+4d 
[55 y+4dy 
0 
Le. 2 ? 
= y? + y- y? +4 
gt 4 
=1+1-—-4+8 
=6 


2 2 


2. Since r? = x? + y?, the desired surface has the cylindrical equation z = 4 — r°. 
Where the graph hits the xy-plane we know z = 0, and hence 0 = 4 — r?, or 
r = 2. A parameterization is thus given utilizing cylindrical coordinates by 


Wr, 0) = (rcos6,r sind, 4 — r°) 
O<r<2, 0€D-«2x 


3. There are several ways to do this. We will utilize spherical coordinates, since 
that is how the problem was stated. Being on a sphere of radius 1 says o = 1. We 
now plug this, and the information 0 = ¢, into the usual spherical coordinates: 


V (0) = (sin0 cos 9, sin @ sin 9, cos 0) 


To get the whole circle the domain should be 0 < 6 < 27. 


0 ð 
4 V-W= — xz? + x 


ax 
55 
=Z 4+ 2472 
j k 
vxw=|2 2 2 
xz 0 x2? 


= (0, 2xz — z, 0) 
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5. Since the region V is cylindrical, this is best done by utilizing a parameterization. 
We parameterize V in the usual way with cylindrical coordinates: 


Wir, 0, z) = (rcosO,r sin0, z) 
T 
O<r<l, DS S es 0<z<1 


To do the integral we will need the determinant of the matrix of partial derivatives, 
which simplifies to r. We can thus integrate as follows: 


151 


| 5m rh asay as [ f f 2 (ie reos + Csin?) (r) dr dO dz 


V 000 





i> 14 
= f | [rie ar do az 


0 0 0 


: T 
=f | | viaa: 
0 01 


NIA 


Q3 — 1) d0 dz 


WIN 


Il 
ox 
Oe is 


T edid 
3 A 


Il 
oN 


JT 23 
aou 

3 ( ) 
6. The derivative of the parameterization is 


d 
d$ = (—2sint,2 cost, 21) 
dt 


rt 


262 Advanced Calculus Demystified 


We now integrate: 


fo 0, x? + y?) -ds = | (0,0, (2cost)* + Q sint)?) - (—2sint,2 cost, 2t) dt 
C 


(0,0, 4) - (—2 sint, 2cost, 2t) dt 


7. First, we compute the partials of the parameterization: 


H = (cos 0, sin 0, 2r) 
ar 
ap . 
— =(- 0, 6,0 
3 (—r sin0, r cos ) 
The cross product of these vectors is 
i j k 
% x ai =| cosd sinô 2r| = (—2r? cos 6, —2r? sin 0, r) 


—r sin rcos0 0 


We now integrate: 
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8. The key to this problem is to notice that W = V f, where f(x, y, z) = xyz. 
Then, using the independence of path of line integrals 


[ [ens 
C 


C 


f (9(2)) - few» 


) — f (0,0, 0) 


9. The surface S bounds a volume V which is parameterized by 
Wr, 0, z) = (r cos0,r sin0, z) 
O<r<l, 0<0<2mr, O0<z<l 


To evaluate an integral over this region we will need the determinant of the 
matrix of partials. We have done this calculation several times. The reader may 
check that the answer is r. We now integrate 


[vs f f [VWa dyaz 
S V 
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INDEX 


A 
algebra, using with partial 
derivatives, 25 
arc length, relationship to line 
integral, 113 
area. See also surface area 
computing, 45 
computing for parallelogram, 105 
computing for region of xy-plane, 
45,119 
determining relative to change of 
variables, 119 
finding for parallelograms, 88, 91 
parameterizing, 65-67 
axes 
labeling as r and z, 52 
labeling when plotting points 
with multiple coordinates, 3 


boundary, determining for 
surface S, 161 


C 

chain rule 
computing derivatives with, 28 
using, 28, 30 





circle of radius one, parameterizing, 
12-13 
components, relationship to vectors, 69 
composition 
of functions of multiple variables, 
29-30 
relationship to parameterized 
curves, 26-29 
visualizing, 27 
constants, relationship to 
derivatives, 24 
continuous function, definition of, 21 
coordinate planes, plotting intersection 
of graphs with, 6—9 
coordinate system 
definition of, 49 
identifying type of, 15 
importance of, 3 
coordinates, plotting points with, 2—3 
Cross product 
determining, 91 
versus dot product, 91 
magnitude of, 92 
rewriting definition of, 92 
of vectors for surface of 
revolution, 11 
curl, geometric interpretation of, 
164—166 
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curl of vector fields 
computing, 150 
formula for, 149 
curl operation, defining relative to 
vector fields, 130 
curves 
finding lengths of, 100-101 
parameterizing, 101, 103 
cylindrical coordinates 
converting to rectangular 
coordinates from, 50 
expressing restriction on domains 
in, 62 
parameterizing domain of 
integration with, 121 
parameterizing surfaces with, 
109-110, 114 
translating to rectangular 
coordinates, 66 
cylindrical equations, graphing, 51—53 


D 


derivatives. See also partial derivatives; 
second derivatives 
computing with chain rule, 28 
determinant 
computing, 90 
for volume calculation, 115 
of a3 x 3 matrix, 90, 115, 117 
relationship to parallelograms, 
88-89 
writing for matrix, 90 
direction, finding for level curve, 78 
directional derivative, explanation 
of, 76 
divergence 
calculating for vector field V, 129 
defining for V, “Div V", 128 
definition of, 167 
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geometric interpretation of, 
171-173 

positive divergence, 172 

taking for vector fields, 166-171 

Zero divergence, 173 

domain of integration 

definition of, 33 

parameterizing with cylindrical 
coordinates, 121 

restricting, 62 

dot product 

versus cross product, 91 

definition of, 72 

finding cosine of angle between 
vectors with, 74 

purpose of, 88 

relationship to geometric 
interpretation of curl, 165 

source of, 77 

taking, 74 

usefulness of, 73 

writing directional derivative 
as, 76 


E 
earth's surface, determining position 
on, 97 
equations 
for cross product, 91 
for curl of vector field, 149 
for directional derivative, 76 
for divergence of vector field, 129 
for parabolas, 52 
of spheres in spherical 
coordinates, 62 
for surface area of graph, 108 
examples 
continuity of functions, 21 
cosine of angle between 
vectors, 74 


INDEX 


cross product, 92 

curl of vector field V, 130 

direction for level curve, 78 

divergence of vector field V, 129 

dot product, 73 

finding area of R, 45 

finding gradient for saddle, 80-81 

functions, 1-2 

Gauss' Divergence Theorem, 
167-171 

geometric interpretation 
of curl, 166 

geometric interpretation of 
divergence, 172-173 

gradient of function, 76 

gradient of f (x, y) at point 
(x, y), 127 

gradient vector, 77—78 

graphing cylindrical equations, 
51-52 

graphing functions, 5-6 

Green's Theorem, 153-155, 
158-159 

integral evaluation, 121 

locating points in rectangular 
coordinates, 50 

optimization, 83-84 

oriented curve relative to line 
integral, 135 

overcoming gravity, 147—148 

parallelogram area 
computation, 92 

parameterization, 64 
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signed volume of 
parallelepiped, 90 

Stokes' Theorem, 162-163 

surface area of sphere, 106-112 

surface integrals, 141—142 

translation from polar to 
rectangular coordinates, 66 

unit vector, 71—72 

vector addition, 70 

vector fields, 126-127 

vector subtraction, 71 

volume bounded by sphere, 116 

volume computations for 
integration, 40—42 


Figures 


contours of f (x, y) and domain 
D, 86 

level curves, 7, 9 

level curves forming saddle, 11 

local and absolute maxima, 85 

paraboloid z = x? + y?, 10 

perpendicular axes drawn in 
perspective, 3 

plotting points, 4 

point and vector, 70 

rectangle, Q, 152 

topographic map related to level 
curves, 8 

weather map showing level 
curves, 8 


formulas. See equations; examples; 
problems 
Fubini’s Theorem, relationship to 
integration, 36 
function of one variable, integral of, 33 
functions 
change of variables related to, 118 
continuity of, 21-22 


parameterization for graph of 
spherical equation, 61 

parameterizations, 66—67 

parameterized surfaces, 59-60 

particle speed moving along 
helix, 96 

rate of change for functions, 76 

rectangular coordinates, 54 
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defining integrals of, 102 

differentiating for variables, 
75-718 

evaluating at indicated 
points, 2 

evaluating limits of, 20 

examining matrices for, 79 

examples of, 1-2 

graphing, 4—5 

integrating, 36 

integrating over parameterized 
surfaces, 113 

integrating via Green's Theorem, 
158-159 

limits of, 19 

mental model of, 1 

picturing as graphs, 2 

plotting for variables, 2 

producing from vector fields, 
128-129 

of three variables, 11—12 

functions of multiple variables. See 
also variables 
composition of, 29-30 
limit for, 18 
Fundamental Theorem of Calculus, 
145, 147 
deducing, 149 
significance of, 145, 147 


G 


Gauss’ Divergence Theorem 

applying, 166-171 

example of, 171—172 
Generalized Stokes’ Theorem, 145 
geometric interpretation 

of curl, 164—166 

of divergence, 171—173 
grad f vector field, producing, 127 
gradient, finding for saddles, 80-81 


Advanced Calculus Demystified 


gradient of function, example 
of, 76-77 
graphing 
functions, 4—5 
level curves, 6—9 
graphs 
formula for surface area of, 108 
picturing functions as, 2 
sketching, 10 
gravity, overcoming, 147—148 
Green's Theorem 
versus Gauss' Divergence 
Theorem, 167 
geometric interpretation produced 
by, 171 
on Rectangular Domains, 
149-155 
relationship to Stokes’ 
Theorem, 160 
versus Stokes’ Theorem, 165 


I 


i vector, defining, 92 
integrals 
cancellation relative to Green’s 
Theorem, 156—157 
computing relative to Gauss’ 
Theorem, 170 
computing surface area of sphere 
with, 109 
computing surface of revolution 
with, 111 
computing volume with triple 
integrals, 44-47 
of curl of W over D, 164 
defining for functions, 102 
determining limits of integration 
for, 41 
evaluating, 121 


INDEX 


over nonrectangular domains, 
38-44 
setting up, 43 
for Stokes’ Theorem, 162 
integrand, determining for volume, 117 
integration 
act of, 33 
over nonrectangular domains, 
38-44 
performing, 36 
relationship to rectangles, 33 
integration order, determining for triple 
integral, 46 


J 


j vector, defining, 92 


K 
k vector, defining, 92 


L 


Lagrange multipliers method, using, 87 
Law of Cosines, explanation of, 73 
level curves 
finding direction of, 78 
graphing, 6-9 
level sets, plotting, 11 
limit properties, relationship to 
integration and rectangles, 35 
limits 
for curve, 101 
calculation for continuous 
functions, 21 
definition of, 18 
evaluating for functions, 20 
of functions, 19 
proving non-existence of, 19—20 
significance of, 17 


269 


surface area as, 105 
volume as, 116 

limits of integration 
determining for outer 

integral, 41 

finding, 45-47 
significance of, 40 

line integrals 
defining, 138 
determining, 103 
interpreting, 133, 135 
physical motivation for, 134 
recalculating, 136 
relationship to arc length, 113 
relationship to vector fields, 

137-138 

representing, 135 
significance of, 146 

local maximum and minimum, 

determining, 81 


M 


matrix 
tracking information in, 79 
writing determinant of, 90 
maxima 
detecting via derivative test, 
79-80 
distinguishing from saddles, 80 
determining local maximum, 81 
finding, 84 
minima 
detecting via derivative test, 
79-80 
distinguishing from saddles, 80 
seeking for f, 87 
determining local minimum, 81 
multiple integral, definition of, 
150-151 
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0 


orientation 

calculating for rectangles relative 
to Green’s Theorem, 158 

checking relative to Gauss’ 
Theorem, 168 

determining for surface S, 161 

example of, 136 

producing, 140-141 

relevance to surface integrals, 140 

significance of, 135 


P 


parabola, equation for, 52 
paraboloid 
definition of, 52 
example of, 10, 59-60 
parallelepiped, explanation of, 89 
parallelograms 
computing area of, 88, 91, 92, 105 
computing signed area of, 89 
parameterizations 
of areas and volumes, 65-67 
computing partials of, 116 
computing volumes with, 118 
determining shapes from, 63-64 
differentiating, 95-96 
finding for graphs of 
equations, 61 
restricting domains of, 62 
using with volume, 116 
parameterized curves 
composition with, 26—29 
examples of, 12-15 
parameterized surfaces, definition 
of, 59 
partial derivatives. See also derivatives 
applying to vector fields, 129 
computing for parameterization, 
116 


Advanced Calculus Demystified 


computing relative to chain 
rule, 28 
example of, 25-26, 31, 75 
for parameterizations, 114 
representing with respect to x, 
24—25 
for surface area of sphere, 
106, 108 
for surface of revolution, 110 
particle speed, determining along 
helix, 96 
path independence of line integrals of 
gradient fields, 147—148 
perpendicular vectors, detecting, 74 
points 
locating in rectangular 
coordinates, 2-4, 50 
locating with spherical 
coordinates, 53 
plotting for functions, 2 
plotting with three coordinates, 
2-4 
plugging into functions, 2 
versus vectors, 69 
polar coordinates 
adapting to three dimensions, 49 
translating to rectangular 
coordinates, 66 
problems 
area of parallelogram spanned by 
vectors, 93 
area of parameterized plane, 67 
calculating curve length, 102 
circles for level curves, 9 
computing integral of f(x, y) 
over parameterized 
curve, 104 
continuous function, 22 
converting to rectangular 
coordinates, 55 
coordinate system, 15 


INDEX 


cosine of angle between 
vectors, 75 

curl of vector fields, 130—131 

curve parameterization, 15 

derivatives, 26 

determinant, 91 

divergence of vector fields, 
129, 131 

domain of functions, 22 

dot product for perpendicular 
vectors, 75 

dot product of vectors, 75 

function evaluation at point, 2 

Gauss’ Divergence Theorem, 
170-171 

geometric interpretation of 
curl, 166 

graphs of equations, 53 

graphs of functions, 10 

graphs of spherical 
equations, 57 

Green’s Theorem, 155 

Green’s Theorem over more 
general domains, 159—160 

integrals of functions, 
122-123 
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parameterizing surface with level 
curves, 65 

partial derivatives, 30 

plotting points on axes, 4 

right handed coordinate 
systems, 4 

second partial derivatives, 31 

shape of cylindrical equations, 
52-53 

shapes with parameterizations, 63 

showing parameterized 
curve, 102 

signed area of parallelepiped, 91 

signed area of parallelogram, 90 

for Stokes' Theorem, 163 

surface integral, 115 

triple integrals, 47 

unit vector perpendicular to 
vectors, 93 

vector fields, 127 

vector fields with curls, 131 

vectors, 72 

volume calculation, 118 

volumes of parameterizations, 67 

writing rectangular coordinates of 
points, 50 


integration, 37-38 

integration and integrals, 44 

intersections of graphs of 
functions with coordinate 


problems, optimizing, 83-84 
Pythagorean Theorem, applying to 
magnitude of vectors, 71 


planes, 5—6 
level curves, 9 R 
line integrals, 138 r and z, labeling axes as, 52 
local maxima, minima, and rectangles 
saddles, 82 applying Green's Theorem to, 156 


relationship to integration, 33 
rectangular coordinates 
converting from cylindrical 


parameterization for volume 
below cone, 67 
parameterizations, 15, 65 
parameterizations for graphs of coordinates to, 50 
equations, 61 converting spherical coordinates 
parameterized sphere, 63 to, 53 
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explanation of, 49 
locating points in, 50 
translating cylindrical coordinates 
to, 66 
translating polar coordinates 
to, 66 
region 
computing area of, 45, 119 
cutting into slabs, 38-44 
finding volume of, 33-34 
parameterizing relative to change 
of variables, 119 
revolution, determining surface of, 
109-111 
right hand rule, relationship to labeling 
axes, 3 
right handed coordinate system, 
using, 3 
river analogy 
applying to line integrals, 134 
applying to surface 
integrals, 139 


S 
saddles 
definition of, 10 
determining, 81 
distinguishing from minima and 
maxima, 80 
example of, 11, 79 
finding gradients for, 80-81 
scalar multiplication, relationship to 
vectors, 71 
second derivatives, shorthand notations 
for, 31. See also derivatives 
shapes 
defining with rectangular 
coordinates, 49 
describing with spherical 
coordinates, 53—55 


Advanced Calculus Demystified 


determining from 
parameterizations, 63—64 
with parameterizations, 63 
sphere of radius 1, equation of, 12 
spheres 
bounding volumes by, 116 
computing surface area of, 
111-112 
determining surface area of, 
106-108 
parameterizing, 106, 116 
spherical coordinates 
converting to rectangular 
coordinates, 53 
describing shapes with, 
53-55 
equation of sphere in, 62 
locating points with, 53 
parameterizing unit sphere by 
means of, 168 
spherical equations, graphing, 55—57 
spiral, parameterizing, 14 
Stokes' Theorem 
applying, 160—163 
geometric interpretation produced 
by, 171 
versus Green's Theorem, 165 
producing geometric 
interpretation of curl with, 
164—166 
stream analogy 
applying to line integrals, 134 
applying to surface integrals, 139 
summation properties, relationship to 
integration and rectangles, 35 
surface area. See also area 
computing, 104—106 
computing for sphere, 111—112 
computing with double 
integral, 109 
determining for sphere, 106-111 


INDEX 


as limit, 105 
relationship to surface 
integral, 113 
surface integrals 
calculating, 114 
defining for vector fields, 
139-142 
definition of, 113 
relationship to surface area, 113 
surface of revolution, determining, 
109-111 
surfaces 
determining boundary of, 161 
differentiating parameterizations 
of, 96 
parameterizing, 97 
parameterizing with cylindrical 
coordinates, 114 


T 


tangent plane, 78 
theorems 
Fubunr's Theorem, 36 
Gauss' Divergence Theorem, 
166-171 
Green's Theorem, 171 
Green's Theorem on Rectangular 
Domains, 149-155 
Green's Theorem over more 
general domains, 156—160 
path independence of line 
integrals of gradient fields, 
147-148 
Stokes' Theorem, 160-166, 171 
triangle, applying Law of Cosines 
to, 73 
trigonometry, applying to dot 
product, 73 
triple integrals, computing volume 
with, 44—47 
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U 
unit sphere, parameterizing relative to 
Gauss' Theorem, 168 
unit vector 
definition of, 71 
perpendicular example of, 93 
pointing in same direction, 93 


V 


vector fields 

applying partial derivatives 
to, 129 

defining surface integrals for, 
139-142 

definition of, 125-126 

differentiating, 166—167 

divergence of, 129 

integrating divergence of, 
166-171 

producing functions from, 
128-129 

relationship to line integrals, 
137-138 

vectors 

adding, 69—70 

computing cross product for, 92 

computing magnitude of, 96 

defining i, j, and k, 92 

detecting perpendicular 
vectors, 74 

defining product of, 72 

finding cosine of angles 
between, 74 

magnitude for surface of 
revolution, 110 

magnitude of, 71, 98, 107—108 

multiplying by numbers, 71 

versus points, 69 

producing, 127 


274 Advanced Calculus Demystified 


producing from cross 
product, 91 

product of, 72 

subtracting, 71 

volumes 

bounding by spheres, 116 

computing for integration, 40—42 

computing with 
parameterizations, 118 

computing with triple integrals, 
44-47 

determining, 117 

above nonrectangular areas, 
38-44 


under graph and above 
rectangle, 36 
parameterizing, 65-67 
positive state of, 117 
of thin slab, 39 
volumes of regions, finding, 33-35 


X 
xy-plane 
example of, 5 
finding area of region R of, 45 
relationship to integration, 33 
xz-plane, finding intersection with, 7 


